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CONVERSION  FACTORS,  LI.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  tins  report  can  be  converted  to  metric  (SI) 
units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  y ards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton, long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

'To  obtain  Celsius  (C)  tem|>erature  readings  from  Fahrenheit  (F)  readings,  use  formula:  C - (5/9)  (F  — 32). 
To  obtain  Kelvin  (K)  readings,  use  forutnla.  K ' (5/9)  (F  - 32)  +273.15. 
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SYMBOLS  AND  DEFINITIONS 


A1  >a2 

al 

B 

c0 

CT 

Fj(t) 

KHij 

KMi  j 

kl,k2 

L 

ij 

U,y,t) 


6\,&2 

n(x,t) 

ni(x,t) 

nx(x,t) 

Xij 

Hj 


P 


O) 

(01,  0)2 


Amplitudes  of  two  incident  waves 

Amplitude  of  sway,  heave,  or  roll  motion  for  i = 1,2,3 
Characteristic  beam  of  breakwater 
Body  contour 

Transmission  coefficient 

Sway,  heave,  or  roll  exciting  forces  or  moment  for  j = 1,2,3 

Hydrostatic  restoring- force  coefficient  for  force  in  the  jth 
direction  due  to  motion  in  the  ith  direction 

Similar  to  KHij  but  due  to  the  mooring  system 

Wave  numbers  of  two  incident  waves 

Incident  wavelength 

Mass  or  moment  of  inertial  when  i = j,  0 when  i ^ j 

Unit  interior  normal  to  body  surface 

Pressure 

Vector  from  center  of  gravity  to  a point  on  the  body  surface 
Sway,  heave,  or  roll  motion;  speed  or  acceleration 
Phase  angle 

Phase  angles  for  two  incident  waves 

Free-surface  elevation 

Wave  surface  elevation  for  incident  wave 

Wave  surface  elevation  for  transmitted  wave 

Damping  coefficient  for  force  in  the  jth  direction  related  to 
velocity  in  the  ith  direction 

Added-mass  or  inertial-force  coefficient  for  force  in  the 
jth  direction  related  to  acceleration  in  the  ith  direction 

Fluid  density 

Velocity  potential 

Frequency 

Frequencies  for  two  incident  waves 


FLOATING  BREAKWATER  FIELD  ASSESSMENT  PROGRAM, 
FRIDAY  HARBOR,  WASHINGTON 

by 

B.H.  Adee , E.P.  Richey, 
and  D.R . Christensen 


I . INTRODUCTION 

Floating  structures  for  use  in  the  attenuation  of  water  waves  were 
introduced  by  Joly  (1905).  Little  was  done  with  the  concept  until  the 
Bombardon  floating  breakwater  was  deployed  to  form  a harbor  during  the 
Normandy  invasion  of  World  War  II.  The  use  of  mobile  harbors  for  po- 
tential military  applications  provided  the  incentive  for  extensive  work 
during  the  postwar  years.  Representative  articles  from  this  period 
include  those  by  Minikin  (1948)  who  discussed  floating  breakwaters  in 
general  terms,  Carr  (1951)  who  used  basic  mechanics  to  predict  trans- 
missivn  characteristics,  and  the  review  of  the  performance  of  the  Bom- 
Bardon  by  Lochner,  Faber,  and  Penny  (1948).  In  1957,  the  Naval  Civil  En- 
gineering Laboratory,  Port  Hueneme,  California,  began  a concerted  ex- 
ploration of  the  existing  knowledge  of  transportable  units  that  could 
serve  as  breakwaters  or  piers.  Results  of  the  study  are  summarized  in 
Naval  Civil  Engineering  Laboratory  (1961),  which  was  an  invaluable 
state-of-the-art  assessment  with  particular  emphasis  on  military  uses 
under  the  rather  severe  site  criteria  of  an  incident  wave  with  a 15-foot 
height,  13-second  period,  minimum  water  depth  of  40  feet,  inshore  trans- 
mitted wave  height  of  4 feet,  and  tidal  range  of  12  feet.  A sequel  to 
the  earlier  study  (Naval  Civil  Engineering  Laboratory,  1971)  surveyed 
concepts  for  "transportable"  breakwaters,  including  over  60  in  the 
"floating"  category.  Although  no  breakwater  system  was  disclosed  which 
would  meet  the  stringent  military  site  criteria  and  transportability 
requirement,  these  state-of-the-art  reviews  sparked  renewed  interest  in 
the  floating  breakwater  for  nonmilitary  applications.  A review  of  de- 
velopments in  floating  breakwaters  was  summarized  by  Richey  and  Nece 
(1974);  Seymour  (1974)  introduced  a new  and  innovative  concept  for  wave 
attenuation  using  a system  of  tethered  floats  which  may  have  application 
over  a wide  range  of  wave  conditions. 

Continually  increasing  pleasure  boat  ownership  has  nearly  exhausted 
the  available  supply  of  moorage  space  in  many  areas.  The  need  for  addi- 
tional moorage  space  in  conjunction  with  escalating  construction  costs 
and  more  stringent  environmental  restrictions  require  careful  scrutiny 
of  alternatives  to  the  traditional  fixed  breakwater  and  excavation  tech- 
niques employed  in  marina  construction.  Productive  time  in  weather- 
dependent,  waterborne  activities  such  as  construction,  logging,  and  cargo 
handling  could  be  increased  if  protective  floating,  transportable  break- 
waters were  used.  Other  uses  in  the  control  of  shoreline  erosion  and 
in  the  emerging  mariculture  industry  may  also  be  found. 
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The  information  on  the  performance  of  floating  breakwaters,  i.e., 
their  wave  attenuating  characteristics,  mooring  line  forces,  and  motions, 
is  contained  primarily  in  reports  of  laboratory  scale- model  tests  with 
monochromatic  incident  waves;  the  few  exceptions  are  the  early  analytical 
work  by  Carr  (1951)  and  the  occasional  piece  of  information  from  a full- 
scale  test  like  that  performed  by  Harris  (1974) . There  is  a need  for  a 
fundamental  analytical  procedure  to  predict  the  performance  characteris- 
tics of  floating  breakwaters  with  arbitrary  cross  section  when  exposed 
to  a given  incident  wave.  This  procedure  could  be  used  to  systematically 
compare  performance  information  available  in  the  literature,  to  examine 
new  design  proposals,  and  either  eliminate  or  reduce  and  systematize 
auxiliary  experimental  studies. 

The  development  of  the  predictive  procedure  was  the  primary  thrust 
of  the  project  with  the  concommitant  field  assessment  of  a full-scale 
floating  breakwater  in  operation  at  Friday  Harbor,  Washington  (Fig.  1). 
The  analytical  model  developed  from  the  two-dimensional,  linearized  so- 
lutions of  the  hydrodynamical  equations  formulated  in  terms  of  a boun- 
dary value  problem  for  the  velocity  potential.  The  model  was  refined 
progressively  by  comparisons  with  results  already  reported  in  the  li- 
terature, by  auxiliary  laboratory  tests,  and  by  the  results  from  the  Fri- 
day Harbor  field  program,  where  measurements  of  incident  and  transmitted 
waves,  mooring  line  forces,  and  acceleration  in  sway,  heave,  and  roll 
were  measured  over  a 6-month  period. 


Figure  1.  Aerial  view  of  Friday  Harbor  breakwater. 


II.  THEORETICAL  ANALYSIS 
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In  the  analysis  of  complex  systems  such  as  floating  breakwaters, 
there  is  a great  need  for  model-scale  experiments  to  predict  their  per- 
formance and  provide  data  for  the  application  of  rational  engineering 
design  principles.  Full-scale  measurements  are  also  extremely  valuable 
in  verifying  scaling  relationships  and  in  providing  confidence  that  the 
data  obtained  from  smaller  scale  experiments  are  reasonable. 

When  one  considers  the  myriad  possible  breakwater  configurations 
which  have  been  proposed  to  date  and  the  different  conditions  which  pre- 
vail at  each  potential  breakwater  site,  the  number  of  required  model 
tests  and  the  attendant  expense  are  very  large.  To  avoid  this  expense 
and  also  to  permit  parametric  studies  aimed  at  obtaining  optimum  break- 
water configurations,  a theoretical  model  was  developed.  The  goal  was 
to  theoretically  predict  the  performance  which  could  be  measured  in  la- 
boratory studies  or  at  prototype  installations. 

The  initial  restriction  imposed  on  the  theoretical  model  was  to 
consider  only  two-dimensional  conditions.  Under  this  restriction  the 
breakwater  is  assumed  to  be  very  long  in  one  direction  with  long-crested 
waves  approaching  so  that  their  crests  are  parallel  to  the  long  axis  of 
the  breakwater.  At  most  breakwaters  where  the  wave  climate  results  from 
wind- generated  waves,  this  condition  would  rarely  be  approached.  How- 
ever, experiments  performed  using  a boat  wake  to  generate  incident  waves 
on  the  beam  and  at  ar.  angle  to  a breakwater  indicate  larger  breakwater 
motions  and  larger  transmitted  waves  when  the  incident  wave  crests 
approach  parallel  to  the  long  axis  of  the  breakwater  (Stramandi,  1975). 

As  a design  tool,  a two-dimensional  theory  provides  information  on  the 
worst  conditions  which  might  be  expected  to  occur.  In  addition,  the 
extensive  two-dimensional  wave-channel  experiments  provide  the  data  need- 
ed to  test  the  theoretical  model . 

Throughout  the  development  of  the  theoretical  model,  every  attempt 
was  made  to  orient  the  model  toward  providing  a useful  tool  applicable 
to  realistic  problems.  To  perform  the  calculations  the  user  need  only 
know  the  incident  wave  frequencies  of  interest,  the  contour  of  the 
breakwater  cross  section  (catamaran-  or  trimaran-type  cross  sections  are 
permitted) , and  the  physical  properties  of  the  breakwater  (these  include 
mass,  mass  moment  of  inertia,  and  the  static  restoring- force  coeffi- 
cients) . 

The  approach  used  here  has  been  to  employ  the  techniques  which  naval 
architects  have  developed  to  deal  with  ship  motion  problems.  Mathema- 
tically, the  hydrodynamic  equations  are  formulated  in  terms  of  a boun- 
dary value  problem  for  the  velocity  potential.  Solution  of  this  complete 
problem  is  presently  impossible  because  the  free-surface  boundary  condi- 
tion is  nonlinear.  An  approximate  solution  may  be  obtained  if  restric- 
tions are  imposed  on  the  boundary  value  problem,  and  the  procedure  of 
linearization  is  applied.  The  restrictions  limit  the  applicability  of 
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the  solution  to  cases  of  small  incident  wave  amplitude  and  small  motion 
response  of  the  breakwater. 


When  using  the  linearized  theory  which  is  presented  here,  one  must 
be  well  aware  of  the  limits  of  applicability  which  are  imposed  on  the 
results  in  order  to  permit  the  formulation  of  a tractable  mathematical 
problem.  Care  must  also  be  exercised  because  these  restrictions  may 
exclude  phenomena  which  occur  in  nature  from  appearing  in  the  mathemati- 
cal analysis.  For  instance,  field  observations  clearly  demonstrate 
the  occurrence  of  mooring  line  force  oscillations  at  periods  greater  than 
those  which  could  be  attributed  directly  to  wind- generated  wave  exci- 
tation. Using  a linearized  approach,  these  long-period  oscillations 
would  not  appear  in  the  analysis.  A theoretical  model  which  includes 
nonlinear  behavior  of  the  system  is  required  if  these  long-period  os- 
cillations are  to  be  included. 

A possible  nonlinear  mechanism  for  the  transfer  of  wave  energy  to 
lower  frequencies  has  been  postulated  and  is  presented  to  supplement  the 
linear  analysis. 


1 . Linear  Theoretical  Model . 

The  problems  involved  in  theoretically  predicting  the  performance  of 
a two-dimensional  floating  breakwater  are  illustrated  in  Figure  2.  Here 
an  incident  wave  approaches  the  breakwater  on  the  beam.  A part  of  the 
energy  contained  in  the  incident  wave  is  reflected,  part  passes  beneath 
the  breakwater,  and  some  is  lost  through  dissipation.  Another  part  of 
the  incident  wave  energy  excites  the  motions  of  the  breakwater.  These 
motions  are  restrained  by  the  mooring  system.  The  oscillating  break- 
water in  turn  generates  waves  which  travel  away  from  the  breakwater  in 
the  directions  of  the  reflected  and  transmitted  waves.  The  total  trans- 
mitted wave  is  the  sum  of  the  component  which  passes  beneath  the  break- 
water and  the  components  generated  by  the  breakwater  motions.  The  total 
reflected  wave  is  composed  similarly. 

In  completing  the  calculations,  the  information  which  is  of  most 
interest  to  the  designer  includes: 

(a)  Total  transmitted  and  reflected  waves  including  their 
components . 

(b)  Wave  forces  on  the  breakwater. 

(c)  Motions  of  the  breakwater. 

(d)  Forces  on  the  mooring  lines. 

For  the  two-dimensional  breakwater,  definitions  for  the  motions  are 
shown  in  Figure  2.  Sway  is  defined  as  the  oscillation  perpendicular  to 
the  long  axis,  or  along  the  x-coordinate  axis.  Heave  is  the  vertical 


dimensional  floating  breakwater. 
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motion  of  the  breakwater  along  the  /-coordinate  axis,  and  roll  is  the 
rotation  about  the  long  axis  or  the  z-coordinate  direction. 

As  long  as  the  problem  is  linear,  computing  the  performance  of  a 
floating  breakwater  may  be  separated  into  three  parts: 

(a)  Formulate  equations  of  motion, 

Calculate  hydrostatic  forces  and  moments. 

Evaluate  hydrodynamic  coefficients  in  equations  of  motion. 

Compute  exciting  forces  on  breakwater. 

Solve  for  the  motions  and  motion-generated  waves. 

Compute  forces  in  the  mooring  lines. 

(b)  Solve  for  the  waves  diffracted  by  a rigidly  restrained 
breakwater. 

(c)  Sum  components  to  obtain  total  reflected  and  total  trans- 
mitted waves. 

When  combined,  these  parts  of  the  calculation  provide  complete 
performance  data  for  a two-dimensional  breakwater. 


a.  Breakwater  Motions.  In  deriving  the  equations  of  motion, 
Newton's  law  is  used. 

mij  ai  = E forces;  (1) 

here : 

= motion  of  the  breakwater  in  sway,  heave,  and  roll  for 

i = 1,2,3,  respectively.  The  dot  above  indicates  differen- 
tiation with  respect  to  time. 


m^  = mass  or  mass  moment  of  inertia  when  i = j and  zero  when  i ^ j . 

Expanding  this  equation  to  include  the  various  forces  in  the  summa- 
tion yields: 


m^  = Fj  (inertial)  + F^  (wave  damping) 

♦ Fj  (friction)  + F^  (hydrostatic)  ♦ F^  (mooring) 

♦ Fj  (wave  exciting) 
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The  inertial  force  (or  added-mass  force)  arises  when  the  breakwater  acce 
lerates,  which  also  accelerates  the  fluid  around  it.  The  motion-gene- 
rated waves  are  moving  away  from  the  breakwater  and  result  in  the  wave- 
damping term.  A term  representing  the  forces  due  to  viscosity  is  in- 
cluded, but  these  forces  are  neglected  in  the  analysis.  Experience  in 
ship  motion  analysis  (Salvesen,  1970)  has  shown  this  to  be  acceptable 
for  all  motions  but  roll,  where  damping  may  make  a more  significant  con- 
tribution than  for  sway  and  heave  motions.  At  present,  the  main  reason 
for  neglecting  the  frictional  forces  is  that  they  lead  to  nonlinear 
terms  in  the  equations  of  motion,  which  make  their  solution  far  more 
complex.  Hydrostatic  forces  arise  because  of  changes  in  the  displaced 
volume  of  the  breakwater  when  it  moves.  In  this  analysis  the  mooring 
forces  are  modeled  as  simple  springs  with  their  contribution  to  the 
damping  and  inertial  forces  considered  small  in  comparison  to  similar 
terms  resulting  from  the  breakwater  motion.  The  wave  exciting  force 
results  from  the  incident  waves  striking  the  breakwater. 

If  we  neglect  the  nonlinear  terms  and  assume  that  the  fluid  is  in- 
viscid,  then  the  equations  of  motion  describing  the  coupled  sway,  heave, 
and  roll  motions  of  the  breakwater  are  of  the  form: 


1 {(mij  + »*ij>  *i  + Xij  ®i  + (KHij  + °i}  = Fj  (t)  (2> 

1=1  for  j = 1,2,3. 

The  symbols  are  defined  as  follows: 

p. . = added-mass  coefficient  with  the  representing  the 

^ added-mass  force  or  moment  in  the  jth  direction  due  to 
acceleration  in  the  ith  direction. 

X. ■ = damping- force  coefficient  relating  damping  force  or  moment 

^ in  the  jth  direction  to  velocity  in  the  ith  direction. 

KH,  • = hydrostatic  spring  constant  relating  the  restoring  force  or 
J moment  in  the  jth  direction  to  displacement  in  the  ith 
direction. 

KKTj  = similar  to  KHjj  but  due  to  the  mooring  system. 


Fj  = exciting  force  or  moment  in  the  jth  direction. 

In  order  to  solve  these  equations,  the  physical  mass  and  moment  of  iner- 
tia, added  mass  and  damping  coefficients,  static  spring  constants,  and 
the  exciting  forces  must  all  be  known.  Mass  and  moment  of  inertia  are 
computed  directly  from  the  specifications  of  the  breakwater  section.  The 
KH  are  derived  directly  from  hydrostatic  considerations  in  Appendix  A, 
while  approximate  values  for  KMjj  are  obtained  by  using  a discretized 
approximation  for  the  mooring  line  as  described  in  Appendix  B.  Potential 


theory  and  the  principle  of  linear  superposition  permit 
the  hydrodynamic  coefficients  and  forcing  function  mjj  , 


derivations  for 
Xjj  and  Fj  (t) . 


Steady-state  solutions  of  the  form: 

ou(t)  = a^  sin  (iot  + 6^)  for  i = 1,2,3  (3) 

are  assumed.  Substitution  of  the  assumed  solution  (eq,  3)  into  the 
equations  of  motion  (eq.  2)  yields  a set  of  linear  algebraic  equations 
which  may  be  solved  for  the  unknown  amplitudes  and  phase  angles  aj,  and 
6^.  Transfer  functions,  , are  then  defined  by  the  a^  and  <5^  since 
the  incident  waves  are  assumed  to  be  sinusoidal. 

b.  Hydrodynamic  Coefficients  and  Waves.  Potential  theory  is  em- 
ployed in  computing  the  reflected  and  transmitted  waves,  hydro- 
dynamic  coefficients  and  the  exciting  forces.  Under  the  assumptions  of 
small  incident  waves,  small  breakwater  motions  and  an  inviscid  fluid, 
the  velocity  potentials  may  be  found  and  the  problem  subdivided  using 
the  principle  of  linear  superposition.  The  total  velocity  potential: 


^total  ^incident  + ^diffracted  + 


(i) 


motion 


for  i 


1,2,3  (4) 


is  the  sum  of  the  incident  wave  potential,  the  diffracted  wave  poten- 
tial and  the  potential  resulting  from  forced  sway,  heave,  and  roll  mo- 
tions. 

The  incident  wave  potential  is  well  known  and  may  be  expressed 
directly.  Obtaining  the  diffracted  wave  and  breakwater  motion  poten- 
tials requires  the  solution  of  boundary  value  problems.  These  problems 
and  their  solutions  are  described  in  Appendix  C.  Appendix  D provides 
the  computer  program  used  to  calculate  breakwater  performance. 

When  the  velocity  potentials  have  been  obtained,  the  free-surface 
elevation  at  any  position  is  found  using  the  linearized  free-surface 
boundary  condition: 

n(x,t)  = - j 4>t  (x.O.t) . (5) 

Here: 

n(x,t)  = free-surface  elevation  measured  from  Stillwater 
level  (y  = 0) , 

g = acceleration  of  gravity, 

$t(x,0, t)  = derivative  of  the  velocity  potential  with  respect  to 
time  evaluated  at  y = 0. 
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\otal  ^ = 


g incident  + diffracted 


♦ 


(i) 


t motion 


(x,0,t) } . 


(6) 


The  fluctuating  component  of  pressure  in  the  fluid  and  on  the  breakwater 
hull  surface  may  be  computed  using  Bernoulli's  equation: 


P(x,y,t)  = - P <f>t  (x,y,t) . 


(7) 


By  computing  pressures  on  the  hull  surface  and  integrating  these  around 
the  contour,  the  forces  on  the  breakwater  may  be  computed.  The  force 
per  unit  length  acting  on  the  breakwater  is  then: 


F (t) 


P n ds. 


(8) 


In  this  case. 


F(t)  = 
-► 

n = 

Co  = 
The  rolling 

M(t)  = 


force  on  the  breakwater, 
unit  interior  normal  vector 
contour  of  breakwater  cross 
moment  is : 

f ->  - 

P r x n ds. 


on  the  hull  surface, 
section. 


where. 


(9) 


r = the  vector  from  the  center  of  gravity  to  a point  on  the  surface 

To  compute  the  exciting  forces  on  the  breakwater  in  linear  theory,  the 
pressure  due  to  the  incident  and  diffracted  waves  is  integrated  over  the 
hull  surface.  These  forces  and  moments  become: 


r 


Fj(t)  - { - P 

Co 

t*t 

incident 

(s,t) 

F2(t)  = { - P 

4 

r 

incident 

(s,t) 

Vi 

0 

Fj(t)  - { - p 

C 

incident 

(s,t) 

o 


+ *t  diffracted(s't)]  " ds}  * T 

* *t  diffracted n ds } J 

+ \ diffracted(s*t)]?X"ds}  ‘ k 

(10) 


Hydrodynamic  coefficients  are  found  using  the  potential  resulting 
from  forced  oscillation  of  the  breakwater.  In  this  case  the  pressure 
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integrated  over  the  surface  has  a component  in  phase  with  acceleration 
and  a component  in  phase  with  velocity.  The  component  in  phase  with 
acceleration  is  normally  referred  to  as  the  added  mass,  while  the  compo- 
nent in  phase  with  velocity  is  the  damping. 

The  hydrodynamic  coefficients  shown  in  this  section  are  derived  in 
greater  detail  in  Appendix  C. 

c.  Mooring  Forces.  At  the  time  the  spring  constants  for  the  mooring 
lines  are  computed,  mooring  force  coefficients  are  also  calculated.  These 
are : 

— = change  in  mooring  line  force  per  unit  displacement  in 
Aai  sway,  heave,  or  roll  when  i = 1,  2,  or  3,  respectively. 

The  forces  in  the  mooring  lines  may  then  be  computed  once  the  motions 
have  been  found. 

3 

Mooring  Force  = £ (A — ) a'(t) 

i=l  Aai  1 

The  description  of  the  linear  system  is  now  complete.  The  block  diagram 
in  Tigure  3 shows  the  relationships  among  the  calculations  which  are 
required . 


Nonlinear  Theoretical  Model. 


Measurements  taken  at  the  Tenakee,  Alaska,  floating  breakwater  be- 
fore this  research  program  was  begun  indicated  the  presence  of  a long- 
period  oscillatory  motion  of  the  breakwater.  These  long-period  motions 
were  manifested  most  clearly  in  the  measured  mooring  line  forces.  Look- 
ing at  these,  one  can  visually  observe  an  oscillation  with  a period  of 
about  60  seconds  superimposed  over  the  expected  shorter  period  oscilla- 
tions. Figure  4 shows  the  results  of  a spectral  analysis  of  the  seaward 
mooring  line  data  after  a low-pass  filter  has  been  applied  (the  tech- 
nique for  performing  the  spectral  analysis  is  given  in  Section  III  of 
this  report) . 

The  linear  theoretical  model  permits  the  system  to  respond  only  at 
the  frequency  of  the  incident  wave.  In  order  to  explain  the  presence  of 
these  long-period  oscillations,  nonlinearities  must  be  included  in  the 
analysis.  To  perform  a mathematically  complete  analysis  including  all 
nonlinear  effects  is  beyond  the  present  state  of  the  art.  However,  in 
the  case  of  the  floating  breakwater,  one  can  show  that  if  two  incident 
waves  are  considered  and  second-order  terms  are  retained,  then  an  excit- 
ing force  is  present  at  the  difference  between  the  frequencies  of  the 
incident  waves.  The  complete  derivation  in  Appendix  E shows  that  the 
nonlinear  pressure  may  be  expressed  as: 
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(k1+k2)y 

- 2u)jU>2AjA2e  cos  [(k^  - k2)x 

- (wj  - o>2)t  + - 62]}, 

where. 


= fluid  density. 


u) ^ , uj 2 = incident  wave  frequencies, 

A1,A2  = incident  wave  amplitudes,  2 2 


k ,k.  = incident  wave  numbers  = , 

* • 8 


6^,62  = incident  wave  phase  angles. 
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Combining  this  pressure  with  the  pressure  obtained  from  the  linear  theory 
and  integrating  over  the  hull  would  provide  additional  exciting-force 
terms  at  zero  frequency  and  at  the  difference  frequency.  Carrying  the 
nonlinear  exciting-force  terms  back  through  the  linear  response  analysis 
should  provide  a quasi-linear  approach.  While  there  is  no  reason  to 
expect  this  to  provide  exact  correlation  with  measured  data,  the  quasi- 
linear  approach  would  at  least  permit  the  natural  phenomena  to  enter  into 
the  mathematical  analysis. 


One  would  expect  terms  to  appear  in  the  second-order  pressure  (eq. 
11)  at  twice  the  incident  wave  frequency  and  at  the  sum  of  the  inci- 
dent wave  frequencies.  Terms  at  twice  each  of  the  incident  wave  fre- 
quencies can  be  derived  by  applying  the  trigonometric  relationships  to 
the  terms  at  zero  frequency.  While  a term  at  the  sum  of  the  incident 
wave  frequencies  does  not  appear  in  the  second-order  incident  wave  po- 
tential, this  term  may  result  when  the  second-order  potentials  repre- 
senting diffraction  or  forced  oscillation  in  calm  water  aie  included. 

A great  deal  more  effort  is  required  in  this  area  to  complete  the 
analysis.  There  is  also  one  other  area  where  a nonlinear,  or  quasi- 
linear,  analysis  should  be  investigated.  This  is  in  the  roll-damping 
coefficient.  Here,  viscous  effects  seem  to  be  important,  and  while  the 
problem  has  not  been  dealt  with  within  the  present  study,  investigators 
have  included  a term  proportional  to  velocity  squared  in  the  equation 
for  roll  motion. 


3 . Results. 

The  computer  program  given  in  Appendix  D has  been  developed  to 
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calculate  the  values  of  hydrodynamic  coefficients,  breakwater  motions, 
and  the  wave  field.  Input  variables  include: 

(a)  The  body  contour,  C , represented  by  a series  of  points  on 
the  contour. 

(b)  The  physical  properties  of  the  body:  mass,  mass  moment  of 
inertia,  and  position  of  the  center  of  gravity, 

(c)  The  mooring  system  spring  constants. 

(d)  The  hydrostatic  restoring  spring  constants. 

(e)  The  incident  wave  frequency,  u> . 

In  this  program  the  exciting  forces  and  moments  appearing  in  the  equa- 
tions of  motion  and  the  fixed-body  parts  of  the  transmitted  and  reflect- 
ed waves  are  found  by  computing  the  forces,  moments,  and  waves  which 
result  when  a rigidly  fixed  body  is  struck  by  a sinusoidal  incident  wave 
of  frequency  w.  Motions  are  found  by  computing  the  steady-state  solu- 
tion to  the  three  equations  of  motion.  The  hydrodynamic  coefficients 
and  the  waves  generated  by  the  body  motions  are  found  by  computing  the 
forces,  moments,  and  waves  which  result  when  the  body  is  forced  to  os- 
cillate in  Stillwater  in  pure  sway,  pure  heave,  or  pure  roll. 

The  physical  properties  used  in  the  performance  calculations  for 
the  various  breakwaters  are  collected  in  Appendix  F. 

a.  Wave  Transmission.  To  assess  the  performance  of  a floating 
breakwater,  one  quantity  which  is  commonly  used  is  the  transmission 
coefficient.  This  is  simply  the  transmitted  wave  amplitude  divided  by 
the  incident  wave  amplitude,  |ryr(x,t) |/ | n j Cx , t) | for  monochromatic  inci- 
dent waves . 

(1)  Proposed  Oak  Harbor  Breakwater.  At  one  time  the  Corps  of 
Engineers  was  considering  a marina  and  floating  breakwater  at  Oak  Har- 
bor, Washington.  Model  experiments  were  carried  out  by  Davidson  (1971) 
to  determine  transmission  characteristics  and  mooring  forces.  The  break- 
water itself  had  a catamaran- type  cross  section.  A comparison  between 
the  theoretically  predicted  and  experimentally  measured  transmission 
coefficient  is  shown  in  Figure  5.  This  figure  as  well  as  the  others 
plotted  in  this  section  and  Section  IV  were  drawn  using  a CALCOMP  plot- 
ter. The  plotting  program  uses  a parabolic  fit  to  determine  additional 
points  between  the  given  data.  Varying  numbers  of  data  points  were  used 
to  describe  each  curve  depending  on  its  behavior.  Data  points  were 
closely  spaced  in  regions  where  thq  theoretical  predictions  indicated 
large  changes  in  curvature.  Wavelength  is  calculated  in  all  the  figures 
using  the  relationship  between  wavelength  and  period  for  waves  in  deep 
water. 

In  this  case,  the  results  compare  reasonably  well  except  for  the 


□ EXPERIMENTAL  MEASUREMENTS.  DAVIDSON  (1971) 


Figure  5.  Transmission  coefficient  for  proposed  Oak  Harbor  breakwater 


predicted  dip  in  transmission  just  above  a B/L  (beam/wavelength)  of 
0.2.  There  is  also  some  difference  at  higher  B/L  ratios. 

The  theory  predicts  that  the  part  of  the  transmitted  wave  which 
would  result  where  the  body  is  rigidly  fixed  is  almost  1 for  a B/L  less 
than  0.1  and  drops  rapidly  at  higher  B/L  ratios  to  the  point  where  it  is 
of  little  consequence  beyond  0.2.  Waves  generated  by  the  breakwater  mo- 
tions play  an  increasing  role  for  B/L  ratios  above  0.15.  Heave  motion 
is  the  major  contributor  to  the  transmitted  wave  in  the  very  narrow  band 
of  B/L  between  0.15  and  0.18  with  a predicted  heave  resonance  at  a B/L 
of  about  0.18.  The  dip  occurs  because  the  waves  generated  by  heave  and 
sway  motions  are  almost  180°  out  of  phase  and  cancel  each  other  out.  At 
B/L  ratios  above  0.25,  sway  motion  assumes  an  increasingly  dominant  role. 
Roll  motions  are  small  throughout  and  generate  only  very  small  waves. 

(2)  Rectangular  Breakwater.  A breakwater  of  rectangular 
cross  section  with  the  same  beam  and  draft  as  the  proposed  Oak  Harbor 
breakwater  was  tested  at  the  University  of  Washington  by  Nece  and  Richey 
(1972).  Results  for  the  water  depth  of  29.5  feet  are  shown  in  Figure 
6. 


Again  the  agreement  is  reasonable.  Further  experiments  with  this 
model  have  confirmed  the  existence  of  the  trough  at  a B/L  of  0.2.  How- 
ever, this  phenomenon  can  be  observed  only  for  very  small  wave  heights. 
For  practical  purposes,  the  dip  may  be  smoothed  over  considerably.  The 
major  discrepancy  is  at  the  high  B/L  ratios  where  the  theory  shows  con- 
siderably greater  transmission  than  is  actually  measured  in  the  model 
tests.  Since  the  transmitted  wave  is  almost  totally  a result  of  sway 
motion,  the  problem  must  lie  in  the  wave  predicted  by  this  motion. 

Over  the  entire  range  of  wavelengths  of  interest,  the  predicted 
results  follow  the  pattern  previously  discussed  for  the  proposed  Oak 
Harbor  breakwater.  The  transmitted  wave  is  almost  completely  a result 
of  fixed-body  transmission  followed  by  regions  of  heave  resonance,  heave 
and  sway  cancellation,  and  finally,  sway  wave  generation  as  the  B/L 
increases . 

It  is  interesting  to  note  that  there  is  very  little  difference  be- 
tween the  open-well  breakwater  and  the  closed  rectangle  of  the  Same 
overall  dimensions. 

(3)  Rectangular  Breakwater  Tested  by  Sutko  and  Haden.  In 
some  recent  experiments  Sutko  and  Haden  (19?4)  have  examined  the  effect 
that  restricting  breakwater  motions  has  on  the  transmission  coefficient. 
They  used  a rectangular  breakwater  model  with  a beam-to-draft  ratio  of 
1.5.  Plexiglas  end  assemblies  were  used  to  restrict  the  breakwater 
motions . 

Figure  7 shows  the  transmission  coefficient  when  the  breakwater  is 
restricted  to  sway  motions  only.  Here,  the  transmitted  wave  contains  a 
component  resulting  from  the  fixed-body  transmission  and  a component 
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resulting  from  the  wave  generated  by  the  sway  motion.  At  very  low  B/L 
ratios  the  fixed-body  transmission  is  the  more  important  component.  At 
a B/L  of  about  0.1  both  the  fixed-body  transmission  and  the  sway-generat- 
ed wave  are  of  equal  importance.  At  B/L  ratios  higher  than  0.215,  the 
wave  generated  by  sway  motion  dominates.  The  agreement  between  theory 
and  experiment  is  quite  reasonable  for  this  case. 

A comparison  when  the  breakwater  is  restricted  to  heave  motion  only 
is  shown  in  Figure  8.  There  is  clearly  a discrepancy  in  this  case  be- 
tween measured  and  theoretically  predicted  transmission  coefficient  near 
the  B/L  ratio  of  0.13.  As  a matter  of  fact,  the  theory  predicts  a heave 
resonance  in  this  region  which  does  not  seem  to  be  supported  by  the  mea- 
sured data. 

In  examining  the  mechanism  used  to  restrain  the  breakwater  motion, 
it  seemed  possible  that  this  apparatus  was  introducing  damping  into  the 
system.  To  test  this  supposition,  transmission  coefficients  were  com- 
puted with  the  calculated  hydrodynamic  damping  increased  by  an  arbitrary 
amount.  The  major  effect  of  increasing  the  damping  was  to  decrease  the 
transmission  near  the  heave  resonance  region.  With  damping  at  three 
times  the  hydrodynamic  value,  the  results  were  quite  close  to  the  experi- 
mental measurements.  Increasing  the  damping  beyond  this  had  very  lit- 
tle additional  effect  on  the  predicted  transmission  coefficient.  The 
scatter  which  appears  in  the  experimental  data  in  this  region  is  a 
further  indication  that  some  nonrepeatable  effect  may  be  influencing 
the  experimental  results.  So  long  as  the  additional  damping  is  included 
in  the  theoretical  calculations,  the  results  compare  well  with  experi- 
mental measurements. 

Figure  9 shows  a comparison  between  model  measurements  when  the 
model  is  unrestrained  except  by  a horizontal  mooring  cable  and  the 
theoretically  predicted  results  without  mooring  restraints.  The  theo- 
retical results  are  characteristic  of  the  rectangular  breakwater  with 
the  dip  in  transmission  near  a B/L  equal  to  0.2.  The  pattern  of  inter- 
actions between  motion-generated  waves  and  fixed-body  transmission  is 
similar  to  the  previous  description.  The  agreement  between  these  re- 
sults indicates  that  the  theoretical  model  may  also  yield  the  correct 
results  when  the  model  is  free  to  heave  only.  At  least  further  experi- 
mental investigation  is  warranted. 

(4)  Alaska-Type  Breakwater.  The  State  of  Alaska  has  embarked 
on  an  ambit ious~program  for  constructing  moorages  using  floating  break- 
waters. As  part  of  a Sea  Grant  project  the  University  of  Washington  has 
been  studying  the  performance  of  this  type  breakwater.  A theoretically 
predicted  transmission  coefficient  and  the  transmission  coefficient  mea- 
sured in  model  tests  are  shown  in  Figure  10.  The  model  tests  were  con- 
ducted using  very  small  incident  waves  (wave  heights  on  the  order  of  0.2 
to  0.3  feet  at  prototype  scale).  Results  for  larger  wave  slopes  were 
not  included  in  the  figure  but  do  show  the  same  trends  with  lower  values 
of  transmission  coefficient.  Theoretical  predictions  without  added 
damping  and  with  double  the  hydrodynamic  damping  are  shown  in  Figure  10. 


□ EXPERIMENTAL  MEASUREMENTS.  SUTKO  ANO  HAOEN  (1974) 


Figure  8.  Transmission  coefficient  for  a rectangular  breakwater 
restricted  to  heave  motion  only. 
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Transmission  coefficient  for  a rectangular  breakwater. 


Transmission  coefficient  for  Alaska-type  breakwater 


Clearly  the  increased  damping  makes  a significant  difference  in  the 
results . 

Some  insight  into  the  performance  of  this  breakwater  may  be  gained 
by  following  the  theoretical  results  as  a function  of  B/L.  At  very  low 
B/L  the  fixed-body  transmission  dominates.  The  trough  at  a B/L  of  0.4 
comes  mainly  from  the  interaction  of  the  roll-generated  wave  and  the 
fixed-body  transmission.  For  the  next  peak  at  B/L  of  O.S  the  roll- 
generated wave  dominates  as  the  roll  resonant  frequency  is  encountered. 
The  next  trough  at  B/L  of  0.65  is  a result  of  all  three  motion-generated 
waves  interacting  with  the  fixed-body  transmission  making  only  a rela- 
tively small  contribution  to  the  transmitted  wave.  The  following  peak 
at  B/L  of  0.7  results  from  interactions  among  the  motion-generated  waves 
which  are  of  about  equal  magnitude.  At  a B/L  of  0.86  the  heave-generated 
wave  dominates  again,  but  as  B/L  increases  beyond  this  the  effect  of 
heave  and  roll  are  rapidly  decreasing  while  sway  motion  is  becoming  the 
dominant  wave-generating  mechanism.  In  the  region  of  B/L  between  0.4 
and  1.0,  changing  the  physical  properties  of  the  breakwater  can  have 
a marked  effect  in  shifting  the  peaks  and  troughs  by  altering  the  heave 
and  roll  resonant  frequencies. 

Experience  with  linear  ship  motion  theory  has  shown  that  the  worst 
agreement  between  predicted  and  measured  motions  occurs  when  rolling  mo- 
tions are  considered  (Salvesen,  1970).  This  discrepancy  is  often  over- 
come by  arbitrarily  increasing  the  computed  roll  damping  to  compensate 
for  the  viscous  damping  which  is  neglected.  As  indicated  in  Figure  10, 
when  damping  is  added  the  theory  gives  a better  prediction  where  roll 
motion  plays  a significant  role.  This  places  a significant  restriction 
on  the  theory  requiring  careful  monitoring  of  predicted  roll  motion. 

Where  the  theory  predicts  large  roll  motion,  additional  damping  will  be 
required  to  obtain  results  comparable  to  measurements. 

Figure  11  shews  the  predicted  fixed-body  transmission  coefficient 
and  the  results  of  model  tests.  Agreement  is  quite  close  except  at  B/L 
of  0.78.  The  peak  in  predicted  transmission  may  be  due  to  a resonance 
of  the  waves  within  the  well  of  the  catamaran  breakwater.  There  is 
another  peak  near  B/L  of  1.4  indicating  the  presence  of  higher  harmonic 
resonances  as  well.  Model  tests  show  at  least  a slight  hump  in  this 
region  suggesting  that  the  theoretical  prediction  clearly  overestimates 
the  effect  of  this  phenomena,  but  that  this  probably  is  occurring  in 
real  life. 

For  the  data  measured  in  the  field,  the  transmission  coefficient  is 
defined  as  the  square  root  of  the  transmitted  wave  spectral  density 
divided  by  the  incident  wave  spectral  density.  Figure  12  shows  the 
transmission  coefficient  derived  from  the  data  obtained  at  the  Tenakee,* 
Alaska  breakwater.  The  theoretically  predicted  transmission  coefficient 
with  the  computed  hydrodynamic  damping  doubled  is  also  shown  for  com- 
parison. Details  of  the  technique  used  in  the  spectral  analysis  of  the 
field  data  may  be  found  in  Section  III. 
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irans mission  coefficient  for  Alaska-type  breakwater 
Tenakee,  Alaska  (record  TK7-23) . ^ DreaKwater> 


It  should  be  noted  that  the  model  for  the  Alaska- type  breakwater 
was  not  built  to  the  correct  scale  to  represent  the  prototype.  Further 
investigation  of  the  physical  properties  of  the  prototype  after  the 
model  tests  were  complete  revealed  that  it  was  heavier  than  originally 
predicted.  The  physical  properties  used  in  making  the  theoretical  cal- 
culations are  correct  for  all  the  comparisons  made  in  this  report.  How- 
ever, care  must  be  exercised  in  comparing  the  model  test  results  and  the 
field  measurements  directly.  The  physical  properties  for  all  the  break- 
waters discussed  in  this  section  are  in  Appendix  F. 

The  first  trough  in  the  transmission  coefficient  curve  results  be- 
cause the  wave  generated  by  roll  tends  to  cancel  the  fixed-body  trans- 
mission. The  sway-generated  wave  is  small  but  cancels  a little  bit  of 
the  heave-generated  wave.  The  total  transmitted  wave  is  then  almost  in 
phase  with  the  heave-generated  wave  at  a slightly  reduced  amplitude. 
Complex  interactions  among  the  components  of  the  transmitted  wave  con- 
tinue to  result  in  oscillations  of  the  transmission  coefficient  up 
through  a B/L  of  0.9.  At  values  of  B/L  above  this,  the  transmitted  wave 
is  primarily  a result  of  sway  motion  except  for  the  peak  at  B/L  equal 
to  1.4  which  results  from  an  increase  in  the  fixed-body  transmission. 
Considering  the  complexity  of  the  breakwater  response,  the  agreement 
should  be  considered  to  be  reasonably  good. 

(5)  Friday  Harbor  Breakwater.  The  computed  transmission 
coefficient  for  the  Friday  harbor  breakwater  is  shown  in  Figure  13.  As 
in  the  case  of  the  Alaska  breakwater  calculations,  the  computations  of 
wave-damping  coefficients  have  been  arbitrarily  doubled  to  reduce  the 
excessive  calculated  motions  in  the  region  of  resonant  motions.  In  this 
figure  the  spacing  of  data  points  varies.  More  points  are  used  to  spe- 
cify the  curve  in  regions  of  rapid  change  so  that  the  plotted  result 
accurately  represents  the  theoretical  prediction. 

In  Figure  13,  the  first  trough  in  transmission  coefficient  at  about 
B/L  =0.5  results  from  heave-  and  roll -generated  waves  canceling  the 
fixed-body  wave  transmission.  This  transmission  coefficient  is  well  be- 
low the  transmission  coefficient  which  would  be  obtained  with  the  break- 
water rigidly  restrained  and  only  fixed-body  transmission  waves  passing 
through.  As  B/L  increases,  there  is  a peak  at  about  0.7.  At  this  point 
the  heave-generated  wave  has  almost  vanished,  and  the  fixed-body  trans- 
mission is  also  small.  The  larger  transmission  coefficient  is  primarily 
the  result  of  a roll-generated  wave  with  a smaller  component  resulting 
from  sway  motion.  The  next  trough  at  a B/L  of  0.9  occurs  as  the  heave 
motion-generated  wave  increases  and  cancels  the  roll  and  sway  motion- 
generated components.  The  fixed-body  transmission  is  very  small  at  B/L 
of  0.9.  As  B/L  increases  beyond  0.9  the  transmitted  wave  is  almost 
totally  the  result  of  sway  motion  of  the  breakwater. 

At  larger  B/L  ratios  there  are  several  oscillations  in  the 
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transmission  coefficient  curve.  In  this  region  one  must  be  careful  of 
the  analysis  because  there  are  certain  "irregular  frequencies"  or  "John" 
frequencies  where  the  approach  adopted  here  breaks  down  mathematically 
(John,  1950) . These  are  described  with  reference  to  the  integral  equa- 
tion technique  by  Frank  (1967) . It  is  extremely  difficult  to  predict 
where  the  first  of  these  irregular  frequencies  will  occur  when  the 
breakwater  cross  section  is  as  complicated  as  the  Friday  Harbor  break- 
water. If  this  cross  section  were  rectangular  with  the  same  exterior 
dimensions  as  the  Friday  Harbor  breakwater,  then  the  first  irregular 
frequency  would  occur  at  B/L  =1.7.  In  practice,  one  may  watch  for  this 
mathematical  phenomenon  by  checking  the  determinant  of  the  matrix  invert- 
ed to  solve  the  system  of  equations.  In  fact,  this  does  decrease  in  the 
region  of  B/L  of  1.7  but  does  not  indicate  a singular  matrix  for  the 
calculation  in  this  region  of  B/L.  Since  this  is  beyond  the  frequency 
range  of  primary  interest,  it  is  best  to  simply  view  the  results  at  B/L 
greater  than  1.7  with  extreme  caution.  The  oscillations  in  the  trans- 
mission coefficient  in  this  region  of  B/L  are  probably  the  result  of 
these  irregular  frequencies. 

b.  Breakwater  Motions.  In  the  wave  channel  experiments  perform- 
ed to  date,  there  has  been  no  attempt  to  compute  the  breakwater  motions. 
While  the  transmission  coefficient  is  the  primary  measure  of  breakwater 
performance,  the  motions  may  be  very  important  to  the  designer,  parti- 
cularly if  boats  are  to  be  tied  to  the  breakwater.  For  the  theoretical 
analysis,  this  is  a critical  intermediate  step  where  extensive  experi- 
mental measurements  used  for  comparison  would  be  invaluable. 

Friday  Harbor  Breakwater.  The  theoretically  predicted 
motions  of  the  Friday  Harbor  breakwater  are  shown  in  Figures  14,  15,  and 
16.  The  motion  response  is  almost  the  same  as  one  would  expect  from  an 
uncoupled  spring,  mass,  dashpot  linear  system.  The  only  unusual 
behavior  is  the  null  response  in  heave  at  B/L  of  about  0.75.  This  null 
occurs  at  a point  where  there  is  a phase  shift  in  the  "added-mass"  force, 
a phenomenon  which  has  been  observed  in  experiments  with  catamaran-type 
cross  sections  (Lee,  Jones, and  Bedel,  1971),  and  is  a result  of  resonant 
wave  conditions  within  the  open  well  of  the  catamaran. 

c.  Mooring  Line  Forces . In  recent  years  a great  deal  of  effort 
has  been  expended  in  understanding  and  predicting  mooring  line  perfor- 
mance, particularly  for  moored  ships  and  drilling  rigs  (e. g.,  American 
Society  of  Civil  Engineers,  1971).  While  many  of  these  analysis  techni- 
ques could  be  applied  to  the  moorings  of  floating  breakwaters,  this  has 
not  been  done  to  date.  There  are  also  very  few  model-scale  experiments 
in  which  mooring  forces  have  been  measured  and  only  a few  cases  where 
good  field  data  are  available. 

Two  techniques  for  calculating  the  spring  consults  for  mooring 
lines  have  been  used.  At  first  the  catenary  equations  were  applied  to 
find  the  change  in  force  per  unit  displacement.  While  this  approach 
leads  to  a fairly  simple  algorithm  for  the  calculation,  there  are  a few 
problems.  In  several  cases  spring  constants  were  needed  when  the  mooring 
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Theoretically  predicted  (twice  hydrodynamic  damping)  sway 
motion  response,  Friday  Harbor  breakwater. 


line  was  too  taut  to  allow  it  to  become  tangent  to  the  bottom  at  the 
anchor.  If  this  condition  occurs,  or  as  it  is  approached,  the  catenary 
equations  no  longer  apply.  For  many  full-scale  installations,  a combi- 
nation chain  and  synthetic  line  anchor  cable  is  used.  This  combination 
anchor  cable  presents  problems  in  attempting  to  use  the  catenary  equa- 
tions . 

Comparisons  between  the  mooring  line  forces  calculated  using  the 
catenary  equations  to  predict  spring  constants  showed  poor  agreement 
with  measured  results  (Adee,  1975).  While  the  general  trends  were  re- 
produced, an  increase  in  the  predicted  spring  constants  of  about  a factor 
of  4 would  have  been  required  to  bring  the  theoretical  prediction  into 
agreement  with  the  measured  results. 

To  overcome  the  problems  encountered  in  using  the  catenary  equa- 
tions, a system  based  on  discretization  of  the  mooring  line  and  static 
equilibrium  was  developed.  This  method  is  described  in  Appendix  B. 

(1)  Proposed  Oak  Harbor  Breakwater.  One  of  the  few  model 
tests  in  which  mooring  line  forces  were  measured  was  performed  by  David- 
son (1971)  for  the  floating  breakwater  proposed  for  Oak  Harbor,  Washing- 
ton. The  model  configuration  with  properties  scaled  to  the  prototype 
is  included  in  Appendix  F.  The  shape  of  this  breakwater  is  basically  an 
inverted  bathtub  with  foam  flotation. 

Applying  the  theory  to  predict  the  mooring  line  force  in  the  seaward 
anchor  line  at  a water  depth  of  29.5  feet,  one  obtains  the  results  shown  in 
Figure  17.  The  mooring- force  coefficient  is  defined  as  the  amplitude 
of  the  force  oscillation  divided  by  incident  wave  amplitude  times  the 
weight  per  unit  length  of  the  breakwater.  In  this  figure,  the  large 
range  of  the  experimental  results  is  directly  related  to  incident  wave 
amplitude.  The  smaller  incident  wave  amplitudes  generally  produce  lower 
measured  mooring  line  forces  per  unit  amplitude  except  at  the  beam  to 
wavelength  ratio  of  0.49.  Since  the  linear  theory  is  mathematically  cor- 
rect only  in  the  limit  as  wave  amplitude  tends  to  zero,  one  would  expect 
the  best  correlation  between  theoretically  predicted  and  measured  results 
for  small  amplitude  incident  waves.  The  results  shown  in  Figure  17  are 
consistent  with  this  expectation.  However,  the  very  large  difference  in 
mooring  line  forces  as  incident  wave  amplitude  increases  indicates  a 
highly  nonlinear  response. 

A potential  explanation  for  the  nonlinear  response  observed  in 
these  experiments  results  from  the  condition  of  the  mooring  lines  at  the 
29.5-foot  water  depth  used  for  the  model  tests.  Under  these  conditions, 
the  mooring  lines  no  longer  maintain  a catenary  shape.  When  the  initial 
tension  in  the  mooring  lines  is  increased  to  this  level,  they  respond 
with  very  large  changes  in  mooring  line  force  for  very  small  displacements 
of  the  breakwater.  Consequently,  small  deviations  in  the  planned  posi- 
tioning of  the  anchors  will  lead  to  large  changes  in  forces  in  the  moor- 
ing line.  This  condition  clearly  should  be  avoided  in  prototype  instal- 
lations where  very  large  mooring  line  forces  are  to  be  avoided. 
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Seaward  mooring  line  force  for  proposed  Oak  Harbor  breakwater 
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A second  possible  explanation  of  the  nonlinearity  results  when  the 
"drift  force"  on  the  breakwater  is  considered.  If  one  carries  the  hydro- 
dynamic analysis  to  second  order,  there  are  terms  at  zero  frequency 
which  yield  a force  on  the  breakwater  in  the  direction  in  which  the  in- 
cident waves  are  traveling.  This  force  has  the  same  effect  as  increas- 
ing the  initial  tension  in  the  mooring  line  and  is  proportional  to  wave 
amplitude  squared.  Increasing  the  initial  tension  tends  to  increase 
the  spring  constants  of  the  mooring  lines  leading  to  larger  oscillating 
forces  as  well. 

(2)  Alaska-Type  Breakwater.  Mooring-force  coefficients 
theoretically  predicted  and  measured  for  the  Tenakee,  Alaska, breakwater 
are  shown  in  Figure  18.  For  the  field  data  the  mooring-force  coefficient 
is  obtained  by  taking  the  square  root  of  the  mooring-force  spectral 
density  divided  by  the  incident  wave  spectral  density  and  then  dividing 
by  the  weight  per  unit  length  of  the  breakwater.  Again, as  with  the  Oak 
Harbor  model  experiments, there  is  good  agreement,  especially  in  predict- 
ing the  peak  in  the  curve  near  B/L  of  0.65. 

One  important  aspect  of  the  mooring  line  problem  which  should  not 
be  overlooked  is  a comparison  between  the  model-scale  results  and  the 
field  measurements.  For  the  Alaska-type  breakwater,  all  the  measured 
results  indicate  the  amplitude  of  oscillation  in  mooring  line  force  is 
in  the  order  of  hundreds  of  pounds,  not  thousands  of  pounds,  as  was  pre- 
dicted for  the  Oak  Harbor  breakwater  in  the  model-scale  tests. 

When  the  mooring  line  tension  data  recorded  at  Tenakee  are  plotted 
as  a function  of  time  as  in  Figure  19,  one  observes  that  there  clearly 
are  oscillations  associated  with  the  incident  waves.  However,  there  are 
also  low-frequency  oscillations  which  are  of  greater  magnitude.  A com- 
plete explanation  of  the  origin  of  these  low-frequency  forces  has  not 
been  developed.  However,  one  possible  explanation  is  that  these  forces 
are  a result  of  breakwater  oscillation  at  the  sway  resonant  frequency. 

Since  the  spring  constant  for  sway  motion  is  very  small,  one  would  ex- 
pect a long  natural  period.  Theoretically  predicted  sway  motion  response 
for  the  breakwater  is  plotted  in  Figure  20.  Predicted  natural  periods 
are  64,  37,  and  29  seconds  for  tidal  conditions  of  mean  lower  low  water 
(MLLW) , +10  and  +20  feet,  respectively.  By  applying  a high-pass  filter 
to  the  field  data,  one  obtains  the  spectrum  of  force  oscillation  shown  in 
Figure  4.  Here,  a peak  is  at  a period  of  about  53  seconds  (tide  height 
= +7  feet).  The  predicted  sway  natural  frequency  is  at  45  seconds  when 
the  tide  height  is  +7  feet,  which  indicates  that  this  explanation  is 
plausible. 

(3)  Friday  Harbor  Breakwater.  The  predicted  performance  of  a 
seaward  mooring  line  on  the  Friday  Harbor  breakwater  is  shown  in  Figure 
21  for  a tide  height  of  +5.33  feet.  The  Friday  Harbor  mooring  lines  are 
different  than  those  at  the  other  breakwaters.  They  are  composed  of  a 
section  of  chain  attached  to  the  breakwater,  followed  by  a length  of 
nylon  rope  and,  finally,  another  section  of  chains  at  the  bottom.  This 
particular  tidal  condition  was  chosen  because  it  is  the  condition  during 
record  FH  7-8  used  later  for  comparison. 
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Seaward  mooring  line  mooring-force  coefficient,  Tenakee,  Alaska 
(record  TK7-23) . 
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Friday  Harbor  breakwater. 


Low- frequency  predicted  sway  motion  and  resonance  are  shown  in 
Figure  22  for  MLLW,  +5.33  feet,  +10  feet,  and  +15  feet  tide  heights. 


+TI0E  HEIGHT  = +15.0  FEET 
ATI0E  HEIGHT  = +10.0  FEET 
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□TIDE  HEIOHT  = 0.0  FEET  (HUH) 


Figure  22.  Theoretically  predicted  long-period  sway  response 
Friday  Harbor  breakwater. 
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1 . Layout . 

The  site  of  the  floating  breakwater  instrumented  in  this  study  is 
located  at  Friday  Harbor,  Washington,  on  San  Juan  Island,  just  east  of 
Victoria,  British  Columbia  (Fig.  23).  The  breakwater  is  25  feet  wide, 
904  feet  long,  anchored  in  approximately  40  feet  of  water,  and  was  in- 
stalled in  October  1972.  The  structure  is  made  of  Polyolefin  flotation 
tanks  linked  together  by  a matrix  of  large  wooden  timbers.  It  is  laid 
out  in  an  expanded  L-shape,  the  inside  angle  being  115°,  with  the  short- 
er leg  (227  ft.)  directed  toward  shore  and  the  longer  leg  (627  ft.) 
toward  magnetic  north.  The  site  itself  is  protected  on  three  sides 
by  San  Juan  and  Brown  Islands  off  the  harbor  entrance.  This  leaves  an 
0.25-mile-wide  channel  into  the  harbor  with  a northeasterly  fetch  of 
about  1.7  nautical  miles.  Southeasterly  winds  can  also  generate  waves 
of  importance  parallel  to  the  shorter  leg  where  the  fetch  is  about  1 
nautical  mile. 


2 . Instrumentation. 

The  shorter  leg  was  instrumented  in  this  study  for  two  reasons: 

(a)  the  most  frequent  winds  are  out  of  the  southeast,  and  (b)  barges  were 
to  be  tied  to  the  longer  leg  during  the  winter  months  for  added  protec- 
tion. However,  the  wave  gages  are  positioned  to  give  the  proper  inci- 
dent and  transmitted  wave  data  for  ail  relative  wind  directions  (Figs. 

24  and  25)  . 

Four  types  of  time-dependent  data  which  are  basic  to  describing 
the  response  of  the  breakwater  were  collected:  (a)  wind  velocity  and 
direction;  (b)  wave  heights  at  key  locations;  (c)  anchor  cable  forces; 
and  (d)  directional  acceleration  and  angular  motions  of  the  breakwater. 
The  locations  of  the  measuring  sensors  are  shown  on  Figure  25.  Signals 
from  the  sensors  were  carried  by  underwater  cable  to  the  recording  system 
which  was  located  in  a small  building  mounted  near  the  center  of  the 
short  leg. 


3.  Wind  Data. 

Windspeed  and  di  ecticn  were  measured  by  Weather  Measure  Corpora- 
tion’s W121  sensor.  Some  additional  circuitry  was  required  to  record 
the  windspeed,  and  the  sensor  was  recalibrated  to  this  circuit.  The 
sensor  was  mounted  on  the  breakwater  at  the  intersection  of  the  two  legs 
at  20  feet  above  the  water  surface. 


4.  Waves. 

Wave  characteristics  were  measured  at  four  locations  with  the  second 
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transmitted  gage  being  used  as  a backup.  Two  spar  buoys  instrumented 
to  measure  wave  elevation  were  located  outboard  of  the  breakwater  and 
positioned  so  that  one  measured  the  incident  wave  field,  and  the  other 
measured  the  incident  plus  reflected  wave  field.  Two  stationary  gages 
were  attached  to  pilings  behind  the  breakwater  to  measure  transmitted 
wave  height.  All  four  gages  were  of  the  resistance  type.  The  spar 
buoys  were  used  outside  the  breakwater  to  help  reduce  navigation  hazards 
and  because  of  the  costs  and  logistics  of  placing  stationary  piling  at 
these  locations. 

The  buoys  were  made  of  two  sections  of  PVC  pipe,  the  lower  section 
being  6 inches  in  diameter  and  15  feet  long,  and  the  upper  section  of 
3-inch  diameter  and  12  feet  in  length  with  the  upper  8 feet  wound  with 
a resistance  wire.  Four  feet  were  exposed  above  the  water  surface,  and 
a 2 . 5- foot-diameter  disc  was  attached  to  the  bottom  to  damp  vertical 
motions.  The  natural  periods  in  heave  and  roll,  respectively,  are  18 
and  14  seconds,  well  above  the  anticipated  maximum  wave  period  of  about 
4 seconds.  See  Appendix  J for  a complete  description  of  the  wave  staff 
and  buoy  designs. 


5.  Cable  Forces . 

Anchor  cable  forces  were  measured  using  a bonded  strain  gage-type 
load  cell  that  was  placed  in  the  anchor  chains  beneath  the  water  surface. 
These  cells  and  the  associated  electronics  were  designed  and  built  for 
this  project.  They  have  an  overall  system  accuracy  of  0.75  percent  of 
the  designed  or  rated  total  load  cell  capacity  over  a temperature  range 
of  10°  Celsius  (design  load  12,500  pounds) , These  load  cells  employ  a four- 
arm  Wheatstone  bridge  circuit  which  has  two  strain  gages  in  each  leg  of 
the  bridge  and  are  self-temperature  compensating.  The  units  are  0-ring 
sealed  and  wired  directly  to  the  bridge  amplifier  circuitry  mounted  in 
the  recording  package. 


6.  Motion  Package. 

Breakwater  accelerations  were  measured  using  three  Kistler  servo- 
accelerometers  (Model  303T) . One  accelerometer,  oriented  horizontally, 
was  mounted  at  the  center  of  the  breakwater  to  measure  the  sway  accelera- 
tion. The  other  two  were  oriented  vertically  and  mounted  at  opposite 
outboard  edges  of  the  breakwater  to  measure  the  vertical  accelerations. 
The  heave  acceleration  was  obtained  by  taking  the  average  of  the  signals 
from  the  two  outboard  accelerometers;  the  roll  acceleration  was  obtained 
by  taking  the  difference  of  these  two  signals  and  dividing  by  the  distance 
between  them.  The  accelerometer  locations  are  indicated  in  Figure  25. 


7 .  Data  Acquisition  System. 

The  data  recording  and  electronic  package  was  built  around  the  Sea 


Data  Corporation's  Series  610  four-track  incremental  digital  cassette 
tape  recorder.  The  complete  package,  which  included  all  the  electronic 
circuitry  for  the  individual  transducers  plus  the  tape  recorder,  was 
housed  in  a watertight,  6-inch-diameter  PVC  cylinder  5 feet  in  length. 

The  system  was  designed  to  be  operated  manually  or  in  a completely 
automated  mode,  thus  requiring  only  periodic  tape  changes  (Fig.  26). 

In  its  automatic  mode,  the  system  was  activated  when  the  windspeed 
reached  or  exceeded  a preset  value  and  stayed  there  for  at  least  1 min- 
ute. At  this  point,  a single  17-minute  sample  of  all  the  inputs  was 
taken.  Each  68  minutes  following  this,  another  17-minute  sample  was 
recorded  if  the  wind  was  still  above  its  preset  value;  if  not,  the  system 
was  shutdown  until  the  windspeed  increased.  Each  17-minute  record  con- 
sisted of  2,048  samples,  taken  at  0.5-second  intervals,  of  all  13  chan- 
nels plus  a clock  channel.  Twenty-five  of  these  records  could  be  record- 
ed on  a single  cassette  tape. 


8 . Data  Processing  and  Analysis. 

The  initial  step  in  the  data  handling  was  to  transfer  the  data  from 
the  individual  cassettes  to  seven-track  magnetic  tape  by  means  of  the 
Sea  Data  reader.  These  tapes  were  then  converted  to  a computer  compati- 
ble format  on  the  University  of  Washington's  CDC  6400  computer.  The 
histograms  for  all  records  plus  the  basic  statistics,  i.e.,  the  minimum, 
maximum,  mean  values  and  standard  deviations  as  well  as  the  transmission 
coefficients  based  on  these  standard  deviations,  were  then  computed  and 
tabulated  (App.  G) . A digital  filter,  with  a cutoff  frequency  of  0.05 
hertz  (Gold  and  Radar,  1969)  was  applied  to  the  transmitted  wave  data 
prior  to  these  tabulations  to  remove  tidal  drift.  The  transmission 
coefficients  given  in  these  summary  sheets  are  a ratio  of  the  standard 
deviations  for  the  transmitted  and  the  incident  wave  gages. 

In  the  initial  conversion,  the  data  were  checked  for  reader  errors. 
These  points  were  smoothed  using  a linear  interpolation  between  the  pre- 
ceding and  the  following  good  data  points.  Following  this,  the  data  were 
checked  for  extreme  values.  Data  points  departing  from  the  mean  by  more 
than  five  standard  deviations  were  smoothed  in  the  same  manner  as  were 
the  reader  errors.  In  no  case  did  the  number  of  errors  warrant  elimina- 
tion of  a complete  record  (greater  than  six  bad  points).  Record  FH  11-1, 
however,  had  bad  data  for  channels  3,  4,  5,  7,  and  8.  This  record  was 
run  manually  while  calibrations  were  being  made,  and  the  affected  chan- 
nels were  not  connected  properly  at  this  time.  The  final  edited  data 
were  then  stored  on  magnetic  tape. 

The  autospectra  for  all  the  wave  data  for  all  records  were  computed 
with  a more  complete  analysis  of  the  force  and  acceleration  data  applied 
to  the  more  desirable  events. 

Digital  filtering  techniques  were  used  prior  to  spectral  analysis 
on  all  the  wave  and  force  data.  The  procedures  used  follow  those  given 
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Figure  26.  Instrumentation  and  recording  package  layout. 
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by  Gold  and  Radar  (1969).  The  first  step  in  the  development  of  this 
filter  function  is  to  assume  an  ideal  filter  response  function. 


Ffc(f) 

Fh(f) 


, 0 < |f|  < fc 

10  , f < If  I < f 

JO  , 0 < |f | < fc 

11  . f < Ifl  < f 


(12) 


where  f£  is  the  cutoff  frequency,  fn  the  Nyquist  frequency,  and  F ^ (f ) 
and  (f)  are  the  ideal  low-pass  and  high-pass  filter  response  functions. 

The  ideal  filter  response  function  is  then  Fourier-transformed  to 
the  time  domain,  giving  the  impulse  response  function,  which  is  trun- 
cated by  using  an  appropriate  window  function  and  transformed  back  to 
the  frequency  domain  givin0  a complex  frequency  response  function.  The 
number  of  points  used  in  representing  the  filter  function  is  allowed  to 
vary  and  the  resulting  convolution  with  the  original  time  series  is 
accomplished  by  using  the  overlap-add  method  of  convolving  smaller  ser- 
ies with  larger  ones.  This  allows  for  more  economical  filtering  proce- 
dures . 

This  gives  three  variables  to  choose  from  in  the  final  filter  func- 
tion design:  the  length  or  number  of  points  used  in  the  filter,  the 

type  of  window  used  to  truncate  the  impulse  response  function, and  the 
number  of  points  to  be  truncated. 


This  procedure  is  analogous  to  spectral  estimation  techniques  except 
for  the  truncation  of  the  impulse  response  function.  The  larger  the 
number  of  points  used  in  the  filter  function,  the  better  the  estimate. 

The  smoother  the  window  function,  the  broader  the  transition  band.  In 
addition,  the  ripple  or  Gibb's  phenomena  is  reduced.  Generally  speaking, 
the  more  points  that  are  truncated  (set  to  zero)  the  better  the  result- 
ing approximation.  In  practice,  the  actual  number  is  determined  experi- 
mentally by  comparing  results  for  different  truncation  values.  This 
results  in  setting  approximately  20  percent  of  the  impulse  response  func- 
tion to  zero.  The  hanning  window  function  was  used  with  128  points  in 
the  filter  response  function  and  38  points  being  set  to  zero  in  the  im- 
pulse response  function.  That  is: 


h(nAt) 

and 


w(nAt)  h(nAt) 


+ cos  it  2jj-) , 1 1.  n ^ 46 

, 45  < n < 83 

- cos  n ~ jj~~) , 83  < n < 128, 
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w(nAt) 


(13) 


where  h(nAt)  is  the  impulse  response  function  and  w(nAt)  is  the  hanning 
window  function.  The  final  filter  response  function  is  defined  as: 

F(fn)  = Z 6(nAt)e"-j2,TnfnAt, 
n 

where  j = and  At  is  the  constant  time  interval  between  samples. 

The  transition  band  or  the  frequency  increment  traversed  by  the 
cutoff  of  the  filter  function  can  be  approximated  by: 

. 10 

A = y2g  = 0.078 

and  the  maximum  stopband  attenuation  for  the  hanning  window  is  55  de- 
cibels. These  values  can  only  be  achieved  through  proper  filter  design. 
The  actual  values  for  the  filters  used  are  A =0.08  and  a maximum  atten- 
uation of  greater  than  55  decibels.  The  ripples  in  the  passband  for  each 
filter  used  were  below  0.01  percent.  These  values  could  be  improved  on 
by  increasing  the  number  of  points  used  for  the  filter  response  function 
estimate.  Also  the  stopband  attenuation  could  be  improved,  at  the  ex- 
pense of  a wider  transition  band  for  a given  size  filter  function  by 
using  the  Blackman  window  function.  However,  the  accuracy  of  the  filter 
response  functions  used  exceeds  that  of  the  measurements  and  is  suffi- 
cient for  this  application. 

After  initial  processing  and  prior  to  all  spectral  calculations  a 
tapered  cosine  data  window  was  applied  to  the  first  and  last  10  percent 
of  the  data  to  reduce  spillover  of  spectral  energy  to  adjacent  frequency 
points.  For  data  stretching  from  n = 1 to  n = N,  the  formulas  for  the 
data  window  are: 


w(nAt)  < 1 


cos  it  q— jyj-)  for  1 ± n <_  0.1N 

for  0 . IN  < n < 0.9N 


cos  n Q-~yN-)  for  0.9N  n N. 


The  data  were  then  transformed  directly  using  fast  Fourier  transformation 
procedures  and  smoothed  by  averaging  adjacent  raw  spectral  components. 
Initial  sampling  was  performed  at  0.5- second  intervals  with  2,048  samples 
per  record,  and  20  adjacent  points  averaged  together  in  the  autospectral 
calculations  to  get  the  final  smoothed  spectral  estimates.  This  gives  a 
frequency  resolution  of  0.0195  hertz  with  40  degrees  of  freedom  per  spec- 
tral estimate. 


All  of  the  wave  data  was  high-pass  filtered,  using  the  filtering 
techniques  previously  outlined  with  a cutoff  frequency  of  0.05  hertz. 

This  was  done  to  remove  the  tidal  influence  on  the  transmit5ed  wave 
staffs  and  to  eliminate  any  possilbe  buoy  motion  in  the  incident  wave  re- 
cords. Also  the  anchor  cable  force  data  were  separated  into  a low-and 
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high  frequency  signal  using  the  same  filtering  procedures.  For  the 
high-frequency  case  this  was  done  to  remove  the  influence  of  the  large 
low-frequency  spikes  in  the  spectra.  A high-pass  filter  with  a cutoff 
frequency  of  0.1  hertz  was  used. 

For  a closer  look  at  the  low-frequency  information  in  the  anchor 
cable  force  data,  a new  time  series  was  generated  from  the  original  re- 
cord by  sampling  every  eighth  data  point.  To  reduce  aliasing  of  the 
higher  frequency  energy  in  the  original  signal,  each  record  was  low-pass 
filtered  prior  to  this  sampling  using  the  filtering  techniques  previous- 
ly outlined  with  a cutoff  frequency  of  0.2  hertz.  The  sampling  of  every 
eighth  point  of  the  original  time  series  gives  a sampling  interval  of 
4 seconds,  a Nyquist  frequency  of  0.125  hertz  and  a record  length  of 
256  points  or  1,024  seconds.  Five  raw  spectral  points  were  averaged 
together  to  give  the  final  smoothed  spectra’  estimates.  This  results 
in  a frequency  resolution  of  0.0049  hertz  with  10  degrees  of  freedom  per 
spectral  estimate. 

A total  of  95  records  was  recorded  at  the  site  from  1330  hours  on 
30  December  to  3 May  1975.  There  were  no  known  equipment  failures  or 
breakdowns  except  for  one  of  the  load  cells  going  off  scale  at  low  tide 
on  the  first  tape  (FH  7,  NW  load  cell  channel  3).  A complete  summary  of 
these  events  is  given  in  Appendix  G.  Also,  Figure  25  gives  the  relative 
locations  of  the  individual  transducers. 

The  wind  direction  in  all  cases  is  referred  to  the  long  leg,  which 
has  a north-south  compass  bearing  (magnetic  declination  in  this  area  is 
23°  east).  There  are  two  wind-direction  windows  of  interest.  For  the 
long  leg,  the  directions  are  approximately  50°  to  95°;  for  the  short  leg, 
130°  to  160°  (Figs.  24  and  25). 

Two  storm  events  were  chosen  for  presentation  and  further  analysis. 
These  events  cover  records  FH  7-6  through  FH  7-12  and  FH  11-8  through 
FH  11-14  (Apps.  G and  H) . They  were  chosen  because  of  their  directions 
relative  to  the  short  and  long  legs,  respectively,  and  because  of  their 
duration  and  magnitude.  Both  events  lasted  for  over  7 hours  with  maxi- 
mum windspeeds  in  excess  of  35  miles  per  hour,  with  all  the  mean  wind 
direction  within  or  close  to  the  desired  wind-direction  windows.  Appen- 
dix H gives  the  pertinent  wave  spectra  and  transmission  curves  for  the 
above  two  events . 

The  average  overall  response  or  transmission  curves  for  the  events 
within  each  wind-direction  window  and  for  all  the  recorded  data,  are 
given  in  Figure  27.  These  plots  were  obtained  by  averaging  the  square  root 
of  the  ratio  of  the  transmitted  to  the  incident  wave  spectras  for  the 
records  indicated  for  each  curve.  Therefort , they  have  the  same  frequen- 
cy resolution  of  0.0195  hertz. 

A puzzling  feature  in  all  the  transmission  response  curves  calculat- 
ed from  field  data  is  the  rise  at  lower  frequency  to  a value  near  one  and 
then  dropping  off  again.  This  can  partially  be  attributed  to  a lack  of 
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FREQUENCY  (Hz) 

Figure  27.  Average  transmission  curves  for  Friday  Harbor  breakwater. 


energy  in  the  incident  wave  spectra  at  lower  frequencies.  This  possi- 
bility can  be  backed  up  in  part  by  data  from  two  similar  projects  under- 
taken in  the  past  2 years  on  a styrofoam-filled  concrete-type  break- 
water (Christensen  and  Richey,  1974).  The  first  is  located  in  relative- 
ly sheltered  water  while  the  second  breakwater  is  located  near  the  open 
ocean  where  swell  becomes  an  influence  and  results  in  a much  broader 
spread  of  spectral  energy  over  the  frequency  band  in  question.  The  first 
breakwater  showed  similar  response  curves  to  Friday  Harbor  while  the 
second  tended  to  approach  one  near  the  lower  frequencies  and  stay  there 
(see  App.  H.,  Figs.  H-10  and  H-ll). 

All  of  the  anchor  cable  data  showed  a very  dominant  amount  of  ener- 
gy at  lower  frequencies.  Appendix  I shows  the  results  of  the  low-fre- 
quency analysis  for  three  of  the  anchor  cable  force  signals  for  record 
FH  7-8.  The  autospectra  for  the  force  gages  show  several  large  peaks  in 
this  lower  frequency  band  (App.  1,  Figs.  1-1,  1-2,  and  1-3).  The  exact 
location  of  these  peaks  varies  for  different  records,  but  in  all  records 
analyzed,  the  dominant  amount  of  energy  in  the  force  spectra  was  con- 
tained in  this  lower  frequency  band  of  approximately  0.015  to  0.05  hertz. 
In  most  cases,  however,  a relatively  dominant  peak  appeared  in  the  56- 
to  63-second-period  range.  The  anchor  forces  measured  were  all  quite  low; 
the  largest  range  was  only  628  pounds.  The  cables  are  spaced  at  50- foot 
intervals . 

The  phase  and  coherency  spectra  for  three  of  the  force  gages  for  re- 
cord FH  7-8  are  given  in  Appendix  I (Figs.  1-4  through  1-7).  They  show 
a strong  linear  relationship  between  the  gages  on  the  same  side  of  the 
breakwater  and  for  the  opposing  gages.  The  forces  in  the  two  anchor 
lines  on  the  same  side  were  in  phase;  the  two  opposing  were  180°  out 
of  phase.  This  indicates  that  the  sway  or  roll  motions  are  dominant  in 
this  frequency  range.  The  accelerations  at  these  lower  frequencies  were 
too  small  to  be  recorded  and  could  not  be  used  to  help  confirm  which 
motion  was  involved.  However,  in  the  overall  frequency  range  (0  to  1.0 
hertz)  the  variances  for  sway  and  heave  were  two  orders  of  magnitude 
greater  than  roll,  which  indicates  that  sway  would  have  to  be  the  domi- 
nant motion  involved  here. 

The  analysis  of  the  complete  frequency  range  for  the  force  data  for 
FH  7-8  is  shown  in  Appendix  J (Figs.  J-2  through  J-8).  The  data  were 
high-pass  filtered  (fc  = 0.1)  for  these  spectra.  A comparison  of  the 
variances  computed  for  the  high-  and  low-frequency  sections  of  the  force 
spectra  showed  that  over  90  percent  of  the  energy  in  all  cases  analyzed 
was  contained  in  the  lower  frequencies.  A summary  of  the  force  data, 
without  any  filtering,  for  all  the  data  collected  in  this  experiment  is 
given  in  a table. 

The  autospectra  for  the  force  gages  (FH  7-8)  are  relatively  spread 
out  (see  App.  J,  Figs.  J-2,  J-3,  and  J-4) , with  the  outside  force  gages 
showing  a greater  response  to  the  lower  frequency  incident  wave  energy 
than  the  inside  gages.  However,  the  outside  and  the  opposing  gages  show 
relatively  high  coherency,  with  the  outside  gages  being  in  phase  and  the 
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outside  leading  the  inside  gage  by  approximately  130°  over  the  fre- 
quency range  of  0.25  to  0.37  hertz.  This  indicates  that  the  forces  are 
relatively  uninfluenced  by  waves  above  approximately  0.37  hertz.  This 
frequency  range  is  also  where  the  transmission  curves  rise  to  near 
unity.  This  agrees  with  the  low-frequency  analysis  and  suggests  that 
the  response  is  similar  over  the  complete  frequency  range  below  0.37 
hertz. 

The  acceleration  force,  autospectral  and  cross-spectral  analysis 
results,  are  also  given  in  Appendix  J for  the  higher  frequency  range  for 
record  FH  7-8.  No  dominant  features  were  observed  in  the  motion  spectra. 
Their  peak  values  and  spread  of  energy  with  frequency  appear  to  follow 
the  general  character  of  the  incident  wave  spectra  in  all  records  ana- 
lyzed. This  implies  that  any  natural  frequencies  in  each  of  the  motions 
is  outside  the  range  of  significant  incident  wave  energy.  The  cross- 
spectral  analysis  shows  a high  coherency  and  zero  phase  shift  between 
the  heave  and  roll  accelerations.  In  both  the  sway  and  roll,  and  the 
sway  and  heave  accelerations,  the  sway  acceleration  leads  by  approxima- 
tely 180°  over  the  range  of  significant  incident  wave  energy  and  then 
tapers  to  near-zero  phase  shift  at  higher  frequencies.  Also,  the  cohe- 
rency is  high  enough  over  the  incident  wave  energy  band  to  imply  near 
linearity  between  all  three  motions. 

These  conclusions  are  b'sed  on  positive  sway  being  outward  from  the 
short  leg  (south),  heave  positive  up,  and  the  positive  roll  to  be  clock- 
wise around  a positive  axis  pointing  westerly  toward  shore. 
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IV.  COMPARISON  OF  THEORY  WITH  FIELD  DATA  FOR  FRIDAY  HARBOR  BREAKWATER 


Although  the  Friday  Harbor  breakwater  has  a very  complex  geometry 
and  does  not  respond  as  a rigid  body  to  the  incident  wave  excitation,  it 
is  important  to  draw  some  comparisons  between  the  theoretical  prediction 
of  performance  and  the  field  measurements.  In  seeking  a "typical  event" 
from  the  enormous  quantity  of  data  gathered,  the  goal  was  to  find  a case 
where  the  wind  was  reasonably  close  to  being  on  the  beam  of  the  short 
leg  of  the  breakwater. 

The  one  striking  item  which  emerges  from  the  data  is  the  similarity 
of  all  the  transmission  coefficients  examined.  These  curves  seem  iden- 
tical no  matter  what  the  wind  direction.  This  was  not  expected  because 
there  were  barges  tied  to  the  breakwater  along  the  entire  long  leg,  while 
there  were  none  along  the  shorter  leg.  A further  investigation  of  the 
reasons  for  the  similarity  is  certainly  warranted. 

The  record  selected  for  comparison  with  the  theory  was  FH  7-8. 

Figure  G-3  in  Appendix  G shows  the  incident  and  transmitted  wave  spectra 
and  transmission  coefficient.  This  record  is  also  listed  in  the  statis- 
tical summaries  of  Appendix  F.  The  spectral  analysis  using  a high-pass 
filter  was  performed  as  described  in  Section  III. 

A comparison  of  the  theoretically  predicted  and  measured  trans- 
mission coefficient  is  shown  in  Figure  28.  So  long  as  the  calculated 
hydrodynamic  damping  is  doubled  in  the  theoretical  analysis,  the  results 
are  quite  good.  As  described  in  Section  II,  the  peak  in  the  transmission 
curve  at  a frequency  of  0.95  hertz  probably  results  from  the  "irregular 
frequency"  phenomenon  which  occurs  in  this  mathematical  formulation. 

Comparisons  of  sway,  heave,  and  roll  acceleration  predictions  with 
measurements  are  shown  in  Figures  29,  30,  and  31,  respectively.  Here, 
the  acceleration  response  has  been  nondimensionalized  by  multiplying  by 
the  beam  or  beam  squared,  as  appropriate,  and  dividing  by  the  accelera- 
tion of  gravity  times  the  incident  wave  amplitude. 

In  the  case  of  sway  acceleration,  the  theory  overpredicts  the  values 
throughout  the  entire  frequency  range.  The  peak  at  0.5  hertz  appears  in 
the  correct  location,  but  the  measured  values  would  need  to  be  doubled 
to  bring  the  curves  into  better  agreement. 

For  heave  acceleration  the  curves  appear  to  be  in  closer  agreement, 
at  least  above  the  frequency  of  0.4  hertz.  Below  0.4  hertz  there  seems 
to  be  little  correlation. 

Roll  acceleration  seems  to  show  the  worst  agreement  of  all.  Here 
again,  the  predicted  accelerations  are  considerably  higher  than  the  mea- 
sured values. 

There  are  several  possible  explanations  for  the  discrepancy  between 
predicted  and  measured  accelerations.  In  the  field,  even  if  the  wind 
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28.  Transmission  coefficient  for  Friday  Harbor  breakwater  (record  FH7-8) 


Figure  30.  Heave  accelerat 
(record  FH7-8) . 


were  blowing  directly  on  the  beam  of  the  breakwater,  one  would  not 
find  the  condition  of  long-crested  waves  impinging  directly  or.  the  beam 
of  the  breakwater.  As  a result,  the  breakwater  is  not  excited  uniformly 
along  its  entire  length.  Therefore,  the  breakwater  itself  provides 
restraint  against  motions  which  are  excited  in  a local  area.  The  con- 
struction of  this  particular  breakwater  is  also  quite  flexible,  which 
allows  for  considerable  internal  damping  of  the  wave-excited  motions. 

The  barges  tied  to  the  long  leg  also  serve  to  restrain  the  motion  and 
provide  additional  damping. 

There  is  a strong  need,  in  this  case,  to  provide  laboratory  data 
on  the  breakwater  motions,  which  could  be  further  correlated  with  the 
theory  and  the  measured  motions. 


If  one  looks  at  the  measured  accelerations  by  themselves,  a consi- 
derable resemblance  in  all  three  degrees  of  freedom  appears.  Further, 
if  these  accelerations  are  viewed  along  with  the  incident  wave  spectrum, 
considerable  similarity  appears  again,  suggesting  that  further  investi- 
gation of  the  measurement  scheme  would  also  be  welcome. 

The  final  comparison  to  be  made  is  between  the  theoretically  pre- 
dicted and  measured  mooring-force  coefficient.  The  theoretical  predic- 
tion and  measured  data  for  the  seaward  mooring  line  is  shown  in  Figure 
32.  The  correlation  appears  to  be  quite  good  in  this  case. 


In  looking  at  the  time  series  of  force  on  the  mooring  lines  and 
the  windspeed,  one  can  observe  a definite  correlation  between  the  wind 
gusts  and  increases  in  the  mooring  force.  This  is  probably  a result  of 
the  large  barges  tied  to  the  structure  which  act  almost  as  sails.  If 
this  is  the  case,  the  increase  in  tension  caused  by  the  mean  wind  on  the 
barges  needs  to  be  accounted  for.  No  attempt  has  been  made  to  do  this. 

The  most  common  method  of  presenting  the  spectral  data  obtained  in 
the  field  uses  a frequency  scale  rather  than  the  nondimensional  beam/ 
wavelength  scale  used  in  Section  II.  In  this  section  the  comparisons 
are  made  using  a frequency  scale.  For  the  Friday  Harbor  breakwater 
(beam  = 25  feet)  the  conversion  is: 


B 
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27iBf2 

g 


4.87f2 


assuming  deepwater  waves. 
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(record  FH7-8)  . 


V.  CONCLUSIONS 


Results  for  the  predicted  transmission  coefficients  were  in  good 
agreement  with  laboratory  and  field  data,  and  they  showed  how  the  in- 
fluence of  fixed-body  transmission,  and  of  sway,  heave,  and  roll  motions 
on  the  transmission  coefficient  changed  with  increasing  values  of  the 
beam  to  wavelength  ratio. 

The  curves  predicting  the  mooring  line  forces  as  a function  of  the 
beam  to  wavelength  ratio  (or  of  incident  wave  frequency)  followed  those 
for  the  measured  responses.  Care  must  be  exercised  in  the  analysis  of 
mooring  line  forces  because  there  is  strong  evidence  of  nonlinear  be- 
havior. 

An  extreme  storm  event  did  not  occur  during  the  sampling  season  at 
Friday  Harbor,  nor  during  two  winter  sampling  periods  on  the  Alaskan 
breakwaters;  however,  the  anchor  forces  measured  were  about  an  order  of 
magnitude  less  than  anticipated. 

The  barges  tied  to  the  long  leg  of  the  breakwater  did  not  noticeably 
affect  the  transmission  coefficients  above  a frequency  of  about  0.3 
hertz,  since  the  curves  for  all  incident  directions  were  approximately 
coincident  above  that  mean  frequency.  Below  the  frequency  of  0.3  hertz, 
it  appears  that  the  barges  may  have  reduced  the  transmitted  energy 
somewhat . 

The  extension  of  the  theoretical  model  to  include  second-order  terms 
showed  the  presence  of  additional  exciting-force  terms  at  zero  frequency 
and  at  the  difference  frequency  of  the  incident  waves.  Additional  work 
on  the  basic  theoretical  model  is  needed  to  incorporate  these  terms  into 
the  calculations  for  mooring  forces.  The  most  appropriate  means  of  veri- 
fying the  role  of  the  second-order  terms  may  be  in  a model  basin,  where 
breakwaters  of  simple  cross  section  and  incident  wave  spectra  having 
only  two  or  three  components  could  be  employed  under  controlled  condi- 
tions. 
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APPENDIX  A 


HYDROSTATIC  RESTORING  FORCES  AND  SPRING  CONSTANTS 


Hydrostatic  restoring  forces  and  spring  constants  are  computed  for 
the  two-dimensional  analysis  under  the  following  assumptions: 

(a)  The  body  rotates  about  the  origin  of  the  coordinate  system 
and  all  forces  and  moments  are  computed  about  that  point. 

(b)  The  body  has  vertical  sides  in  the  region  of  its  waterplane. 

(c)  All  motions  are  small. 

1 . Sway  Motion. 

In  the  horizontal  plane  the  body  is  in  neutral  equilibrium.  There- 
fore, there  are  no  hydrostatic  restoring  forces  and 

KHn  = KH12  = KH13  = 0.  (A- 1) 


2.  Heave  Motion. 


Vertical  displacement  of  the  body  results  in  a change  in  the  buoy- 
ant volume  of  the  body  and  consequently  a change  in  the  buoyant  force  on 
the  body.  Since  this  force  must  be  perpendicular  to  the  waterline, 
there  is  no  change  in  the  horizontal  force  as  a result  of  vertical  dis- 
placement and 

KH21  = 0.  (A-2) 

If  one  considers  a small  vertical  displacement,  6y,  there  is  a resulting 
change  in  volume: 


6V  = - <SyAw  (for  <$y  ♦ upwards) . 

Here,  A^  is  the  waterplane  area.  The  vertical  force  then  is: 

F = KH226y  = - PgAw6y, 
or 

m22  * pgAw  * pglxb  " xaJ‘  (A-3) 

In  this  equation  xa  and  X5  denote  the  sides  of  the  body  as  shown  in  the 
Figure  in  this  appendix.  Since  the  vertical  force  may  be  regarded  as  act- 
ing at  the  centroid  of  the  waterplane  area,  xc,  the  moment  may  be  expressed. 
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(A-4) 


M = K236y  = - pgAwxc6y. 

Substituting  for  x£  and  Aw  yields: 

^23  = - pSVc  = ' PgAw  7 lxa  + *b]  = 7 p^xa2  ' xb]  * 


3.  Roll  Motion. 

The  analysis  of  roll  motion- induced  forces  and  moments  is  compli- 
cated by  the  fact  that  the  body  is  assumed  to  rotate  about  the  origin 
of  the  coordinate  system  and  not  the  centroid  of  the  waterplane. 

The  problem  is  illustrated  in  the  figure.  Here,  line  2,  the  water 
line  after  rotation  through  and  angle  60  must  pass  through  the  inter- 
section of  the  y'  coordinate  axis  and  the  initial  waterline.  Equations 
for  lines  1 and  2 may  then  be  obtained. 

Line  1:  y = c. 

Line  2:  y = mx  + b. 

The  slope  of  line  2 is: 

m * ^7  = - tan  60. 

Ax 

Line  2 must  also  pass  through  the  point  P so  that: 

x = + c tan  60 
P 

and 


These  equations  yield  the  relationship: 
b = c(l  + tan260) . 

To  find  the  force  acting  on  the  body  as  a result  of  the  rotation, 
the  net  lost  or  gained  volume  is  needed. 

(mx  ♦ b - c)  dx 

[(-  x tan  60  c(l  + tan260)  - c]  dx 
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= - i [x  2 - x 2]  tan  60  + c[x,  - x 1 tan260. 

2 1 b a J b aJ 

By  applying  the  "small  angle"  approximation  and  neglecting  terms  of  the 
order  of  602.  Then, 

6V  5 i [x  2 - x 21  50 
2 1 a b J 

and  the  force  is: 

F ^ 1 pgtxa2  ‘ xb2l66* 

The  x and  y components  of  the  force  are: 

F = F cos  60  s pg[x  2 - x 2]  60  cos  60 

X i,  3.  D 

and 

F = F sin  60  = i pg[x  2 - x 2]60  sin  60. 
y 2 61  a b J 

Again  applying  the  small  angle  approximation  one  finds: 

Fx  SIpg[xa2  - Xb2]66 


and 

F ; 0. 

V 

The  hydrostatic  spring  constants  coupling  roll  to  sway  and  heave  are 
then: 

KH31  = 0 (A- 5) 


and 


*^32  = 2 pg^xa  ' xb  ^ 


(A- 6) 


To  obtain  the  moment  induced  by  roll  motion  compute: 

Moment  of  Gained  Volume  * - (i  x 2 tan  60) x ) ° i x ^ 60, 

Z c l J a J a 


Moment  of  Lost  Volume  5 j x^"*  60 


and 


Moment  of  Original  Volume  = Wy^  60. 


In  this  formula, 

W * weight  per  unit  length 
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and 


y^  = distance  the  center  of  buoyancy  is  below  the 
center  of  gravity. 

The  total  moment  then  is: 

M = 3*  (xb3  ' Xa3)69  + Wyb  66  * 
and  the  spring  constant  becomes: 


' ^ <‘b3  - ‘a3)  * “V 


(A-7) 


Expressed  in  traditional  naval  architecture  terminology,  this  reduces 


KH33  = WGM, 


where 


GM  = metacentric  height. 


Collected  Results. 


(A-8) 


l^ll  * ffll2  = mi3  = m21  = ^31  = ° 
^22  = Pg[xb  * Xa] 

^23  = ^32  =7pgtxa2  - \ ] 

*^33  = T^^b  ' xa  ^ + Wyb  * 


(A-9) 
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APPENDIX  B 


y'4 

MOORING  ANALYSIS 

1.  Purpose  of  the  Program. 

Computer  program  BRKMOOR  computes  the  forces  and  moments  imparted 
by  a pair  of  mooring  cables  on  a floating  breakwater  section.  BRKMOOR 
also  computes  the  changes  in  the  mooring  cable  tensions  and  the  spring- 
constant  values  for  the  moorings  as  the  breakwater  moves  in  sway,  heave, 
or  roll. 


2.  Program  Description. 

Program  BRKMOOR  is  written  primarily  in  FORTRAN  IV  although  FORTRAN 
II  print  statements  are  used. 

The  program  consists  of  the  main  program  BRKMOOR  and  the  subrou- 
tines LINE2,  CHAIN,  NYLON,  EQULIB,  SPRING,  and  LTERPS . 

BRKMOOR  calculates  the  forces  in  a mooring  cable  by  using  a discre- 
tized approximation  to  the  cable.  The  cable  is  divided  into  the  number 
of  segments  specified  in  the  input  data.  Each  segment  may  be  of  a diff- 
erent material  or  size.  Each  segment  is  in  turn  divided  into  a specified 
number  of  sections.  The  cable  is  considered  to  be  made  of  these  sections 
with  the  weight  of  each  section  concentrated  at  the  node  at  the  bottom 
of  the  section.  Connecting  each  node  i a straight  but  elastic  section. 

The  main  part  of  the  program  specifies  15  different  angles  at  the 
attachment,  ranging  from  nearly  vertical  to  nearly  straight  to  the 
farthest  reasonable  anchor  position.  A first  guess  at  a top  tension  is 
made. 

LINE2  then  sums  down  the  cable  computing  forces  and  coordinates  of 
each  node  starting  with  the  initial  angle  and  initial  tension.  The  po- 
sition of  the  end  of  the  cable  is  compared  with  the  specified  water  depth 
at  the  anchor.  The  initial  tension  is  adjusted  and  the  summation  repeated 
until  the  cable  ends  at  the  proper  depth.  Control  then  returns  to  the 
main  program. 

LINE2  calls  the  subroutines  NYLON  or  CHAIN  to  compute  the  strain  of 
the  cable  section  of  the  appropriate  material.  If  other  materials  are 
used  new  subroutines  should  be  written  for  strain  computation,  along 
with  the  appropriate  calling  expression  in  LINE2. 

At  each  angle  the  cable  forces  at  the  attachment  and  the  anchor 
position  are  stored  in  arrays.  EQULIB  then  computes  the  breakwater 
equilibrium  position  for  the  specified  conditions. 

SPRING  is  called  by  EQULIB.  SPRING  computes  the  change  in  mooring 
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cable  tensions  with  breakwater  displacement  in  sway,  heave,  and  roll 
and  the  spring  constants  of  the  moorings  on  the  breakwater. 

LTERPS  is  a linear  interpolation  subroutine  which  computes  the 
slop,  , and  the  interpolated  value  of  Y for  a given  X and  an  array 

of  X vsv  Y values.  LTERPS  is  called  by  EQULIB  and  SPRING. 

« 

3.  Type  of  Computer  and  Peripherals. 

BRKMOOR  was  written  for  use  on  the  CDC  6400  computer.  It  uses  about 
40,000g  words  of  memory.  No  peripherals  other  than  the  card  reader  and 
line  printer  are  required. 


4.  Input  Data. 

The  input  to  BRKMOOR  is  as  follows: 


Card  #1 
Card  #2 


Card  #3 
Card  #4 


Title  card.  Format  (8A10) . 

80  alphanumeric  characters  max. 

Breakwater  geometry  card,  Format  (5F10.0). 

YCG  = Vertical  location  of  breakwater  CG  relative  to 
water  surface. 

XCAB(l)  = x coordinate  of  cable  #1  attachment  to  break- 
water (the  CG  is  at  X = 0 and  cable  #1  is  de- 
fined as  the  cable  with  its  anchor  in  the  +x 
direction) . 

YCAB(l)  = y coordinate  of  cable  #1  attachment  to  breakwater. 
XCAB(2)  = x coordinate  of  cable  #2  attachment  to 
breakwater. 

YCAB(2)  = y coordinate  of  cable  #1  attachment  to  breakwater . 
Number  of  desired  conditions  Format  (12). 

(Also  number  of  condition  cards  to  follow) 

Condition  cards,  Format  (4F10.0). 

(One  card  for  each  condition) 

FEXT  = Force  applied  to  the  breakwater  not  due  to  moor- 
ings in  x direction  (could  be  due  to  wave  action, 
tide,  wind,  etc.  force  in  pounds). 

SEP  = Anchor  separation  in  horizontal  direction  (feet). 

TENS1  = Nominal  tension  in  cable  #1  (lb.). 

TENS2  = Nominal  tension  in  cable  #2  (lb.). 

It  should  be  noted  that  only  the  following  condition 
combinations  are  possible: 

SEP 

SEP+FEXT 

TENS1 

TENS1+FEXT 
TENS  2 

TENS2+FEXT 

TENS1+TENS2 
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Card  #5  - Tide  Card, Format  (I1,9X,5F10.0) . 

NT IDE  = Number  of  tide  values  to  follow  (max  = 5). 

TIDE  = Tide  position  in  feet  relative  to  that  at  which 
the  anchor  depths  are  given. 

Card  #6  - Cable  #1  Parameters,  Format  (I2,8X,2F10.0) . 

NSEG  = Number  of  different  segments  (types  of  cable  ma- 
terials) from  which  the  cable  is  constructed. 

DEPTH  = Depth  of  water  at  the  anchor  (feet) . 

BSLOPE  = Slope  of  bottom  in  region  of  anchor  (feet/feet). 

Card  # 7 - Cable  segment  properties  Format  (I5,5X,2F10.0,A10,F10. 0). 

One  card  for  each  of  the  number  of  segments  listed  in  card 

6 parameter  NSEG. 

NSECT  = Number  of  sections  into  which  it  is  desired  to 
divide  the  cable  segment. 

ALSEG  = The  length  of  this  cable  segment. 

WPF  = Weight  per  foot  in  water  of  the  cable  material  in 
this  segment. 

MATL  = Material  name  (as  the  program  now  stands  this 

must  be  CHAIN  or  NYLON  (Name  must  begin  in  column 
31). 

DIAM  = Diameter  of  the  nylon  rope  or  of  the  chain  link 
in  inches. 

Card  # 8 and  #9  - Same  as  cards  #6  and  #7  only  as  applies  to  cable. 

#2. 

Table  B-l  illustrates  the  input  cards  for  a test  case.  All  the 
read  statements  for  the  program  are  in  the  main  program  along  with  com- 
ments and  explanations  of  input  requirements. 


5.  Mathematical  Procedures  and  Program  Limitations . 

The  basic  cable  computations  which  take  place  in  LINE2  require  some 
explanation.  As  was  stated  previously,  the  weight  of  each  cable  section 
is  co-sidered  to  be  concentrated  at  the  bottom  of  the  section.  fn  order 
to  find  the  shape  of  the  cable,  summations  of  forces  are  computed  for 
static  equilibrium  at  each  node.  At  each  node  we  know  the  tension  in 
the  cable  section  above  the  node  as  well  as  the  angle  of  that  section 
with  the  horizontal.  Figure  B-l  illustrates  the  cable  about  the  ith 
node . 

If  the  angle  is  taken  to  be  the  angle  from  the  horizontal,  then 
the  angle  <^+1  can  be  computed  as  follows: 
x Ti  sin^i  + WA 

♦i*l  = tan  hrf  cos  -H 1 * (B-l) 

where  T^  = tension  in  section  i, 

W^  ■ weight  of  section  i concentrated  at  node  i. 
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MEASURED  CHAIN  TEST  3/11/76 
0.  -!•  °* 


I 

TEST 

0. 

Cfc 

0. 


1 

01 

00030 

01 

00030 


Mm 


CASE  - 

1. 

58.02 

53.21 

36. 
42  • 

54. 


36. 

30. 

7.167 

29.33 

.722 

CHAIN 

.25 

7.167 

29.33 

.722 

CHAIN 

.25 

Table  B-l.  Example  input  for  program  BRKMOOR. 


I 
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Cable  Section  i 


Node  i 


i + h — Node  i + 1 


Weight  " T. 


Figure  B-  1.  Cable  sections  about  node  i and  free 
body  diagram  of  node  i. 

This  new  angle  is  then  used  to  compute  the  tension  in  the  next 
section : 

1,  cos 

T.  = — i . (B-2) 

i+l  cos  *i+1 

LINL2  computes  the  angle  and  tension  of  each  section  starting  from 
the  top.  At  each  section  the  angle  is  compared  with  the  slope  of  the 
bottom.  When  the  angle  <p  is  parallel  or  more  positive  than  the  bottom 
then  <£  is  set  to  the  slope  of  the  bottom. 

The  x and  y coordinates  of  each  node  are  computed. 


i+l 

•H 

X 

II 

+ LliXT. 

i + l 

costi+1 

i + l 

= Y. 

sin<i>.  , . 

l 

F.XT.  , 
l + l 

i+l  * 

where  X^  = x coordinate  of  node  i 
Y^  = y coordinate  of  node  i 

* length  of  section  when  under  tension. 


At  the  last  node  the  /-coordinate  is  compared  with  the  depth  of 
the  anchor.  If  there  is  a difference  the  initial  tension  value  is  ad- 
justed. Guesses  at  the  first  and  second  tensions  arc  made.  From  then 
oil  a secant  (discrete  form  of  Newton  Raphson)  iteration  method  is  used 
to  compute  the  subsequent  initial  tension  values.  An  error  of 
0.0001*depth  is  allowed.  In  most  cases  4 or  5 iterations  yield  the 
desired  accuracy.  Some  important  values  are  printed  for  each  iteration 
to  aid  in  troubleshooting. 

Within  EQUL1B  and  SPRING  interpolation  is  required  to  find  the  va- 
lues of  tension  forces  and  x coordinates  which  are  between  the  points 
computed  by  BRKMOOR  and  LINL2.  The  linear  interpolation  routing  LTERPS 
was  chosen  over  higher-order  interpolation  schemes  because  of  the  asymp- 
totic nature  of  the  tension  versus  X values.  If  values  are  requested  be- 
yond the  ends  of  the  computer  arrays,  they  can  be  extrapolated,  but  a 
warning  message  will  be  printed  by  EQUILIB. 

An  iterative  procedure  is  required  within  EQULIB  if  the  anchor 
separation  condition  is  selected.  Again  the  secant  iteration  method 
is  used.  EQULIB  prints  out  values  at  each  intcration  which  can  aid  in 
troubleshooting  but  which  can  normally  be  ignored. 

Subroutine  CHAIN  computes  the  strain  in  a chain  using  the  basic 
elastic  properties  of  a steel  bar  with  a total  area  equal  to  the  area  of 
both  parts  of  the  links,  and  a factor  of  6 to  allow  for  the  deformation 
characteristics  of  the  links.  This  factor  of  6 came  from  a finite  ele- 
men t computation. 

Subroutine  NYLON  computes  the  strain  in  a nylon  rope  using  a power- 
function  fit  of  the  form: 

B 

e = AX  , 

where  e = Strain, 

A = 0.02052, 

B = 0.2237, 


T = Tension  (pound) , 

D = Diameter  of  rope  (inches) • 

This  function  was  determined  using  a least-squares  power- function 
fit  of  experimental  data  provided  by  Sampson  Cordage  Works  for  their 
2-in-l  nylon  braided  rope. 

An  experimental  verification  test  was  conducted  as  a check  of  the 
program.  A chain  was  suspended  from  a spring  scale.  Measurements  were 
made  of  the  length  of  the  chain,  its  weight  and  the  tension  in  two 
geometrical  configurations.  The  program  gave  computed  values  of 
the  tension  very  close  to  those  measured. 
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6. 


Flow  Chart. 


Figure  B-2  illustrates  the  flow  chart  of  BRKMOOR  and  its  subrou- 
tines . 


7 .  Program  Comments  and  Glossary  of  Terms. 

The  program  listing  contains  many  comments  which  aid  in  following 
the  logic  of  the  program.  The  important  variable  names  are  explained  as 
well  as  the  input  requirements. 


8.  Run  Time  and  Memory  Size. 


BRKMOOR  requires  about  40  seconds  on  the  CDC  6400  to  compile  and 
compute  results  for  one  value  of  the  tide  parameter.  Each  additional 
tide  value  requires  about  30  seconds  additional  time.  These  values  are 
for  cables  divided  into  50  sections  each.  Time  should  be  somewhat  pro- 
portional to  the  total  number  of  cable  sections.  The  number  of  test 
conditions  has  much  less  effect  on  time  than  does  the  tide.  As  stated 
previously,  a central  memory  of  about  40,000  octal  is  required. 

9 .  Run  and  Card  Deck  Setup  Procedures  and  Special  Operation 
Instructions . 


In  order  to  run  the  FORTRAN  source  program  deck  on  the  University 
of  Washington  CDC  6400,  the  following  deck  is  required: 


BMOOR.T40.  Job  card 

ACCOUNT  (Account  no.,  password) 

FORTRAN . 

LGO (LC=6000)  LC  = line  count  value;  depends  on  how  many 

tides  and  conditions  are  run 

7/8/9 

FORTRAN  DECK 

7/8/9 

DATA  DECK 
6/7/8/9 


10.  Sample  Output  Data. 

Example  output  from  program  BRKMOOR  is  shown  in  Table  B-2;  a listing 
of  program  BRKMOOR  is  shown  in  Table  B-3. 


LINE  2 


6 


LOOP  THROUGH 
EACH  CABLE 
SEGMENT 


PRINT  X,Y, 
TENSION  I 


\YES  , 
MATERIAL 
_ NYLON 


-""NOT  E =lVE  s / 
LAST  I'ERATION 


^ IF  \ 
MATERIAL 
_ CHAIN  _ 


-^ENO  OF\NO 
CABLE  AT  BOTTOM  J — 
\ DEPTH 


COMPUTE  LENGTH 
OF  SEGMENT 
UNDER  TENSION 


SET  NOTE=l 


COMPUTE  X * Y 
COORDINATES 
OF  EACH  NODE 

» 

COMPUTE  X & Y 
DIRECTION  FORCES 


COMPUTE 
NEXT  GUESS 
AT  INITIAL 
TENSION 


/IF\ 

LAST 

ITERATION, 

. (NOTE=l)  . 


PRINT  FORCES 
& COORDINATES 
ALONG  THE  CABLE 


COMPUTE  SLOPE 
OF  CABLE 


''ilOTTOM''' 
SLOPE  i CABLE 
\ SLOPE  / 


CABLE  PARALLEL 
TO  BOTTOM  & 
ASSUMED  TO 
REST  ON  BOTTOM 


LAST  \ 
.CABLE  SEGMENT 

7 > 


Figure  B-2.  Continued 
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SPRING 


GIVEN  X AND  Fx 
COMPUTE  DFXDX 


CALL 

LTERPS 


COMPUTE 

- / 

T,  DTDX 

.( 

Fy  DFYDX 

. 

PRINT  THE  ABOVE 

r 

CALCULATE 
FORCES  & MOMENTS 
ON  BRKWTR 


CALCULATE 
SPRING  CONSTANTS 
DUE  TO  MOORINGS 


CALL 

LTERPS 


) 


PRINT 
FORCES,  MOMENTS,! 
SPRING  CONSTANTS 


C 


RETURN 


3 


Figure  B-2.  Continued 
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IN  TNI  FOLLOWING  TABLES  - X 

REL.  TO  CG»  Y 

REL.  TO  WATFP 

SURFACE 

TEST  CASE  - 

- NFASURFD  CHAIN 

TEST  3/11/76 

800RINC  LINE 

NllNBER*  1 TIDE 

■ -.000 

Y 

X 

TOP  TENSION 

FORCEX 

FORCEY 

-7.167 

29.605 

5.202 

.930 

-5.189 

-7.167 

25.299 

5.802 

.956 

-5.722 

-7.167 

25.893 

6.9C7 

1.579 

-6.211 

-7.167 

26.387 

7.168 

2.326 

-6.770 

-7.167 

26.809 

8.101 

3.257 

-7.917 

-7.167 

27.179 

9.278 

9.919 

-8.158 

-7.167 

27.507 

10.776 

5.899 

-9.019 

-7.167 

27.819 

12.709 

7.809 

-10.021 

-7.167 

28.099 

15.259 

10.338 

-11.216 

-7.167 

28.366 

18.719 

13.777 

-12.669 

-7.167 

2B.619 

23.557 

18.601 

-19. 959 

-7.167 

28.839 

30.675 

25.695 

-16.755 

-7.167 

29.089 

91.699 

36.687 

-19.809 

-7.167 

29.290 

61.608 

56.999 

-29.690 

-7.167 

29.910 

121.396 

119.819 

-39.276 

-7 .239 

29.599 

5.251 

.939 

-5.233 

-7.739 

23.258 

5.898 

.963 

-5.768 

-7.239 

23.839 

6.999 

1.595 

-6.295 

-7.739 

26.396 

7.290 

2.353 

-6.897 

-7.239 

26.769 

8.190 

3.292 

-7.999 

-7.239 

27.192 

9.380 

9.968 

-6.298 

-7.239 

27.977 

10.900 

5.967 

-9.122 

-7.739 

27.790 

12.833 

7.888 

-10.122 

-7.239 

28.077 

15.917 

10.999 

-11.336 

-7.239 

28.397 

18.906 

13.919 

-12.795 

-7.239 

28.602 

23.805 

18.797 

-19.606 

-7.239 

28.899 

31.005 

25.971 

-16.935 

-7.239 

29.077 

92.115 

37.058 

-20.010 

-7.239 

29.277 

62.505 

57.266 

-25.050 

-7.239 

29.395 

125.036 

118.305 

-90.970 

-7.095 

29.658 

5.160 

.926 

-5.192 

-7.09J 

25.328 

5.757 

.998 

-5.678 

-7.093 

23.995 

6.327 

1.555 

-6.133 

-7.093 

26.927 

7.080 

2.301 

-6.696 

-7.095 

26.839 

8.016 

3.223 

-7.339 

-7.095 

27.203 

9.179 

9.372 

-8.071 

-7.093 

27.536 

10.657 

5.833 

-8.918 

-7.095 

27.837 

12.578 

7.732 

-9.922 

-7.095 

28.122 

15.085 

10.229 

-11.092 

-7.093 

28.386 

18.508 

13.626 

-12.525 

-7.093 

28.636 

23.319 

18.910 

-19.305 

-7.095 

28.879 

30.353 

23.923 

-16.379 

-7.093 

29.102 

91.256 

36.302 

-19.601 

-7.095 

29.303 

60.738 

35.697 

-29.391 

-7.093 

29.926 

117.873 

111.530 

-38.152 

Table  B-2.  Example  output  from  program  BRKMOOR. 
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TEST  CASE 


MEASURED  CHAIN  TEST  3/11/76 


MOOR  INC-  LINE 

NUMBER*  2 

TIOE*  - 

.000 

Y 

X 

TOP 

TENSION 

EflRCEY 

FORCEY 

-7.1/7 

-26. 

605 

5.202 

-.630 

-5.186 

-7.167 

-25. 

296 

5.802 

-.956 

-5.722 

-7.167 

-23. 

693 

6.607 

-1.576 

-6.211 

-7.167 

-26. 

387 

7.158 

-2.326 

-6.770 

-7.167 

-26. 

806 

8.101 

-3.237 

-7.617 

-7.167 

-27. 

176 

9.278 

-6.619 

-e.i;e 

-7.167 

-27. 

507 

10.776 

-5.899 

-9.019 

-7.167 

-27. 

816 

12.706 

-7.809 

-10.021 

-7.167 

-28. 

099 

15.256 

-10.338 

-11.216 

-7.167 

-28. 

366 

18.716 

-13.777 

-12.666 

-7.167 

-28. 

619 

23.357 

-18.601 

-16.656 

-7.167 

-28. 

859 

30.673 

-23.695 

-16.755 

-7.167 

-29. 

089 

61.696 

-36.687 

-19.809 

-7.167 

-29. 

290 

61.608 

-56.666 

-26.690 

-7.167 

-29. 

610 

121.366 

-116.816 

-39.276 

-7.239 

-26. 

566 

5.251 

-.636 

-5.233 

-7.239 

-23. 

238 

5.868 

-.963 

-5.768 

-7.239 

-23. 

839 

6.696 

-1.595 

-6.295 

-7.239 

-26. 

366 

7.260 

-2.353 

-6.867 

-7.239 

-26. 

769 

8.190 

-3.292 

-7.699 

-7.239 

-27. 

162 

9.380 

-6.668 

-8.268 

-7.239 

-27. 

677 

10.900 

-3.967 

-9.122 

-7.239 

-27. 

790 

12.833 

-7.888 

-10.122 

-7.239 

-28. 

077 

15.617 

-10.669 

-11.336 

-7.239 

-28. 

367 

18.906 

-13.919 

-12.795 

-7.239 

-28. 

602 

23.805 

-18.797 

-16.606 

-7.239 

-28. 

866 

31.003 

-23.971 

-16.935 

-7.239 

-29. 

077 

62.115 

-37.058 

-20.010 

-7.239 

-29. 

277 

62.505 

-37.266 

-25.050 

-7.239 

-29. 

393 

125.036 

-118.305 

-60.670 

-7.093 

-26. 

658 

5.160 

-.626 

-5.162 

-7.095 

-25. 

328 

3.757 

-.968 

-5.678 

-7.093 

-25. 

965 

6.327 

-1.353 

-6.133 

-7.093 

-26. 

627 

7.080 

-2.301 

-6.696 

-7.095 

-26. 

839 

8.016 

-3.223 

-7.339 

-7.093 

-27. 

205 

9.179 

-6.372 

-8.071 

-7.093 

-27. 

536 

10.657 

-5.833 

-8.918 

-7.095 

-27. 

837 

12.578 

-7.732 

-9.922 

-7.093 

-28. 

122 

15.085 

-10.226 

-11.092 

-7.095 

-28. 

386 

18.508 

-13.626 

-12.525 

-7.093 

-28. 

636 

23.316 

-18.610 

-16.305 

-7.093 

-28. 

876 

30.353 

-23.625 

-16.579 

-7.095 

-29. 

102 

61.256 

-36.302 

-19.601 

-7.093 

-29. 

303 

60.738 

-55.667 

-26.361 

-7.093 

-29. 

626 

117.875 

-111.530 

-38.152 

Table  B-2 
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BMOORUO  ////////////////////////////////////////////////////////////////////// 
BUNT  VERSICN  FEB  74  B 13i04  04/09/76 


PROGPAM  ERKNQORIINPUT, OUTPUT, PUNCH, TAPE  5-TNPUT, T APE fc -OUTPUT) 

C 

C PROGRAM  BRKHOOR  COMPUTES  THE  FORCES  AND  SPRING  CONSTANTS  THAT  A PAIR 
C OF  MOORING  CABLES  IMPART  ON  A FLOATING  BREAKWATFR  SECTION 

C 

C« 

C INPUT 

C FIRST  CAPO— TITLE  - 80  ALPHANUMERIC  CHARACTERS 
C BREAKWATER  6F0HETRY— 

C NUMBER  OF  TEST  CONDITIONS 
C TEST  CONDITIONS— ONE  CARO  FOR  EACH  SET 

C TIDE  CARD— NUMBER  OF  TIDE  CONDITIONS  AND  THE  CONDITIONS 
C FOR  FIRST  CABLE— NUMBER  OF  SEGMENTS  ANCHOR  DF°TM  AND  BOTTOM  SLOPE 
C FOR  EACH  OF  ABOVE  CABLE  SEGMFNTS  —CARD  WITH  SEGMENT  PROPERTIES 
C REPEAT  —NUMBER  OF  SEGMENTS  AND  THEIR  FROFERTTFS  FOR  SECOND  CABLE 

c 

3 COMMON/ONE /NSE6( 2 )*NSECT(2»5),WSECT(2» 5 ),  MATL ( 2, 5 > , D I AN ( 2, 5 ) , 

2 ALSECT(2,5> 

3 C0HH0N/TV0/WY!2,3),EX(2,3,20>»FX(2,?,?0).FY(2,3,20),TENS(?,3,20), 

2 FEXT(9),SEP(9),TENSH9),TENS2(9)»NANGLE.NC0ND,TITIE(8> 

3 C0MMrN/F0UR/YCG»XCABC2).YCABI2l»TTPE(f I.ITIOF 

3 DIMENSION  WPF(2,5),ALSEG(2,5),  PFPTH! 2 > , BSLOPE ( 2 > , ALNOM ( 2) 

3 PI-3. 1415926535 

C***READ  A TITLE  CARO  — 60  CHARACTERS  "A* 

5 READ  3. TITLE 

12  3 FORMAT! BA10) 

12  PRINT  16, TITLE 

20  16  F0RMAT!lMl,5X,BA10///> 

C***REAO  AND  ECHO  THE  BREAKWATER  GEOMETRY 
20  2 READ  5.YCG,XCAB(1),YCABI1),XCAB(2),VCAB(2> 

C YCG-Y  COORDINATE  OF  CG  RELATIVE  TO  WATER  SURFACE 
42  5 FORMAT! 5F 10.0) 

42  OUTPUT, YCG.XCAB III , YCAB! 1 ) , XC ABI 2) ,YCAR<2» 

C XC AS « I ) -X  COORD  OF  CABLE  I ATTACHMENT  RELATIVE  TO  CG 

C YCAB ! I ) *Y  COORD  OF  CABLE  I ATTACHMENT  RELATIVE  TO  WATER  SURFACE 

C NOTE— CABLE  NUMBER  1 IS  THE  CABLE  WITH  ITS  ANCHOR  IN  THE  + X DIRECTION 
C«*MNPUT  THE  NUMBER  OF  DESIRED  CONDITION  CAROS  MAX  NUMBER-9 
72  READ  10,NC0ND 

100  10  FORMAT! 12) 

C***READ  AND  ECHO  DESIRED  CONDITIONS 
100  00  17  IC0ND-1,NC0ND 

102  READ  15,FEXTIICONO),SEP!ICCND),TENS1IICOND),TENS2!ICOND) 

121  15  F0RHATI4F10.0) 

121  17  OUTPUT,  TEXT (ICO NO), SEP! ICCND ) , TENS1! ICOND ), TENS 2 1 ICONO ) 

C FEXT-EXTERNALLY  APPLIED  FORCE  ( HORI 20NTAL  DIRECTION)  LB. 

C SEP  .ANCHOR  SEPERATION  IN  THE  X DIRECTION  FT, 

C T ENS 1 -TENSION  IN  CABLE  1 (B. 

C TENS2-TENSI0N  IN  CABLE  2 LB. 

C INPUT  SFP,  OR  TENS1  OR  TENS2  CR  BOTH  TFN$1  AND  TENS2 
C«*PREAD  AND  ECHO  TIDE  CONDITIONS 
C NTIOE-NUMBER  OF  TIOE  CONDITIONS  MAX-5 

C TIOE-TIPE  POSITION  RELATIVE  TC  NOMINAL  OFPTH  MEASUREMENTS  FT. 

151  READ  20,  NTIOE,(TIOE(I),I-1,NTIDE) 

166  20  FORMAT! II, 9X,5F10.0) 


Table  B-3.  Listing  of  program  BRKMOOR. 


93 


PUNT  VERSICN  FEB  7*  8 


13*0+  04/09/76 


166 

2C1 

703 

717 


717 

773 

773 

??« 


730 

764 

764 

326 

376 

331 

332 
343 
357 
37? 

377 

403 

405 

406 
411 


411 


471 

423 

431 

443 

444 
44  5 
446 
456 
461 
463 
470 


OUTPUT. NTI0E.TI0E 
C LOOP  THP0U6H  THE  TWO  CABLES 
00  65  1-1,2 

C**4 INPUT  THE  CABLE  PROPERTIES  ANO  BOTTOM  DEPTH  AND  SLOPE 
PE  AO  2?.K$EGM)»DEPTHU).RSL0PE(n 
27  FORMAT, 12, 8X.2F 10.0) 

C NSEG-  NUMBER  OF  CABLE  SEGMENTS  OR  MATFPJALS 
C DEPTH-  DEPTH  OF  THE  HATER  AT  THE  ANCHOR  FT. 

C BSLOPE-  SLOPE  OF  THE  BOTTOM  (FT  RISE/FT1 
PRINT  25 

75  FORMAT! // /5R» * I NUMBER  SECTIONS  SEGMENT  LENGTH  WT  PER  F00T4 
7 4X4  MATERIAL  DIAMETER*/) 

NS-NJEGU) 

00  30  J-l.NS 

C***FOR  EACH  CABLE  SEGMENT  INPUT 

C NSECTO-  NUMBER  OF  SECTIONS  INTO  WHICH  CABLE  SEGMENT  I IS  OIVIOED 
C ALSEGin-LENGTH  OF  CABLE  segment  I FT. 

C WPF (II-  HEIGHT  PER  FOOT  IN  WATER  OF  CA8IF  SEGMENT  J LB/FT 

C MATL-MATER1AL  OF  CABLE  SEGMENT  FITHFP  Nrt ON  CR  CHAIN 

C MUST  BE  LEFT  JUSTIFIEO  IN  CATA  FIHO 
C OIAN-OIAMFTER  OF  ROPE  OR  CHAIN  LINK  INCHES 

RE  AO  40.NSECT! It J)  » At.  SEG( I.J).WFF(I.J).  M ATI ( I, J > . 01  AM ( I . J ) 

40  FORMAT! I5.5X, 2F10 .0. A10. f 10.0 ) 

30  PRINT  50, J,NS6CT(I,J»,ALSFG(I. J),WPF(I,J),  **ATL(  I»  J » ,0  IAM  1 1,  J» 

•0  F0RMAT(X»I5»8X»I5»8X,FlO.?»4X»F10«2,4X»A10»  5X.F6.3) 

C4P4FIN0  THE  NOMINAL  LENGTH  OF  THE  CABLE,  LENGTH  AND  HEIGHT  OF  SECTIONS 
ALN0N(Il-0. 

00  60  J- 1 ,NS 

ALN0M(II-AIN0M<I)4AISEG( I,J) 

ALSECT(I,JI-ALSEG(I, J)/NSECT(I, J) 

<0  WSECT(I,J»-WPFtl,J)*ALSECT(I,J) 
t 5 CONTINUF 

BSL0PE(2>— 8SL0PEI7) 

C*44L00P  THPOUEH  THE  TIDE  POSITIONS 
00  400  IT  IDE -1,NTI0E 
C***LOOP  THROUGH  THE  CABLES 
00  1*0  1-1,2 
PRINT  70 

70  F0RMAT(1H1,5X*N0TE— IN  THE  FOLLOWING  TABLES  X ANC  Y APE  MEASURED 
7PELATIVF  TO  THE  CABLE  ATTACHMENT*/ / I 
C D-V  DIRECTION  SEPFR AT  ION  BETWEEN  ANCHOR  ANn  ATTACHMENT 
D-OEPTH(I)+TIDE(ITIOEI+YCAB<I) 

C*4*C0MPUTE  INITIAL  ANGLES  TO  BE  USEO 
C NANGLE-NUMBER  OF  ANGLES  USF3  MAX-70 
NANGLS-15 

PHIMIN- AS  INI  0/ ALNOM( 1 1 ( 

0ELPHI-(PI/2.-PHIMIN»/FL0AT(NAN6LF4l> 

PHIONF— PI/2. 

C444C0HPUTE  A FIRST  GUESS  FOR  THE  INITIAL  TENSION  FOR  STEEPEST  ANGLE 
AL SUM-0. 

T2FR0-0.  -f 

OAF-0+1  ALNOMd  I-0I4BSLOPEI  I) 

NS-NSEGfll 
00  90  /-I ,NS 
NSS-NSECT ( I. J I 
00  40  K-I.NS? 


Table  B-3.  Continued 
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SUNT  VEPSItN  FEB  7*  B 13iO*  0W09/76 


ALSUM-AL$UN*ALSECTII,J) 

I F ( SI  SUM  .GT.  DAF>  GO  TO  95 
90  TZER0-T7fR0*VSECT< I,  J» 

95  CONTINUE 

C CONFUTE  the  NOMINAL  AND  PERTUSBED  DEPTHS 
ORFF-D 

DELD-DREF/100. 

DPLUS-ORFF*OELD 

DMINUS-OPFF-DELD 

C***LOOP  THROUGH  THE  INITIAL  ANGIES 
00  100  K-l.NANGLF 
PRINT  97,I,K,NANGLE,TIDE!ITIDE> 

*7  FORMAT! /X*CABLE  NUMBER  *11*  INITIAL  ANGLE  NO.  *12*  OF  *12 
? * TIDE  • *F 5 ,2/ ) 

RHIONE-PHIONE*OElPHI 

C***LOOP  THROUGH  THE  NOMINAL  DEPTH  AND  PEPTURBED  OEPTHS 
DO  100  J-1,3 
IPRINT-0 

C**P*****TO  SKIP  THE  PRINTING  OF  EACH  CATINAPV  - INSERT  A GO  TO  1111 
IF  C J ,E0.  1)  IPRINT-1 
1311  IFIJ  ,E0.  1)  D-DREF 
I F < J .EO.  2)  D-DPLUS 
IF ( J .EO.  3)  O-DMINUS 
WY!I,J)-YCABII>-D 

OUTPUT, J,K,D,YCABI II .WY( I > J ) * I ON  E 

CALL  LINF2CI,PHI0NE,TZER0,r,BSLCPFm,Y,Y,F0RCEX,F0RCEV,lPRIET» 

EX (I, J,K)«XCABII»-X*!-1)**1 

FXiI,J,K)-FORCEX*i-»-l)**I> 

FYII.J.KJ-FORCEY 

TENS!I,J,KI-TZERO 

C WY ( I , J ) ■ Y COORD  OF  THE  ANCHOR  TO  NO.  1 CABLE— HATER  SURFACE-ORIGIN 
C EX(I,J,K)-X  COCRD  OF  ANCHOR  RELATIVE  TC  CG  OF  BREAKWATER 
C TEMS-TENSION  AT  ATTACHMENT 
C FX-FORCE  AT  ATTACHMENT  IN  X DIRECTION 
C FY-FORCF  AT  ATTACHMENT  IN  Y DIRECTION 


712 

100 

C 1 

716 

150 

720 

72* 

102 

72* 

726 

733 

103 

733 

7*6 

105 

7** 

7A6 

7A7 

1013 

110 

1013 

120 

1020 

1023 

125 

1023 

1025 

160 

END  OF  CABLE  LOOP 
CONTINUE 
PRINT  102 


X RFL.  TO  CG,  Y RFL . TO  HA 


2TER  SURFACE*//! 

DO  160  1-1,2 
PRINT  103, TITLE 
FORMAT!  5X» BA10  > 

PRINT  105,1,TIDE(ITIDE) 

FORMAT! ///5X*N00RING  LINE  NUMBER-  *11,*  TT0E-*F6.3// 

2 10X,*Y*1AX,1HX,BX,*T0P  TENSI 0N*7X*F0PCEX*7X*FGRCE Y*/ ) 

DO  120  J-1,3 
DO  120  K-1,NANGLE 

PRINT  110,WYII, J),EX(I,J,K),TENS(T,J,K|,FXIT, J,K),FY(I,J,K) 


PRINT  125 


CALL  EOULIB 


Table  B-3.  Continued 
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• ... 


HINT  VEP5ICN  FES  74  S 


13<04  04/09/76 


C FND  OF  TIDE  LCO» 
1026  *00  CONTINUE 

1031  STOP 

1033  ENO 


Table  B-3.  Continued 


RUNT  VERSION  FES  74  B 13«C4  04/09/76 


SUBROUTINE  LINE2IK.PHI0NE.TZER0* DEPTH*  B SLOPE* X* Y*  FORC  ER» FORCE Y» IP ) 
19  C0MH0N/0NE/NSE6I2 ) *NSECT !2*S)*WSECTI2»9)*  MATLI2* 9 > * DI AH  12. 9 >* 

2 ALSECTI 2*9  > 

C**THE  INPUT  TO  SUBROUTINE  LINE 
C FMIONF. INITIAL  ANGLE  OF  CABLE 
C TZERO-INITIAL  GUESS  OF  TENSION  AT  TOP  OF  CABLF 
C**SUBROUTINE  LINE  CONFUTES 
C TZFRO-  TENSION  AT  CABLE  TOP 
C FORCER-FORCE  IN  R DIRECTION  AT  CABLE  TOP 

C FORCEV-FORCE  IN  Y DIRECTION  AT  CABLE  TOP 

C R-HORIZONTAL  SEPERATION  BETWEEN  TOP  ANO  BOTTOM  OF  CABLE 

C Y-VERTICAL  SEPERATION  BETWFEN  TOP  ANO  BOTTOM  OF  CABLE 

C***60  DOWN  THE  CABLE  SECTION  BY  SECTION  COMPUTE  TENSION* ANGLE* 

C ERTENDED  LENGTH. R AND  Y COORDINATES 
19  PI-3.14199 

16  NITER-0 

17  MNITER-29 

21  NOTE-O 

22  T-TZERO 

23  192  NITER-NITER+1 

29  IF  1 IP  .EO.  C)  GO  TO  193 

27  IF! NOTE  .EO.  0)  GO  TO  193 

31  PRINT  195 

39  199  FORMAT! //9R* I J R V TENSION  LSECT* 

2 6R**LERT  PHI-OEGREES  FCRCEY  FORCER*/) 

39  193  V-0. 

37  R-0. 

43  PHI-PHIONE 

44  NSS-NSEGIK) 

47  198  DO  200  I-1»NSS 

91  NS-N$ECT(K*I) 

96  00  200  J-l.NS 

97  PHt0-PHI*lB0./PI 

61  IF IMATL !K» I ) .EO.  9HNYL0N  )60  TO  169 

67  IFIMATLIK* I)  .EO.  9HCHAIN  )G0  TO  160 

79  160  CALL  CHAINIDIAMIK*I)#T»S TRAIN) 

109  GO  TO  170 

111  169  CALL  NYLON !DIAMIK» I)*T*S TRAIN) 

121  170  ALERT-ALSECTIK*I)*I1.*STRAIN) 

134  R-R*ALERT*COSI PHI ) 

149  Y-Y*ALERT*SINI  PHI ) 

192  TCOS-T*COSCPHI) 

196  TSIN-T*SIN(PHI) 

162  IF! IP  .EO.  0)  GO  TO  1B9 

170  IFINOTE  .EO.  0)  GO  TO  1B9 

172  PRINT  1 90*1* J*R*Y*T*ALSECT(K*I)> ALERT. PH TD.TSIN.TCOS 

231  180  FORMAT! R»2I9»8F10. 3) 

231  189  IF  1 1 .EO.  NSS  .ANO.  J .EO.  NS)  GO  TO  200 

247  SL0PE-!TSIN4NSECT(K»I))/TC0S 

297  IF ISLOPE  .CE.  BSLOPE)  SLOPE-BSLOPF 

262  PHI-ATANI SLOPE) 

266  T-TCOS/COSIPHI ) 

271  200  CONTINUE 

302  FCRCER-TZ£RO*COSI PHIONF ) 

311  FORCE Y-TZERC*S INI PHIONE ) 
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t 


u 

320  OUTPUT, NITER. Y.X.TZERO 

C»**THE  SECONO  GUESS  OF  INITIAL  TENSION  IS  CONFUTED 
350  IF (NITER  .GT.  1)  GO  TO  220 

353  TZOLD-TZERO 

354  TZERO-TZERO*ABS(DEPTH/Y> 

365  VOLD-Y 

367  T-TZERO 

370  GO  TO  152 

C*94THE  SUBSEQUENT  INITIAL  TENSIONS  ARE  CONFUTED  USING  SECANT  ITERATIC 
370  220  RELER-ABS(1.*V/DEPTHI 

402  IF (NOTE  .EO.  1)  GO  TO  300 

405  IF (NITER  .GE.  NNITER  .OR.  RELER  .IF.  .0001  I NGTE-1 

424  DEROLD-DEPTMVOLD 

426  DERR-0EPTM*Y 

430  T-TZOlD-DEROlD*(TZERO-TZOlC)/)DERR-DrRClD) 

437  IF(T  .IE.  0.)  T-TZERO/2. 

443  YOLD-V 

445  TZOLC-TZERO 

446  TZERO-T 

447  GO  TO  152 

447  300  RFTUPN 

450  END 


JNT  VERSION  FEB  74  B 13104  04/C9/76 


\ SUBROUTINE  CHAIN(D.T» STRAIN) 

6 PI-3.14139 

7 E-30.E6 

11  AREA-0 *D*PI/2. 

C C-ELONGATION  FACTOR  — C-6  FOR  OVAL  CHAIN 

14  C-6. 

15  $TRAIN-C*T/( AREA*E) 

21  RETURN 

22  END 


UNT  VERSION  FEB  74  B 13«04  04/09/76 


SUBROUTINE  NYLON(D» T. STRAIN) 

j 

6 

X-T/(0*D) 

10 

A-. 02052 

* 

12 

B- .2217 

13 

STRAIN-A*X**B 

20 

RETURN 

21 

END 
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SUBROUTINE  EQUIIB 

2 C0HMPN/TW0/wyi2.3>,EXI2,3»20>.FXI2.3,2C>.FYI2,3»20),TENSI2,3,20>» 

2 FEXTI9)»SEPI9)»TENSH9)»TENS2I9) , NANGLE.NCOND,  TI Tit ( 8 ) 

2 COMMON/ THREE /XI2)»FI2) 

2 DIMENSION  SE  PDIF  I 3 ) , FOI 3 ) 

C****EOUl IB  FINOS  THE  BREAKWATER  EQUILIBRIUM  POSITION 


C***LOO*  THROUGT  THE  TEST  CONDITIONS 


2 

00  100  IC*1#NCONO 

A 

TFISEPIIC)  .NE.  O.t  GO  TO  20 

7 

IFITENS1I1C)  .NE.O.IGO  TO  10 

13 

IFITENS2IICI  .NE.O.IGO  TO  12 

17 

PRINT  155 

23 

155 

FORMAT C //X*NO  INITIAL  CONDITIONS  S»ECIFIEr*> 

23 

GO  TO  ICO 

C***FOR  THE  CASES  WHERE  INITIAL  TENSION  IS  GIVEN  THE  FOLLOWING 

IS  USEO 

2* 

10 

T.TENS1IIC) 

27 

1*1 

31 

J«2 

32 

GO  TO  1A 

33 

12 

T«TENS2(ICI 

36 

r«2 

AO 

j*i 

A1 

1A 

OUTPUT. I, WANGLE, T 

57 

IFIT  .GE.  TENSII. 1.1)1  GO  TO  18 

70 

PRINT  16 

7A 

GO  TO  100 

75 

18 

IF(T  .GE.  TENSII. l.NANGLE)  I PRINT  17 

111 

16 

FOR NAT  1 //5X* GIVEN  TENSION  CLOSE  TO  0»  LESS  THAN  WEIGHT  OF 

VERTICAL 

2 MOORING  LINE*/5X*N0  FURTHER  EVALUATION  ATTFNPTEO  • //> 

111 

17 

F0RMATI//5X*GIVEN  TENSION  TOO  GREAT  FOR  FVALUATION  WITHOUT* 

2 *EX7RAPOlATION*/5X*USE  RESULTS  WITH  CAUTION*//) 

111 

CALL  LTERPS  1 I . l.NANGLE. TENS. EX. T. XI I ) .DUMMY I 

123 

OUTPUT, XII), T 

137 

CALL  LTERPS  I I, l.NANGL E, TENS.FX.T, F 1 I ) .Dummy) 

151 

OUT PUT, F II) 

162 

IF  1 1 .EO.  1 .AND.  TENS2IIC)  .NE.  0)  GO  TO  12 

177 

IFITENS1IIC)  .NE.  0 .AND.  TENS2IIC)  .NE.  0)  GO  TO  AO 

21A 

FI  J).-  FID-FEXTIIC) 

22A 

OUTPUT.FI J) 

23A 

CALL  LTERPS  IJ, l.NANGLE, FX, EX, F IJ) . X IJ ) .DUMMY) 

250 

OUTPUT, X 1 J ) 

C 

NOTE—  FID-X  DIRECTION  FORCE  ON  CABLE  I , XII). X CODEC  OF 

ANCHOR 

261 

GO  TO  AO 

C***FOR  THE  CASE  WHERE  ANCHOR  SEPERATION  IS  GIVEN 

C 

MAKE  A FIRST  AND  SECOND  GUESS  AT  FORCE 

262 

20 

TA>INAN6LE*ll/2 

270 

EPS" SEPIIC)*. 0001 

27A 

DO  30  11.1,2 

275 

XI1).EXI1,1, IA) 

305 

F 1 1 ) >FX 1 1, 1, IA  ) 

315 

OUTPUT.FfD.II.FEXTIIO.SEPI  IC),1A 

3AA 

F0IIII.FI1) 

351 

F 1 2 ) — F I 1 l-FEXTIIC) 

361 

OUTPUT, FI2) 

371 

CALL  LTERPS  12, l.NANGLE, FX, EX, F 1 2) , XI2 ), DUMMY ) 
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40$ 

413 

421 

466 

476 

477 

24 

503 

505 

505 

26 

507 

30 

c*** 

511 

513 

514 

540 

552 

557 

567 

570 

32 

605 

613 

621 

657 

667 

672 

677 

704 

711 

716 

34 

720 

36 

724 

37 

724 

38 

726 

750 

39 

775 

42 

775 

43 

775 

40 

777 

100 

1002 

1002 

ASEP*X<1)-X(2> 

SEROIFCII )-SEP<IC)-ASEP 

OUTPUT, II. IA.X<1),X<2),F!1),F<2),*EP!IC>»ASEP/SEPDIF III) 
IF<ABS<SEPDIF<  ID ) .GT.  EP!)  GO  TO  24 
GO  TO  40 

IF< SFPDIF < II ) .61.  0.)  GO  TO  26 
IA-1 

GO  TO  30 

IA-NANGLE 

CONTINUE 

JSE  SECANT  INTERPOLATION  FOR  THE  SUBSEOUFNT  FORCE  TRIALS 
HN-20 

00  34  K-l.NN 

FO<3)-FO<1)-SEPDIF<1>*<FO<2)-FO<1))/<SFPOIF<2>-SEPDIF<1>> 

IF  <F0( 31  .LE.  0.)  FO  < 3 ) -FO  < 2 ) /2. 

F<l)-F0<3) 

F<  21  •**F  1 1 )-F  EXT<  IC  ) 

DO  32  1-1.2 

CALL  LTERP5  <1. 1. NANGLF. FX, EX. F < I > , X < I ) .DUNHV ) 

ASEP-XIl )“X<  2) 

SEPDIF<3)-SEP<ICI-ASEP 

OUTPUT,*, X<1),X<2),F<1),F< 2), ASEP»SEPDIF<3) 

IF <ABS<SEPDIF< 3))  .LE.  EPS)  GO  TO  3P 
IF  < K .EO.  MN)  GO  TO  36 
FO  < 1 ) -F0<  2 ) 

FO  <2 ) -FO ( 3 ) 

SEPDIF<1)-SEPDIF<2> 

SfPDIF<2)-SEPDIF<3) 

CONTINUE 
PRINT  37 

FORMAT! /5X*MAX  NUMBER  OF  ITERATIONS  REACHED*/ ) 

DO  39  1-1.2 

IF  < ABS<  F < I ) ) .GT.  ABS<FX<I,1.NANGLE)>)  PRINT  42 
IF<ABS<F(I>>  .LT.  ABS<FX< 1.1.1)))  PRINT  43 

FORMAT! //5XPANCH0R  SEPERATION  TOO  6REAT  FOR  EVALUATION  WITHOUT  EXT 
2RAP0LATING— USE  RESULTS  WITH  CAUTION*//) 

FORMAT! //5X4ANCH0R  SEPERATION  TOO  LITTLE  FOR  EVALUATION  WITHOUT  EX 
2TRAP0LATI0N--USE  RESULTS  WITH  CAUTION*//) 

CALL  SPRING! IC ) 

CONTINUE 

RETURN 

ENO 
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6 

6 

6 

6 

6 


6 

7 

10 

27 

35 

54 

62 

107 

133 

140 

147 

174 

221 

223 

230 

230 

261 


261 

265 

265 

270 

337 


337 

345 

352 

360 

365 


SUBROUTINE  SPRIN6IIC) 

COMMON/TWO/W Y!  2»3>»EXI2»3»20)»FX!2»3*2!))»FV!2»3»20)»TENSI2»3*20)» 

2  FEXTI9)»SEPI9)»TENS1!9),TENS2I9> »NANGLE»  NCOND,  TITLE  < 8) 
CPMN0N/TH*EE/X»2>,FI2> 

COMMON/FOUR /YCG,XCAB!2)*YCAB!2),TTDF 15 ),ITIDF 

DIMENSION  DFXDX!2>»0FYDX!2)»0FX0YI2)»DFY0Y!2)»D!3)»FXX12,3)» 

2 FYX!2,3)*FV!2)*0TDX!2)»DTDY!2>»DTDRI2),T!2»3) 

REAL  KMH»KN12»KH13»KH21»KN22»KH23»KM31»KN32,KM33 
C**** SUBROUTINE  SPRING  COMPUTES  THE  BREAKWATER  SPRING  CONSTANTS 
C***COMPUTE  THE  SPRING  CONSTANTS  FOR  EACH  CABLE 
C HORIZONTAL  FORCE  AT  EQUILIBRIUM-FII)  FOR  CABLE  I 
C VERT  FORCE  AT  E OUILIBRIUM-FV! I)  FOR  CABLE  I 
DO  14  I-l»2 
DO  12  J-1,3 

CALL  LTERPS  1 I, J,NAN6LE» EX, FX, X ( II  * FXX <1# J ) ,0F > 

I F ( J .EO.  1)  DFXDXtlJ—  DF 

CALL  LTERPS  ! I * J.NANGLE, EX, TENS* XU ) »T ( I , J ) ,DT > 

IF < J .EO.  1>  DTDXII)— DT 

12  CALL  LTERPS  !I*J»NANGLE»EX»FV»XII)»FYXII»J)*DIJ)) 
DTDYII)-!TII,3)-Tll»2>  >/!WY!I,2)-WY!I»3>) 

OFYOXIII  — 0(1) 

FV« I » -FVX  €1*11 

DFYDY!II-!FYX!I,2)-FYXII*3n/!WYII.2>-NY!I,3) )*l-l.) 
0FXDYII)-IFXXII,2)-FXXtI,3n/!WYlt.2)-WYI  [,3IIM-1.I 

14  CONTINUE 
PRINT  16, TITLE 

It  FORMAT! 1H 1,8 AlO) 

PRINT  15,TIDEIITIDE»,FEXT« IC > , SEP! IC ) . TENS1 ! IC ) , TENS2« IC ) 

15  FORMAT! ///X*FOR  THE  CONDITIONS  —*/ 

1 5X*TIDE  • *F5 .2/ 

2 5X*EXTERNALL Y APPLIED  HORIZONTAL  FORCE,  FEXT-  *F1C.3*LB.*/ 

3 5X*H0«IZ0NTAt  ANCHOR  SEPf RATION,  SEP-  *F10.3*FEET*/ 

4 5X*N0MINAL  TENSION  IN  CABLE  1 -*F1©.3*  LB.*/ 

5 5X*N0MINAL  TENSION  IN  CABLE  2 ■*F10.3*  LB.*//) 

PRINT  IB 

18  FORMAT! /5X*C ABLE  NO.  FX*10X*FY*IIX,1HX,11X*Y*9X*DFXDX*7X, 

2 *DFX0V*7X*DFY0X*7X*0FY0Y*,5X,*TENSI0N*//) 

DO  20  I>1,2 

20  PRINT  25,  I,FII),FV!I),X!I),WY(I,1),DFXDX(I),0FY0X!I),DFXDYII), 

2 DFYDVI I ),T 1 1, 1 > 

25  F0RNATI4X , 11,4 i 2X, F10.3 ) ,81 X,E 11.3  D 

C«**NOW  CALCULATE  FORCES  AND  SPRING  CONSTANTS  FOR  THE  BREAKWATER 
C S-SWAY  MOTION  RX  DIRECTION  FE'T 
C H-HFAVE  NOTION  *Y  LIRECTION 
C R-ROLL  NOTION  COUNTERCLOCKWISE  RADIANS 
C FS-FORCES  CAUSING  SWAY  DUE  TO  THE  MOORING  LINES 
C FH-FORCES  CAUSING  HEAVE  DUE  TO  MOORING  LINES 
C ENR. MOMENTS  CAUSING  ROLL  DUE  TO  NOORINC  LINES 
C CHANGE  YCAB  TO  BE  OIST  TO  CC  IN  Y OIRTCTION 
YCAB11 )«YCAB I 1 l-YCG 
YCAB!21-YCAB!2)-YCG 
FS-FI1)*F12) 

FH-FVI 1) rFVI 2) 

EMR-FVIl)*XCAB«l)*FVI2)*XCABl2l-F!l ) *YCAB 1 1 )-F « 2 )* YC AB I 2 ) 
C«**CALCULATE  CHANGE  IN  TENSIONS  WITH  BREAKWATER  MOTIONS 
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413 

414 
432 
435 
435 
440 
461 

461 

471 

500 


530 

537 

546 


576 

626 

656 


756 

767 

767 

1001 

1001 

1013 

1013 

1025 

1025 

1031 

1031 

1032 


00  26  1-1.2 

26  DTDR<Il-CTOY<Il*XCAB(Il-OTOXtIl*YCA8(I) 

BRINT  27 

27  F0RHAT(//5X*CABLE  NO.  DTDX*8X*DTD V*8X*DTDR*/ / ) 

00  28  1-1.2 

28  88 INT  29,I,CTDX<I1,DTDY<I1»DTDRCI1 

29  F0RHAT(9X.I1.3(XE11.4li 

C S881NG  CONSTANTS  SWAY  DIRECTION 

KH11-<OFXOX(1>+OFXOX<2> I*<-1.  I 
KN12-(OFYDX<11*OFYOX<2I l*(-l.l 

KH13-(DFYDX( 11*XCAB( 11  *OFYDX(  2 1 *XC  AM  2 1-OF  XflX  ( 11*YC  A6  ( 1 1-DF  XDX  t 2 1 ♦ 
2 YCABI2! >*(-1.1 
C SPRING  CONSTANTS  HEAVE 

XH21- (DFXDY (1>*DFXDY(2>>*(“1.> 

KH22-(0FYDY(1>*DFYDY(2>>*(-1.> 

KH23-(OFYOY(1>*XCAB(1>*OFYDY(2>*XCAB(2> 

2  -0FX0rm*YCABU)-DFX0YI2)*YCAB(2))*l-l.> 

C SPRING  CONSTANTS  80U  DIRECTION 

XN31-(0FXDY(1)*XCA6(1>*DFXDV(2>*XCAR(2>-DFXDX(11*YCAB(1> 

2 -0FX0X<2>*YCAB(2>1*(-1.I 

KH32-(0FY0Y(11*XCAB(11*0FYDY(21*XCAP(2) 

2 -OF YDX( 1 1* YC  AB( 1 >-DFYDX(2 >*YCAB(21I*(“1.1 
KN33-tXCAR(l>**2*DFYDY(l>*XCAB(2>**2*DFY0Y(2> 

? ♦YCAB(11**2*DFXDX(1>*YCAB(21**2*DFXDX(2> 

3 “XC AB ll)*YCAS(l)*(DFYOX( 1) ♦OFXOY( 1 1 1 

4 -XCAR(21*YCAB(2>*(DFYDX(21*0FXDV(2>>  >*(-1.1 
PRINT  30.FS.FH.ENR 

30  F0RHAT(///5X*FQRCES  AND  HOPENTS  ON  BREAKWATFP  AT  EQULIBRIUH  DUE  * 

2 PTO  HOOPING  LINES*/10X*FS»  PF12.4/10XPFM-  PF12.4/10XPH8-  *F12.4) 

PRINT  32.M11.KH12.KP13 

32  FORPATI //5X* SPRING  CONSTANTS  SWAY  0IRECTI0N*/10X*KHI1  • *£12.5/ 

2 10X*K«12  • P€12.5/10X*KH13  >*£12.51 

PRINT  34.KM21.KN22.KH23 

34  F0RHAT(/5X*SPRING  CONSTANTS  HEAVE  0IRECTICN*/10X*KH21  - *£12.5/ 

2 10X*KN22  • *£12.5/10X*KH23  >*E12.51 

PRINT  36/XH31.KH32.KP33 

36  FORHATt /5X*SPR ING  CONSTANTS  ROLL  DT*ECTI0N*/10X*KH31  - *£12.5/ 

2 10XAKH32  > *E12. 5/10X«KN33  >*E12.5//1 
PRINT  38 
38  F0RNATI1H1I 
RETURN 
END 
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13 

SUBROUTINE  l TERRS  < l > J»N, X, Y> XX, YY,PYOx ) 
OINENSIPN  XC2,3,20),Y(2,3,20) 

13 

NNO-N-i 

15 

PO  10  K* 1 , NNQ 

16 

20 

IF  ( XX  .EO.  X(I,J,l>)  GO  TO  30 

27 

IFUBS<XX>  .IT.  ABS<X(I,J,l)>>  GO  TO  20 

55 

10 

CONTINUE 

60 

20 

OYOX-(V(I, J,LI-YII,J,K))/(X(I, J»l)-x(I, J,K) ) 

111 

YY-YII, J,KJ*(XX-X< I,J,K)»*PYOX 

131 

RETURN 

132 

30 

IF  It  .EO.  N)  GO  TO  20 

136 

N-L61 

136 

PVOX-»Y«I, J,N»-Y(I,J,K))/«XU, J,N)-X<I,J»K) ) 

167 

VY»YU»J#l> 

176 

RETURN 

176 

END 

5 
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APPENDIX  C 


LINEAR  HYDRODYNAMIC  COEFFICIENTS 


The  linear  theoretical  model  used  in  solving  the  floating  break- 
water problem  has  been  discussed  extensively  by  Frank  (1967) . He  de- 
veloped the  approach  to  solving  the  boundary  value  problem  which  has 
come  to  be  known  as  the  "Frank  close-fit  method".  The  reader  is  re- 
ferred to  the  original  reference  for  a complete  presentation  of  the 
method . 

In  this  approach,  the  classical  linear  boundary  value  problem 
requires  that  Laplace's  equation  be  satisfied  throughout  the  fluid  do- 
main: 


V2$(x,y,t)  = 0 for  y < 0.  (C-l) 

The  free-surface  boundary  condition  is  applied  on  the  undisturbed  free 
surface: 


4'tt(x,0,t)  + g$v  = 0 for  y = 0. 


(C— 2) 


The  body- surface  boundary  condition  requires  that  no  fluid  flow  through 
the  body  surface: 


VHx.y.t)  • n 


Vi(s)  • n(s) . 


(C-3) 


The  bottom  boundary  condition  for  infinite  depth  is  of  the  form: 
y^L  '•yCx.y.t)  = 0. 


(C-4) 


In  addition  there  is  a radiation  condition  specifying  that  the  waves 
travel  away  from  the  body. 

Because  the  problem  is  assumed  to  be  linear,  the  velocity  potential 
may  be  decomposed  and  several  boundary  value  problems  considered.  If 
this  is  done  the  total  potential  becomes: 

<t  * 4^  ♦ ^ ♦ *3  ♦ *4  ♦ $5.  (C-5) 

Here, 

4^  = potential  representing  pure  sway  motion  in  calm  water, 

4>2  * potential  representing  pure  heave  motion  in  calm  water, 

<t>3  = potential  representing  pure  roll  motion  in  calm  water. 


i>  = potential  representing  the  waves  diffracted  by  a fixed 
body, 

$5  = incident  wave  potential. 

Another  velocity  potential  may  be  defined: 

= potential  for  total  fixed-body  problem, 

so  that 


*6  = *4  + V 

Using  this  decomposition  of  the  velocity  potential,  the  boundary 
value  problems  may  be  expressed  as:  I 

V2<ti(x,y,t)  = 0 for  y < 0, 

<t>.  (x,0,t)  + g$.  = 0 for  y = 0, 

Xtt  v 

lim  4> . (x,y,t)  = 0,  (C-6) 

7 y 

and 


V4>.  • n 
i 


= V^(s)  • n(s)  for  i = 1,2,3 


or 


V4K  • n 


0 for  i = 4,6. 


These  boundary  value  problems  are  solved  directly  using  the  Frank  method 
which  distributes  singularities  over  the  hull  surface.  These  singulari- 
ties satisfy  the  radiation  condition,  Laplace's  equation,  the  free- 
surface  boundary  condition  and  the  bottom  boundary  condition.  To  satis- 
fy the  body  boundary  condition  requires  the  formulation  of  a set  of 
linear  equations  whose  solution  reveals  the  strength  of  each  singularity 
distributed  on  the  body. 


Once  the  velocity  potential  is  found  the  pressure  may  be  found  from 
Bernoulli's  equation: 


P(x,y,t)  = - p*t(x,y,t) . 

The  force  on  the  body  surface  is: 
P(s)  n(s)  ds, 

9 

and  the  moment  is: 


(C— 7) 


(C-8) 
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P(s)[rxn]  ds. 


(C-9) 


«« 


M 


The  added-mass  and  damping  coefficients  are  found  by  considering 
the  cases  i = 1,2,3.  The  forces  and  moments  computed  using  these  po- 
tentials may  be  separated  into  components  in  phase  with  acceleration 
and  velocity.  The  component  in  phase  with  acceleration  yields  the 
added-mass  coefficients  and  the  component  in  phase  with  velocity  yields 
the  damping  coefficients.  Exciting  forces  and  moments  are  computed 
when  the  case  i = 6 is  considered. 

Special  Symbols  for  Appendix  C. 

n(s)  = unit  interior  normal  vector  to  the  body  surface 

s = indicates  arc  ldhgth  along  body  contour 

CQ  = body  contour 

P(s)  = pressure  on  body  surface 

V(s)  = velocity  of  body  surface 

$ = total  velocity  potential 
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APPENDIX  D 


FLOATING  BREAKWATER  ANALYSIS 

1 . Purpose  of  the  Program. 

Computer  program  BRK2D  performs  a performance  analysis  for  two- 
dimensional  floating  breakwaters  of  arbitrary  cross  section.  This  ana- 
lysis includes  predictions  of  the  hydrodynamic  coefficients,  the  dyna- 
mics and  mooring  line  forces. 

2.  Program  Description. 

Program  BRK2D  is  written  using  both  FORTRAN  II  and  FORTRAN  IV 
statements . 

The  program  consists  of  the  main  program  BRK2D  and  the  subroutines 
COEFF,  COMP,  PHYSCL,  POTOUT,  DYNAMC , MORTEN,  CPV,  LNEQF . 

The  subroutines  COEFF  and  COMP  calculate  the  quantities  needed  to 
formulate  the  linear  equations  for  the  velocity  potential.  COMP  calls 
on  LNEQF  to  solve  these  linear  simultaneous  equations. 

Subroutine  PHYSCL  calculates  the  physical  quantities  including  ad- 
ded-mass  and  damping  coefficients  and  surface  elevations  per  unit  ampli- 
tude of  motion. 

CPV  is  a subroutine  which  evaluates  the  Cauchy  principal  value 
integral  in  the  Green  function. 

LNEQF  is  a packaged  subroutine  to  solve  simultaneous  linear  equa- 
tions using  the  Gaussian  reduction  method. 

3.  Type  of  Computer  and  Peripherals. 

BRK2D  was  written  for  use  on  the  CDC  6400  computer.  It  uses  about 
55000g  words  of  memory.  No  peripheraJs  other  than  the  card  reader,  line 
printer  and  card  punch  are  required. 

4.  Input  Data. 

The  first  cards  in  the  data  deck  are  label  cards  for  the  output. 
These  are  shown  in  the  example  input  in  Table  D-l  for  the  example  and 
are  not  included  here.  Following  these  cards,  the  input  for  BRK2D  is: 

Card  #1  - Title  card.  Format  (8A10) . 

80  alphanumeric  characters. 

Card  #2  - Logical  control  card,  Format  (5110,615). 

N = Number  of  straight  line  segments  used  to  fit  the  hull. 
NW  * Number  of  points  on  the  free  surface  where  wave 
height  is  to  be  computed. 


NWAVEL  = Number  of  wavelengths  at  which  computations  are 
to  be  performed. 

ISYM  = 1 for  symmetric  section. 

= Anything  else  for  non- symmetric  section. 

ISKIP  = 1 Do  not  solve  equations  of  motion, 

2 Do  not  solve  potential  problem  (read  in 
coefficients  from  data), 

= Anything  else  solve  potential  problem  and  equa- 
tions of  motion. 

LC  = Number  of  body  segments  which  represent  spaces  be- 
tween multiple  hull  configurations  (1  to  5). 

JC  = Designates  the  segment  number  for  segments  repre- 
senting spaces  between  multiple  hulls. 

Card  #3  - Parameter  card.  Format  (5F10. 3,3A10) . 

AREA  = Crossectional  area  of  immersed  body. 

B = Characteristic  beam  as  specified  by  BTITLE. 

D = Distance  below  free  surface  to  origin  of  users 

coordinate  system  (all  motions  are  referred  to  that 
point) . 

ROE  = Fluid  density. 

GEE  = Acceleration  of  gravity. 

BTITLE  = Specifies  B. 

Card  #4  - Beam/wavelength  specification.  Format  (10F8.5). 

BOL  = Beam/ wave length  ratios  for  computation  (up  to  10 
different  ratios  may  be  used). 

Card  #5  - Offset  cards.  Format  (2F10.3). 

There  must  be  N+l  cards  giving  the  offset  points.  In  the 
version  of  the  program  used  here,  N must  be  less  than  or 
equal  to  23  because  of  dimension  statements. 

R(  1 , 1)  = X-coordinate  of  offset  point. 

R(2, I)  = Y-coordinate  of  offset  point. 

Card  #6  - Hydrostatic  spring  constants.  Format  (9F8.3). 


This  is  read  in  subroutine  DYNAMC. 

RKHYD 

(l.D 

= 

KHn. 

RKHYD 

(1,2) 

= 

KHi2. 

RKHYD 

(1,3) 

= 

KHi3. 

RKHYD 

(2,1) 

= 

KH2i. 

RKHYD 

(2,2) 

= 

KH22- 

RKHYD 

(2,3) 

= 

kh23. 

RKHYD 

(3,1) 

= 

KH3i. 

RKHYD 

(3,2) 

= 

K»32. 

RKHYD 

(3,3) 

= 

KH33* 

Card  *1  - Physical  properties.  Format  (6F10.3,3F5.2,I5) . 

This  is  read  in  subroutine  DYNAMC. 

AREA  = Crossectional  area. 

B = Characteristic  beam. 

XG  = X-coordinate  of  the  center  of  gravity. 

YG  = Y-coordinate  of  the  center  of  gravity. 

RMASS  = Mass  per  unit  length  of  breakwater. 

RINERT  = Mass  moment  of  inertia  per  unit  length  of 
breakwater. 
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DAMP(l)  = Added  damping  in  sway.  In  the  equations  of 

motion  sway  damping  will  be  1+DAMP(1)  times  the 
computed  hydrodynamic  damping. 

DAMP (2)  = Added  damping  in  heave. 

DAMP (3)  = Added  damping  in  roll. 

tfPUNCH  = 0,  punch  data  cards  containing  computed  trans- 
mission coefficient,  motion  response  and  mooring- 
force  coefficient. 

= Anything  else,  do  not  punch  data  cards. 

Card  #8  - Mooring  spring  constants,  Format  (9F8.3). 

This  is  read  in  subroutine  DYNAMC. 


RKMOR(l.l) 
RKM0R( 1 , 2) 
RKM0R( 1 , 3) 
RKM0R(2, 1) 
RKM0R(2, 2) 
RKM0R(2,3) 
RKM0R(3, 1) 
RKM0R(3,2) 
RKM0R(3 , 3)  = 
Mooring- line 
This  card 
DELT (1,1) 


KMU. 

KMi2. 

KMl3- 


= KM 


21* 
km22. 

km23. 


= KM 
= KM 


31' 

32* 


is 


KM33. 

Card  #9  - Mooring-line  response  parameters,  Format  (6F10.2). 

is  read  in  subroutine  MORTEN. 

= AF/Ao^  for  shoreward  mooring  line.  This 
the  change  in  mooring  line  force  per  unit 
displacement  in  sway. 

DELT(1,2)  = AF/Aa2  for  shoreware  mooring  line. 

DELT(1,3)  = AF/Ao^  for  seaward  mooring  line. 

DELT(2,1)  = AF/Aa^  for  seaward  mooring  line 

DELT(2,2)  = AF/Aa2  for  seaward  mooring  line. 

DELT(2,3)  = AF/ACI3  for  seaward  mooring  line. 

Note:  The  last  3 cards  (#7,  #8,  and  #9)  provide  the  information 

needed  for  the  dynamic  analysis.  If  it  is  desirable  to 
perform  calculations  varying  the  data,  these  cards  may  be 
repeated  with  different  input  data.  There  is  a limit  of 
25  different  sets  of  data.  In  the  example  data  shown  in 
Table  D-l,  there  are  3 different  conditions  used. 


5.  Mathematical  Procedures  and  Program  Limitations. 


The  mathematics  has  been  described  in  the  report  and  Appendix  C. 

The  main  limitations  are  that  at  most  23  offset  points  may  be  used 
to  describe  the  shape.  This  has  been  found  to  be  very  adequate  for  the 
configurations  considered  thus  far.  Little  change  in  the  results  occurs 
when  more  than  15  points  are  used.  Computer  time  increases  about  as 
the  square  of  the  number  of  points. 

A listing  of  the  program  is  given  in  Table  D-2. 

6.  Flow  Qiart. 

A flow  chart  is  given  in  a figure  of  this  appendix. 
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8012/0* 

NU21/0* 

8U23/(08*S1 

8U12/(08*S) 

1*880*12/00 

1*880*21/00 

1*880*23/(00*8) 

1*880*32/(00*8) 

FY/OF 

GEM  8Y  HE»VE/HS*V£ 
tfcFLECTFO  BY  FXD  B0Y/cT* 
REFLECTED  ♦ INCIOEMT/ET* 

D13ENSI0NAL  FREOUENCY  - H? 
0*8* IMG  00  • AREA*»1E*W 


8U11/Q8 
8U13/ I 08*8 ) 

81)22/08 
8U31 / ( 08*8 ) 

8033/(08*8*8) 

1*880*11/00 
1*830*11/(00*8) 

1*880*22/00 
1*880*31/(00*8) 

1*880*33/(00*8*11 
EX/OF 

8Z/( 3F*B) 

GEN  BY  SWA  Y/SW*Y 
GEN  ttY  R0LL/R0LL(RAD)*8 
INC 10ENT/E  TA 
TRANS  BY  FXO  BOY  /El* 

8E48/VAVELENGTH 
AOOEO  8*SS  08  • *RE**»UE 
NAVE  FOTCES  UF«AREA»R0E*FT»*w2PHASE  rel  to  et*  *t  xo  - o*r, 
RhaSE  «EL  TO  BODY  80TIUN  - OEGNAVE  FIELD  - A8PL ITtJOF  TATTift 


POSITION  - X/WAVELENGTH 
Sn«Y  AHPlITUDE/ET* 

ROLL  A8RLITU0E IR*D)*B/ET* 

GEN  SY  RESULTANT  $W«Y/ET* 

Gfc N BY  resultant  lOLL/ETA 
TOTAL  REFlECTEO/ET* 

NOTION  RESPONSE 

OAK  HARBOR  BREAKWATER  - CORPS  IF  ENGINEERS  TESTS 


DIMENSIONAL  POSITION  - * 
HEAVE  A8PLITU0E /ETA 


GEN  BY  RESULTANT  HF  AVr /c  TA 
TOT*L  TRANS8ITTE0/FT* 


17  NAY  1073 


23 

0 10 

0 0 t 1? 

12.6 

10 • 0 0.0  l 

.9905  32.PFIJU  8'*8 

.1  .130290  .180  .216311 

.250  .280  .31729*  .371  .629 

-5.0 

G.O 

-3.0 

-1.23 

-3.il 

-2.30 

-3.0 

-3.75 

-3.0 

•:.co 

-A. 381 

-5.00 

-A.5BI 

-3.73 

-1.223 

-3.75 

-3.223 

-2 . 30 

-3.223 

-1.23 

-6.381 

-1.23 

-6.331 

0.00 

6.383 

O.Co 

6.5si 

-1.23 

1.223 

-1.25 

3.221 

-2.30 

3.223 

-3.75 

6.183 

-1.75 

6.383 

-5.00 

3.0 

-3.00 

3.0 

-3.75 

3.0 

-2.50 

1.0 

-1.25 

3.0 

0.0 

o.o  o 

•0  0.0  0.0 

66.5  0.0  3.0  0.0  1165. 

12.6 

10.0  0.0 

-2.36  25.1  6P1.  0.  9. 

9. 

0.  0.  0. 

0.  0.  9.  0.  0. 

12.6 

1C.0  0.0 

-2.16  23.1  671.  0.  1. 

116. 8 -3.26 

166.2  -3.732  13. 

21  -3.372  139.9  P.961  2*1.3 

-1376.  610. 

6 -1607.  1172. 

280.9  1713. 

12.6 

10. G 0.0 

-2.36  23.1  621.  1.  1. 

118.8  -5.26 

166.2  -3.732  10. 

21  -1.372  139.9  2.06A  261.3 

-1376.  -.10. 

6 -1607.  1172. 

280.9  1713. 

0. 

0. 


Table  D-l . 


Example  input  for  program  BRK2D  (Oak  Harbor 
breakwater) . 
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8RK20HU  /////////////////////////////////////////////////////////////////////# 


BUM  VERSION  FEE  79  3 17*1?  09/23/76 


PROGRAM  bR<?0  I INPUT.OUTPUT.pUNCh.TAPFB-TNP'JT.TABPE-OUTPUT) 


C***IATL5T  REVISION  ♦**♦*  27  AUGUST  147B 

C ***PRCIGR AM  BRK20  COMPUTES  THE  FIRST-DROE*  R*SPONS=  OF  AN  OSCILLATING 
C CYLINDER  ON  OR  NFAR  THE  FREE  SURFACE  of  AN  TOEAL  fluid  of 
C INFINITE  DEPTH 

3 COMMON  RI12<2S.25),<UJ6<25,25)»  POT<25.»5).  HOWI  25, 25 ).  =F  C ?«.*■) , 

1FK25.6),  RII25.25).  RJI25.25).  R<<?5,9),  »L(?5.9), 

2RMU ( 3* 3* 10) » RLAM(3.3*lO)»  F1I3.10).  0*L e R ( ». 10) . HW* ( 2* , 6 , 10) . 

3 DELHI 25.6. 10)  , XOL (25,10) 

3 COMMON / ONE /X( 25I»Y(25),XR(25)»YB(?5)»ANG<?5)»DFL<25)»VV(?5I 

1.FEINI2S),  eTTN( »5) , rni)RN<25,3)»  JC<5) 

3 C0MM0N/0NF2/CC3(25)»SS3(25) 

3 COMMON  /TWO/  N.NNW,  NaAVEL.  ISYM.  TSKT®.  NO.  P I E . GAMNA  • M , Tw  , T ® 

3 COMMON /THRt E / WAVF.LdO),  NNIIO).  • IKl’n.Tl 

3 COMMON/SI X/XN(5)»CN(5) 

3 COMMON  /SEVEN/  AREA.  R,  0.  ROE.  GFF.  BTTUEO).  TI TLE ( s ) 

3 COMMON  / EIGHT/LBLMUI3.3.3).  LBLAMI 3, 3. 3 ) , LB L F 8 ( 3,  3 > . LBL HV« < 7,  3 > , 

1LBL (10.3).  DE  G ( 3.10) 

3 COMMON  /NINE/  L8LRARI3.3).  LBLHWR('.I).  1*1*15.3) 

3 COMMON/TEN/0ELT(2.3).F3RI2,lO),PHASI?.l0).FO3N0(2.1O).PH*TDC2,?O) 

3 REAL  K 

3 DATA  XN/.26356O3H71S.A.913AO3U59107, 3. 59992877109), 

1 7. 0359 ICOO 5659.12.690000899 ?TR/ 

3 DATA  CN/. 521755610561, .396666811083. .0759929996917. 

1 .00361175367972.. 00002336997339/ 

3 PlE«ATAN2(0..-2. ) 

7 TP-2.*PIE 

11  GAMMA-0. 57721566 

C***BFGIN  READING  INPUT  DATA  AND  PRINTING  FC-H  CHECK 

12  3000  FORMAT  (641)) 

C*****KEAD  LABLES  FOR  PRINT  OUT 

12  PEAD  3000.  I ( (LBLMUII, J.L),  L • 1.3).  J « 1.3),  I • 1.3) 

36  READ  3000.  ( < < L6L 91 < I . J » L > . L • 1.3).  J • 1.3).  T « 1,3) 

63  READ  3090,  I (LBLFB(J.L),  L * 1.3).  J • 1.3) 

103  READ  3000.  1 1 L8LHWB (J. L ) . L • 1.3),  J • 1.7) 

123  READ  30U J , ((  LBLU.l).  L • 1,3).  J > l.»0> 

193  READ  3000.  1 1 LBLRAR I J.  L > , L • 1.3),  .)  > 1,3) 

163  READ  1000,  (ILBLHnRU.L).  L • 1,3).  J • 1,5) 

293  READ  3000,  (LBLRI 1»L  >»  L > 1.3) 

220  READ  20,  TITLE 

226  20  FORMAT  I6A10) 

226  PRINT  30.  TITLE 

239  30  FORMAT  IlHl,  8410///) 

219  READ  50,  N,  NW,  NwAVEl.  ISYM,  ISKTP.  1C. 

256  50  FORMAT  <5110,  616) 

C N • NUMBER  OF  8IRA1GHT  LINF  SFGNENTS  TO  8f  USEO  TO  FIT 

C THE  HULL NOTE.  THERE  MUST  3F  N*t  OFFSET  »OTNTS 

C NW  • NUMBER  OF  POINTS  IN  FRtF  SUBFILE  WHERE  WAVE  HEIGHT  *S 

C TO  Be  COMPUTED.  THIS  I?  TN  ADDITION  TO  THE  COMPUTATION  OR 

C WAVE  HEIGHT  9.9  WAVELENGTHS  ON  FTTHFB  SIDE  OF  THE  *!)DY 

C WHICH  IS  PERFORMED  AUTOMATICALLY 

C NWAVEL  • NUMBER  OF  WAVELENGTHS  AT  WHICH  COMPUTATIONS  ABF  TO 

C BE  PERFORMED 

C ISYM  . I FOR  SYMMETRIC  SECTION 

C ■ ANYTHING  FLSE  FOR  NON-«YMWeTBIC  SECTION 


Table  D-2.  Listing  of  program  BRK2D 


RUNT  VERSION  FEB  76  B I7«12  36/23/76 


256 

2eO 

262 

266 

267 

273 

Hi 


Hi 

362 


6U 


73 


362 

467 

BO 

367 

472 

606 

100 

c 

606 

613 

600 

c 

617 

633 

113 

c**«* 

633 

635 

637 

662 

666 

653 

655 

656 

6t  5 

676 

116 

5o3 

501 

505 

507 

516 

525 

112 

ISKIP  ■ 1 on  NUT  SOLVE  EoUATTDNS  3C  -9TI0N 

■ 2 Oil  NOT  SUL Vt  POTENTIAL  PROBLEM  (READ  I*  TOFFS) 

• ANYTHING  ELSE  $ JIVE  CnFFFfCTENTS  AND  DYNAMJCA 

MUMS  Ed  OF  BUOY  SEGMENTS  WHICH  ogp»c«eNT  FREE  SURFACE  RETWFFI 
CAT  Art  Aid  AN  HULS.  SEGMENT  N'JMBFRS  SPECIFIED  «Y  JC<5 

NKm  • NW  ♦ 2 
NC  ■ H - 1C 
Nwl  • 25  - n - > 

IF  (NW  ,LT.  01  NW  • 0 
IF  (NR  .GT.  NU1)  NW  ■ NW1 

PRINT  60.  N.  N.,  NR  A VC L>  ISVN,  ISR*»»  LC»  JC 
FORMAT  (13XPNU«BER  QF  SEGMENTS  •*.  T6// 

10X6NUN8ER  JF  FUE-SURFACf  STATIONS  • ♦.  167/ 

ICA6NUM8ER  IF  WAVELENGTHS  T6// 

10X*1$YM  «*»  16//  10X6ISKIP  •*.  T6 // 

10X6NC  • 615.  *JC  •*  515  /) 

RtAO  7C . AREA.  I),  0.  R13E.  GEE.  (BTTHFm,  T . 1,3) 

FORMAT  (SFlo.j,  3A10) 

AREA  ■ CROSScC  r IONAL  AREA  OF  tnnEPSFH  B10Y 
B * CHARACTERISTIC  LFN3TH  AS  S®ECIPIFO  BY  PTITLE 
0 • OlSTANCfc  BENEATH  SURFACE  1*  ORTGtn  if  USFPS  COORDINATE 
SYSTFN  (♦>.  ALL  NJTIJnS  REF®® EO  TO  THAT  POINT  ANO  BODY 
SHAPE  SPECIFIER  IN  THAT  SYSTFN 
ROE  ■ FLUID  DENSITY 
GFE  • ACCELERATION  JF  GRAVITY 
PRINT  60,  AREA, 3,  (ITITLEd),  I • 1,3).  D,  ROE,  GEE 
FORMATUOX  ‘AREA  • ♦,  F10.3  //,  10*  *B  ■ *C10.3  .5*.  3A1G 

mox  *c  • ♦ fid. 3 //  io*  ‘Fluid  density  ••  fio.5// 

lwX* ACCELERATION  OF  GRAVITY  ■*,  F\ 0 .3 / ) 

IFIISKIP  .E3.  2)  GO  TO  403 
READ  ICO.  (3  0L  < 1 ) . I • 1.  NWAVeL> 

FORMAT  (13F8.5) 

30L  • «F AN/WAVFLENGTH  RATIO  FOR  COMPUTATIONS 
03  C3C  I • l.NWAVEL 
WAVeL(I)  • R/30L ( I ) 

WAVfcLII)  • DIMCHSIONAL  VAVELeNSTH  DC  INCIDENT  WAVES 
PRINT  UO.  (ROLII).  I • 1.  NkAVEl) 

FORMAT  (lOXPPEAN/WAVElENGTH  RATIO*  1e  TNC TDFNT WA VFS • / / ( 23 *10= l l . 6 
) I 

INITIALIZE  OUTPUT  VARIABLES 

00  113  IL  • 1.13 

WN(1L>  • 0.0 

BOL(IL)  • u.u 

DC  116  I • 1,3 

FB(I.IL)  • 0.3 

OELFbl  I,  IL  ) • w.i> 

DO  116  J • 1,3 
RMUIl,  J,  ID  • 0.0 
R L AM ( I , J,  IL)  • 0.0 
CONTINUE 

DO  112  I • 1,  26 
XPL1I.II)  • C.G 
DC  112  J • 1,6 
HW8 ( I, J , 1 L I • 0.3 
OELW(I,J,IL>  *0.0 
CONTINUE 


Table  D-2.  Continued 


I 


112 


HUN r VUSION  *ce  7*  B 1/112  U*/23/76 


531 

533 

53* 

5*1 

551 

jj3 

if-} 

55/ 

561 

603 


6J3 


606 


631 

533 

6*2 

6*7 

5*7 

651 

553 


671 

672 
67* 
705 
717 
725 
732 
7*0 
7*7 
756 
772 

1007 

1015 

1023 

1031 

103* 

10*1 

10*5 


10*5 

1113 

1113 


113  CONTINUE 

;,****C0MPUTE  1 B/wA  Vt  L > A *0  N0NDIMENSI3NAL  WAVF  NO. 

DO  115  11  • 1.  NVAVEl 
SUL  I IL  > « B/WAVEL1  IL1 
115  WN1ILI-  TP  * 9/  W A VEL I IL  I 
C •**S [ 1 0 IN  OFFSETS  OF  CYLINDER 
NLP  • N ♦ 1 
NfP®  • NU®  ♦ 2 
NTOP  ■ NUPP  ♦ NW  - 1 
Nl  ■ N ♦ 1 ♦ NW 

PEAO  130.  (SKI.  I).  3IIM),  I-l,  NUP1 
130  FORMAT  (2F10.5) 

C KI (1. 1 1. PI  12.  II  • OIPENSIONAL  X,  Y COORDINATE  OF  OFF 9 F T 

C POINTS,  RESPECTIVELY 

IF  (NR  .LT.  1)  G 3 TO  1 95 

C***R£AD  IN  AOOiriONAL  POINTS  ON  THE  FRET  SURFACE  WHERE  WAVE  HFTGMTS 
C ARE  TO  BE  CONPUTcO.  TWU  STORAGE  LOCATIONS  OUST  BE  LEFT  PLANK 

C FOR  THE  P0SITI04  * WAVELENGTHS  FROM  THE  BODY 

READ  130.  <Rl(i.l>.Rl(2»n»I*NUPP.NT0P) 

C Rill, II. RI12.II  ■ COORDINATE*  0®  POINTS  ON  FREE  SURFACE-  WMFRF 

C WAVE  HEIGHT  IS  TO  BE  COMPUTED.  THIS  t«  TRUE  FOR  I .GT.  N ♦ R 

CwwwnON-DIHENSICNALIZE  OFFSETS 
OC  lev.  I • NUPP, NTOP 
XIII  • RIll.n/B 

ifi  Yin  • -l.w£-C8 

IB5  CONTINUE 

00  1*0  I • l.NUP 
min  • Riu.n/B 
1>V  YB1I1  . RII5.II/B 

C...  COMPUTE  MIDPOINT. ANGLE  ANO  LENGTH  OF  STR A TGMT-L l NF  SEGMENT*. 
C«****ANO  COMPONENTS  OF  NOAAAL  TO  BODY 
JF  • 0 

00  200  J • 1 » N 
XI J)«0.5*IXBI JlTXBI J»l)  I 
Y( J)-O.S*f Y*l JI+Y31 J*1  ) I 
T1"YA( JTlI-YBl Jl 
T2-XA1JF1 l-XBl J) 

ANG 1 J I • ATAN2 1 T1.T2I 

CC3IJI-C0S1 ANG1JI I 

S*3( JI-SIN1ANG1 J 1 1 

VVl JI-XI J)*CC3I JI*Y(J)*SS3( J) 

CEUJ  I • SOR T 1 T2**2  ♦ T1**2I 
RNORMU  ,1)  • -SS31JI 
RNORMU  ,»l  ■ CC31JI 
RNORMU  ,31  • VVIJI 
230  CUNTINUE 

PRINT  33.  TITLE 
PRINT  250 

250  FORMAT  levx*c VLINOFR  C s OMF ER Y» / / r QXWDt MENS  TONAL  OFFSETS*. 

1 11X*N0N-D1MENSMN»L  OFFSETS*.  SX*« T non tm t s OF  SEGMENTS*// 

2 6»*l*.  l6X*X*.  9X*Y*»  1HX*X*.  ,»X*Y*,  10X*X*,  *X*Y*. 

3 18X*SL  0PE*.  * X*LE  NOTH* / 1 

PRINT  270,  II,RI(l,II,RII2,II.Xdin.YP(TI,y|n,YII),ANG(II, 

1 DELI  II , 1*1. Nl 

270  FURMAT  IX,  16.  AI10X2FIU.3I I 
NL  ■ N ♦ 1 
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1115 

PRINT  270.  N(.»  RIU.nl).  Rtd.Nll.  y<unl>.  YMNLI 

1H3 

PFINI  231 

1117 

281 

FORMAT  ( / / 1 0 X » *2051  T IONS  FOR  HAVF  HFTG-TT  C AL  CUL  A T 1 ON  S * / > 

il«i7 

280 

FORMAT  <//> 

US  7 

IF  (NW  .IT.  1)  G)  TO  290 

119? 

PRIN1  271.  (I.RI(l.Il.RI(2,n#X(I>,Y(tl.!*N'l»R,NT0Pl 

1201 

271 

FORMAT!  X.  16,  2U0X,  2F10.3I) 

1201 

PRINT  280 

1710 

2°0 

M-N*Nn*2 

c«**«*TRANSFER  TO  COORDINATE  SYSTEM  IN  eo'F  Rlipc^r 

1213 

00  2P5  I • l.M 

.2  IS 

r b < i » * r « ( 1 1 - o/8 

UZZ 

285 

Till  ■ Y< II  - 0/R 

1233 

TP ( NUP ) • Y3(NU®>  - 0/8 

C COMPUTE  FACTORS  OF  I AND  K INDE PEND ANT  OF  FSFO'IFNC* 

IZsZ 

CALL  COcFF 

1213 

Y?URF  • -1.01-03  * 3 

c • • • 

START  FREQUENCY  ITtRATION. 

1215 

00  301  1L  • 1.  NMAVEL 

1217 

K • WNI  IL  1 

c 

ALL  POTENTIALS  INITIALIZED  TO  ZERO. 

C • *** • F t < I.J)  • NONOIMENSIONAL  AM?L 1 TUOe  0=  POTfNTTAL  AT  POINT  T 

DUE  Tn 

C****«MOD£  J.  (ASSOSIATED  -ITT  COS(wT)  ».  FTII.J1  IS  SIMILAR  TO 

FF(I, J) 

C«**®*0|.T  ASSOSIATED  WITH  SlN(WT).  J ■ 1,?,  3. 1.5,1  I«PLY  RESPECTIVELY 

c+***< 

*S»AY . HEAVE.  ROLL.  DIFRACTED.  INCIDENT  ANN  OVERACTED  ♦ INCIDENT 

1252 

00  1 1-1.25 

1253 

DO  1 J • 1.6 

1251 

Ft  d.  J»-6. 

1261 

FKI.JI-O. 

1265 

1 

CONTINUE 

C***AOD  POINTS  TO  THE  OFFSET  ARRAY  FOUR  HA VF l c NGTHS  FROM  T3E  ORIGIN 

C 

ON  THE  FREE  S'J»FACF 

1271 

XIN.ll  • 1.0  * HAVEL ( J L ) / b 

1300 

XIN.2)  - -1.0  ♦ WAVfUILl/3 

1307 

YIN-1)  . YSUPF 

1312 

V ( N *2 ) • YSUPF 

C***COPPUTt  INCIDENT  WIVE  POTENTIALS  AND  normal  VELOCITUS 

1316 

DO  1C2  1 • l.M 

1317 

D01U1  IC  • 1.5 

1320 

101 

I F ( I • E 0.  JCIIC) > GO  TO  102 

1326 

EY  • EXP(K*Y( I) ) 

1335 

CKX  • COS ( R* X ( I ) ) 

1311 

SRX  • SINLRAXlIll 

1353 

I F ( I ,CT.  N)  GO  T1  1j3 

1356 

FEINU)  - £Y«(SS3(  I)*SKX  ♦ CC3d)*'‘«X) 

1370 

F1INII1  --EY«( SS3( l >«C*X  - CC3  C II  * T<  * 1 

1103 

103 

FEtI.51  • EY*CKX*(1./KI 

1111 

Fid, 5)  • £Y*5kX*(l./x) 

1121 

102 

CONTINUt 

1127 

A K • K/B 

1131 

WRF  • SOR T ( GE  E *4K ) 

1136 

WF  • HR F / TP 

1113 

» T • 1.0/HF 

1111 

PRINT  30c  • K.  hRF,  «F . WT,  HAVEL ( TL ) 

1160 

300 

FORMAT  «//.  ♦ -AVE  NUMBER  • K •*.  '9.3,  *Y*CTRC"LAP  FRfO'JENCY  • *, 

1 F9.5.  '/‘FREQUENCY  F9.8,  3Y*®FRno  FK.5. 
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l*6u 

C 

TK*2.J/K 

FJkST-ORDER  POTENTIALS  on  CYLINDER  ARE  civet  CAlCULATEO. 

1*62 

Cat* 

CALL  CONPIK) 

FIAST-ORDFR  PHYSICAL  CuANTITIES  ARF  CALCULATcp. 

1*6% 

CALL  PMYSCLOO 

1*66 

301 

CONTINUE 

1*71 

1009 

FCRHAT  (1CFB.6I 

1*71 

1SKIP1  • 2 

C6****PUMCh  RESULTS  OF  POTENTIAL  SOLUTION  on  CARPS 

1*72 

PUNCH  2(  » TITLE 

1*00 

PUNCH  60.  N.  NW.  N * -*.  V f L . ISYN.  ISRYPl.  tC.  JC 

1622 

PUNCH  To#  AREA.  H , D » ROE . &FE.  »TTTLE 

16*2 

PUNCH  1009.  (VNULI.il  • l.lu),  (BTL ( TL ' » U • l.ltl 

166* 

PO  310  I * 1.3 

1666 

PUNCH  lul/9  . (F6d.IL).  IL  • l.ll).  ( 9eL  F*  ( T.  IL ) . U • 1 

1613 

00  310  J • 1.3 

161* 

110 

PUNCH  lu09»  (RMU< I. J. IL ) . 1L*  1.10),  ( 9 L A* ( T . J , I L ) . U ■ 

1662 

DO  320  I • l.NNV 

166* 

PUNCH  1009,  (XOL(l.lL).lL  • 1.10) 

1673 

CO  320  J * 1.6 

1672 

320 

PUNCH  10C9.  (HV8(  I, J.  IL  ) . IL  ■ 1.131  » ( PE l Of  I . J » It  > » II 

1710 

303 

CONTINUE 

1730 

116 

CONTINUE 

1730 

CALL  DYNAHC 

1731 

302 

CONTINU* 

1731 

STOP 

1733 

END 
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2 


z 

z 

2 

2 


2 

2 

2 

2 

5 

27 

27 

31 

56 

60 

115 

117 

113 

135 

173 

173 

ns 

177 

217 

221 

253 

253 

256 

260 

260 

262 


272 

312 


312 

116 

117 


SUBROUTINE  DYNAMC 

COMMON  RI12I25.25)»RK:6(25,25I»  PPT(25.»5>.  -tOWI 25, 25 ) . FF ( 26, 6) , 
1F1J25.6I.  RK25.25).  RJI25.25).  *K»>5,A).  RLI25.AI, 

2RNUO.3.10).  RLANI3. 3»lu>,  F9<3,1D>,  DELF»<3.10),  HWB ( 25»6 . 10) , 

3 CELWl 25.6.10)  , XOl (25,101 

COMMON /ONE  /X(25I,Y(25)>  X8 ( 25 ) » VBJ 2« > » 4NS t 75 ) , OEL ( 25 ) » VV « 25 ) 
1,FE1N(25>,  fiin(25>,  RNJRM(25,3>.  JT  t *5 1 
COMMON  /T'mO/  N.NN4,  NW4VEL,  ISYM,  TSMI»,  NC  , P IE , GANNA, H, TK , TP 
COMMON/ THRE E 7 NAVELUO),  WN(10),  80Lt  ‘ . U 

C OMMON/ FOUR /R AR (3,1U),0ELR(3,10),  HWR (76,3,101,  OELWR < 25, 3, 10 >, 

2 HWT<2S»10>.DE12T(25,1U)»RKHVO<3.R».  RKMOR (3,3),  RKT813.3). 

3 XG.  YG,  RMASS,  RlNtRT,  0AMP<3> 

COMMON  /SEVEN/  AREA,  8,  0,  ROE,  GCF,  ATTTLEl 3),  TITLE  t R ) 
COMMON/TeN/nELT(2,3),FORt2,10),PNAR(?,lO),eORNOt2,lGI,»HA«OI2,10) 
DIMENSION  A ( 6, 6 ) , C(6),  ERASE(6) 

IF  < ISMP  .NE.  2>  GO  TO  100 
C*****RtAD  POTENTIAL  CJEFS  IF  ISKIP  ■ 2 

READ  1009,  (WNUU.Il  • 1,10),  inilTU,  TL  • 1,10) 

1009  FORMAT ( 10FB, 3 ) 

DO  110  I • 1,3 

READ  1009,  (FbU.ILI,  IL  • 1,10),  nFlFMT.TM*  IL  » 1.101 
DO  110  J • 1,3 

110  READ  1009,  (RNU( I. J, IL  I.  IL-  1, 10),  (RLAK  T,  J.  IL  >,  IL  • 1,10) 

DO  120  I • 1 . NNW 

READ  1009, I XOL I I , IL),  IL  • 1,10 

DO  120  J • 1.6 

*20  READ  U39,  (Hwbt  I,  J,  IL).  IL  • 1.10  , ( OeL  W(  I • J , IL  ) . IL  • t.TD) 
100  CONTINUE 

C *****QL'TPUT  POTENTIAL  COEFS 
call  POTOUT 

I F ( ISKIP  .£0.  1)  GO  TO  1 AC 
C«****REAO  OIME NS 10NAL  JYORPSTATIC  SPRING  CONSTANTS 
READ  luca, ( (RRMYOl I,J) » J- 1 » 3 ) ,1-1,3) 

C*****STAAT  LOOPING  THROUGH  OIFFERENT  DYNAMIC  CONFIGURATIONS 
DO  1 AO  K1  ■ l»5v 

READ  1010,  AREA.  R,  XG,  YG,  RMASS.  » T NFR T. D AMP ( 1 ) , DA*® ( 2 ) . D AMP d ) , 
1 NPUNCH 

,01w  FORMAT  (6F10.3.  3F5.2,  16) 

IF  (EOF, 5)  121.  122 

121  STOP 

122  CONTINUE 

OU  5 K2  • 1,3 


5 IF! DAMP  I K ?)  ,E0.  -C.O)  DAMP  ( K2 ) • D.D 

C XG, YG  ■ COORDINATES  OF  THE  CENTER  OF  GRAVITY  OF  THE  BODY 

C NOTE.  MOMENTS  AND  MO-EMT*  OF  TSERTIA  ARE  COMPUTED 

C ABOUT  THE  CENTER  OF  GRAVITY 

C 04MP  AODS  CORRECTION  FOR  VISCOUS  0®  NONLINEAR  DAMPING 


c*»***KtAO  DIMENSIONAL  MOORING  SPRING  CONSTANTS 
READ  1008, << RUMOR ( I , J ) , J • 1 , 3 ) , I • 1, 3 ) 

1008  FORMAT  (9F6.3) 

C**P*RNONDlM£NTIONALIZt  SPRING  CONSTANTS  ,MAS>,  MOMNNT  OF  INERTIA 
C *****ANO  CG  COORDINATES 

0 • AREA*  ROE*  GEE/B 
00  130  1-1,3 
DO  130  J • 1.3 
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320 
331 
3*6 
3 72 
*02 
*06 
*12 
*13 

*15 

*16 

*21 

**5 

*71 

531 

552 

576 

623 

633 

6*3 

65? 

65* 

655 

666 

705 

723 

736 

753 

771 

773 


77* 

1005 

1006 
1032 
10*7 

1050 

1051 

1052 

105* 

1072 

1111 

1132 

1155 

1163 

120* 

1225 


I F ( RKMQR ( I,  J ) .60,  -O.0I  RKMOR (I » J > • 3.0 
RKTB(I»J)«  (RKMYDd.J)  ♦ RKMOR ( 1. J )) /O 
130  IFd.EO.  3 .OR.  J • EO  • 3)RKTBd.J>  »RKTR  ( I,  .))  /8 
PKTBI 3» 31  • RKT3(3»3)/B 
RMASSB  » RMASSM ARfcA*ROF > 

RINEFB  • RINERT/(AREA*ROE*B*0) 

*08  • *0/3 
YOB  • YG/8 

C*****ST  ART  WAVELENGTH  LOOP 
00  150  IL  • 1.10 
IF  (IL  .GT.  NWAVELI  GO  TO  160 

C*****SET  VALUES  IN  N JND INENT I ON AL  1 2E 0 ALGE  R*  A TC  FOUATIONS  OF  NOTION 
A 1 1. 1 ) • RKTBUfl)  - WMdL)  ♦ (RMASSB  ♦ RMU(l,l,IL>> 

A ( 2 . 2 1 ■ RKT8I2.2)  - WN(IL)  * (RMASSB  ♦ »WU(2,?,IL>) 

A ( 3. 3) • RKTB  ( 3,  3 ) -YGBARNASSB  -WN(TL)  ♦ (RTNER8  ♦ RNU(3.3.Tl)  ♦ 

1 ( XG8**2  .YGB**2>  • RMASSB) 

A ( 1 , 2 ) ■ *K  T8 ( 1 » 2 ) - VN(IL)  ♦ RMU(1.?»TL> 

A(lf3l  • RKTB  d » 3)  - WN (IL ) • (RM'J  ( l , 3,  T L > -YGB*RMASSB) 

A ( 2 . 3 ) • RKTB ( 2. 3 ) - WN(IL)  ♦ < R MU  ( ? » 3 » I L ) .*GB*Rma$$B) 

A(2.l)  • Ad.  21 
A ( 3* 1 ) • A ( 1.  3) 

A(3. 2)  • M2.3) 

DO  20  I • 1.3 

DC  10  J ■ 1.3 

A(1. 3.  J.3 1 ■ Ad.JI 

Ad,  J.3)  -RLANd.  J.IL)*S0RT(NN(IH) 

C ADO  CORRECTION  FOR  VISCOUS  DAMPING 

IF  (I  .EO.  J)  All, J.3)  ° <1.0  . OANP(I) )*A( I, J.3) 

10  A ( 1.3,  J ) ■ - Ad,  J.3) 

Cd)  • FSd.ILI*  SINIDELFBI  I,  IL)  I 
20  CII.3I  > PBd.It)  * C0S<DELF8(I»Il)) 

SCALE  ‘ 1. 

NN  ■ 6 

C*****SCL VE  ALGEBRAIC  EOS  OF  MOTION.  8M).  B(»).  B(3)  ■ ANPLITUOFS 

C ♦•♦•♦OF  COS(wT).  6 ( * ) . 3(51.  B ( 6 ) ■ AMPLITUDES  OF  StN(wT)  FOR  SWAY,  HEA 

C*****ANO  ROLL  AT  CENTER  FO  USERS  COORDINATE  SYSTEM. 

LL  • LNEQF<6. NN.l, A. C. SCALE. ERASE) 

DO  30  I • 1.3 

C*****AMPLlTUOE  ANC  PHASE  OF  RESPONSE 

RARd.IL)  • S0RT(C(I)**2  ♦ Cd.3)**2) 

30  OELRd.lL)  • AT  AN2  (C(I)»Cd*3)) 

OC  *0  I • 1,  NNW 
AW  • 0. 

BW  ■ 0. 

90  DO  50  J • 1.3 

C***P*RESULTANT  WAVE  AMPLITUDE  AND  PHASF  FOR  SWAY,  HEAVE  ROLL. 
HaRIUJ.ILI  • HWBd.J.IL)  * RAR(J.TL) 

OF  LWR  ( I,  J , IL  ) * DELWd.J.lLI  ♦ DELR(J.U) 

AW  > AW  . HWRII.J,  IL)  * SIN( DELW* ( T»  J. Tl ) ) 

50  BW  » BW  ♦ HWRd.J.IL)  • COS  ( OELWR  ( I > J • IL  ) > 

I F ( XOL ( I . IL  > . L T . 0.)  GO  TO  70 

C***p*TOTal  REFLeCTEO  WAVE  (VECTOR  ADOITIONI 

AW  • AW  . HwB ( I, 6, IL  ) * SIN( OELW( T, *. TL ) ) 

BW  • BW  . HW6( 1.6. IL)  * C3S( OELWd • 6*  Tl ) ) 

GO  TO  *5 

C****pTOTAL  TRANSMITTEO  WAVE ( VECTOR  ADDITION) 
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1227  7t  AM  • AN  ♦ MVeiI,6,lL)  * SINIOFLWd.R.ID  » 

1266  BW  • 8W  ♦ rlMB(I.6,IL)  * COS  I DEL*  ( T . 6 . IL  ) ) 

1267  65  HNT(I»IL)*$3RT(AW**2  ♦ 8W**2) 

1307  DEIMTI 1,  IL  ) "ATAN2  ( Am,  PM) 

1317  60  CONTINUE 

1321  GO  TO  150 

C*****SET  OUTPUTS  FOR  H .GT.  nmAVEL 

1322  160  DO  170  I • 1,3 

1326  DIRU.IU  • 0.3 

1330  170  DEL* I i, 11 ) > 0.0 

1337  00  180  I ■ 1,25 

1363  HWT 1 1, I L > • 0.3 

1366  CELWT (1,11)  • O.U 

1351  DO  180  J • 1.3 

1352  HxR (I,J, IL)  • 0.0 

1361  160  DEC  WR  ( I , J • II. ) ■ 0.0 

1376  150  CONTINUE 

C*****  OUT  PUT  DYNAMIC  RESUITS 
1376  CALL  OYNOUT 

1177  NMOR  ■ 0 

1600  DO  133  IP  ■ 1.3 

1602  DO  139  10  • 1,3 

1603  IF  IRKMORI IP. 10)  .EO.  0.0)  GO  TO  139 

1610  NHOR  • NMOR  ♦ 1 

1612  139  CONTINUE 

1616  IF  INMOP  .N6.  0)  CALL  MORTEN 

1621  IF  INPUNCH  ,NE.  0)  GO  TO  160 

1623  PUNCH  2000 

1627  20 00  FORMAT  1*1111111111*) 

1627  PUNCH  2005,  <60L< 10 1 . HWTI 1 . 10) ,R AR 1 1, TO) , » AP I 2 , 1 0 ) , R AR ( 3, T C) . 

1 10  > 1,10) 

1666  2005  FORMAT  I5F10.6) 

1666  PUNCH  2000 

1672  PUNCH  2010,  (50L(I0)»FJRNDIl»I0),r39N0(P,T0),I0»l,)0) 

1517  2010  FORMAT  I F 1 0. 6, 2E 20. 6 ) 

1517  160  CONTINUE 

1521  RETURN 

1522  END 
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2 


2 

2 

2 

2 8 


22 

26 

32 

34 

37 

37 

45 

53 

53 

53 

70 


73 

72 

103 

114 

126 

154 

203 

215 

225 

236 

246 

250 

306 


306 

310 

311 


SUBRCJTINE  HORTEN 

C***SUeROUTINE  HORTEN  COHojTcS  FORCES  IN  THE  HOOR  TNG  LINES 

COMMON/FOUA!RAR(3.10).OELR(3.1U).  HWR(?5.3.m.  QELWR ( 25, 3, TO  ) , 

2 HWT(25#lO).0clWT(25.10l,RKHY0(3,3>,  *»HO»(3,3l,  RKTB(3.3), 

3 Xt i»  YG » RMASS.  R iNt  R T*  0AMP(3) 

COHHON  /Se  VEN!  ARE*.  3.  0.  ROt.  GFF*  BTTTLe<3).  TITLE  1 3 > 
C0MN0N/TEN/0ELT(2,3).FdR(2,lDlf PMAS(2.l3t.FnRND(2,10>.BHA*0(2»10) 
READ  10.  ((DELT(I»J)»J-l*3)»I-l»2) 

10  FORMAT  (6F13.2) 

C OELTU.lt, 1-1,3  * CHANGE  IN  FORCE  IN  SHOR'WAPD  MOORING  LINF 

C PER  UNIT  DISPLACEMENT  TN  SWAY,  HEAVE  AND  ROLL 

C 0ELK2. 11.1-1*3  - CHANGE  IN  FORCE  TN  SCAVARO  MOORING  LINE 

C PER  UNIT  DISPLACEMENT  TN  SWAY.  HEAVE  ANO  ROLL 

CAB  • 1.0/<R0E*GEE*ARE4) 

CONS  • 183.0/AC0SI-1.GI 
00  100  J ■ 1.2 
PRINT  20 

20  FORMAT  ( !//!2QX*M008ING  LINE  MODEL  RESULTS*!! 

IF  (j  .tO.  1)  PRINT  IB 
IF  (J  .EO.  2)  PRINT  IS 

18  FORMAT  (30X*SHQREW4R0  MOORING  LINE*!) 

19  FORMAT  ( 30X*SEAWARD  MOORING  LINE*!! 

PRINT  30.  (DELT1 J.KI.R-1.3) 

30  FORMAT  (•  CHANGE  IN  FORCE  PER  UNIT  MISPLACEMENT  IN  SWAY.  HEAVE* 

i ♦ ano  roll,  respectively  ■*,  M'lo.A!/) 

C***COMPUTE  FORCES  IN  MOORING  LINES  ANO  PHASE 
00  50  I ■ l.lo 
AA  • RAR(i»I)*OELT<J.l> 

AB  • RAR(2.t!*0ELTlJ.2( 

AC  ■ RAR(3.i:*0ELT( J.31/B 

IS  " AA*STN(DELR(1.I!)  ♦ A6*SINt DEI R C 2. T ) ) ♦ AC *S IN ( OE LR ( 3 . I ) ) 

TC  ■ AA*C0S(DEL»I1.I!)  ♦ A8*C0S ( DEL» C 2. t ) I ♦ AC*COS ( DE LR ( 3. T ) ) 
FORU.JI  • SOR  I ( TS*  I S ♦ TC*TC  ) 

PHASU.I)  • ATAN21TC.TS) 

FORNDU.I)  • CAB*FOR(  J.  I ) 

PHASOIJ.il  • CONS*PHAS  I J.  I I 
50  CONTINUE 
C ***PR INT  RESULTS 

PRINT  BO.  (FOR(J.n.I-1.10).(PHA$6(  J.T1.T-1.10). 

1 (FORNDU.I). 1-1*10) 

83  FORMAT  ( 3X*M0QR  I ,MG  LINE  RE  S PONSE  */ 5X»  FOR  CE  AMPl  ITUOE  !ETA*»  ll». 

1 10E10.3/!X*i»HASE  REL  TO  ETA  AT  KO  - OeO  *.  1CE10.4!! 

2 5X30HF0RCL  AMPLITJ0E/R3E*G*AREA*ET*.  10E10.3) 

100  CONTINUE 

RETURN 

END 
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2 


2 


2 

j 

2 

6 

10 

1) 

26 

34 

41 

45 

53 

66 

101 

105 

111 

112 

120 

124 

131 

144 

157 

163 


167 

201 

214 

22 

232 

235 

237 

241 

243 

245 


307 


351 

353 

354 

356 

357 

361 

362 
364 


SUBROUTINE  COEFF 

C THIS  SUBROUTINE  CALCULATES  THE  PART?  0*  IIT.J)  AND  K(I,J) 

C WHICH  ARE  INDEPENDENT  OF  FREQUENCY  N'lHBER  K. 

CriHMON  RI12<  25,25)  , RK 1 6 ( 2 5 » 2 5 > , P0T(»«,»5),  HOW(  2 5.  25  I . F E I 25. 6 ) . 
1 F I ( 25 , 6 > , R I (25,25),  RJ(25»25)»  »K(25.4),  01(25,4), 

2R  HU ( 3, 3 , 10  I , RLAMI 3.3,10),  FB(3,10>,  0ELFR(3,10),  HwB (25,6,10), 

3 Of LW(25.4, 10)  , XOL (25, 10) 

C0HH0N/0NE/X(2SI. Y(?5>,XB(25l,Y&(?B),4MG(»5t,0EL(25)  ,VV(25) 
1,FEIN(25),  FIIN(25),  RNtlRHI  25, 3 ) , JC(5) 

COMMON  /TwQ/  N, NNW,  NWAVEL,  ISYM,  TSXIP,  NC,  P IE , GANN A, M, TR, TP 
COMMON/ONE2/CC3(25),SS3(25) 

N2  • N/2 

00  A I * 1,  M 

IF ( I .GT.  N)  GO  TO  7 

IF(ISYM  .ED.  1 .ANO.  I ,GT.  N2)  GO  TO  7 
K 11  * XIII  - XB  ( 1 1 
Yll-Y) I )-YB(l ) 

X21  • XU  ♦ XR(1) 

Y21-YII)»Y3(1) 

PPl-ALOGI X11**24Y11**2) 

P01»AL0G( X11**24Y21**2) 

TP1*ATAN2(Y11»X11) 

T0I-ATAN2( Y21 , Dll ) 
od  1 J ■ 1,  N 
xi2»x( i >— xa ( j*ii 
Y12*Y(I)-YB(J»1) 

Y22  *Y( I ) »YB 1 J*1 ) 

PP2*AL0G<  X12**2»Y12**2) 

P02»AL0G(X12**2*Y22**2) 

TP2«ATAN2( Y12.X12) 

T02-ATAN21 Y22.X12) 

C CORRECTION  FOR  DISCONTINUITY  IN  ATAN2  AT  PTF 
I F ( XI 1 .GT.  0.  .OR.  X12  .GT.  0.)  GO  TO  A 
IFITP2  .GT.  0.  .AND.  TP1  .LT.  0.)  T»t  • TP1  ♦ TP 
IF( TP2  .LT.  0.  .AND.  TP1  .GT.  U. ) TP1  . TP1  - TP 
a C 3 • CC 3 ( J ) 

S3*  SS3( J ) 

Al-PIE 

IF ( I- J ) 2, 3, 2 

2 A1-TP1-TP2 

3 A2* T02-TC 1 

A5*C3*(-XR( J* 1 ) ♦ <8 ( J)-X12*0.5*PP2*XH *0. 5*»»1 
1*Y12*TP2-Y11*TP1)*S3*( YBU)-YB( JtT )-Xl2*TP2 
1-Y12*0.5*PP2+X11*TP1*Y11*0.5*PP1) 

A6*C3M-XB<  J*1)TXB(  J)-X12*0.5*P02  + XU*0.  5P»0I 
1*Y22*T02-Y2J*T01)-S3*(-YB( J)+YB( J*l)-X1»*T02 
1-Y22*0.5*P02*X11*T01*Y21*C.5*P01) 

4 X11-X12 
Y11-Y12 
Y21-Y22 
PP1*PP2 
P0l*P02 
TP1-TP2 
T01-T02 

K 1 12 ( I , J ) * A1  - A2 
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CQEFF 


•unr  ygrsic*  fee  7*  * i5»i*  03/i4/7e 

372 

1»J)  • A*  - AC 

17  7 

GJ  TO  1 

<*oo 

7 

00  200  L • i*N 

*02 

*11211*11  • C.C 

407 

200 

««56  U0D  • 0.0 

416 

1 

CQNTlNJt 

423 

RtTUB** 

424 

|N0 

Table  D- 
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SUBROUTINE  COMP(K) 

C THIS  SUBROUTINE  COMPUTES  THE  CO?Fr TC I ENTS  TfpENOENT  ON  K 

C AND  CALLS  ON  LNEOF  TU  SOLVE  THE  SIMULTANEOUS  EQUATIONS 

C FOR  THE  VELOCITY  POTEMIAtS  F£  ( ,JJ\  ANT  Et(r,JJ).Fne  FE ? 

6 COMMON  R I 1 2 ( 25 « 2 3 ) » RR5b ( 25 » 25 ) • PnY(?«.»?».  HOW ( ? 5 , 25 ) . FE ( 25, A ) , 

1FK25,6I,  *1(25.25).  KJ(25.25>,  K(»5.4),  H(25»4>, 

2RMU(3,3.10),  RLAMI3.3.10I,  *6(3,101.  0FLpM(3.101»  HWB t 2 5 , 6 , 10  ) . 

3 OELW< 25,5,10)  , X0L(25,10> 

6 COMMON /ONE / X ( 25 ) • Y ( 25 ) , XB 1 25  > , Vt ( 75 ) , ANG ( 25  t , OF L ( 26 ) , VV ( 25 ) 

l,Ft INI25) * F l IN ( 25 1 . PNORM ( 25  » 3 ) • JC C M » 

6 CUMM0N/0NE2/CC3I25) 3) 26) 

6 COMMON  /TaO/  N.NNW,  N a A VE  L » ISYM.  TX<T».  N r,  PIE, gamma, M*TK*T0 

6 DIMENSION  A(50.50>,9('U,4>,ERASE(5T> 

6 REAL  K 

6 N2  ■ N/2 

12  00  1 I • 1 , M 

13  DO  6 12*  1*4 

14  RM  I , I 2 ) • 0.0 

21  6 RL ( 1, 1 2 > • U.O 

27  DO  4 IC  ■ 1,5 

31  4 IF  ( I .EO.  JC(IC)  I 60  TO  1 

37  I F ( ISYM  .NE.  1)  60  TO  3 

42  I F ( I . 6T.  N?  .AND.  I ,LE.  N)  CU  TO  0 

55  B Xll  • X ( I ) - X5(l> 

63  X21  * Xll  ♦XBll) 

66  Y21-Y( I)*Y3(1> 

74  P01-ALOGI Xll**2*Y21**i ) 

107  T01»ATAN2(Y21»Xil) 

113  CALL  CPV(Xll,Y21,E2l.Cll»Sll.A911,AlOH,t<) 

124  C21-C11 

126  S2i*SU 

127  A921* A911 

131  A102 1* A 101 1 

132  DO  7 J ■ 1,  N 

135  X12*X( 1)-XB( J + 1 » 

143  Y22* Y( I ) ♦ YB ( J*1 ) 

147  PQ2-ALOG(X12**2*Y22**2) 

163  TQ2*ATAN2(Y22,X12) 

167  S3«S53(J) 

172  C3*  CC3 ( J ) 

175  CALL  CPV«X12,Y22,E22.C12,S12,A912,A1T1?,K) 

206  DO  13  IC  • 1,5 

211  13  I F ( J .EQ.  JC(ICI)  GO  TO  41 

217  A3.A1O11-A1012 

221  A4«E2l*Sll-622*S12 

225  A7-S3* (0.5*(PQ1-P02I *4912-4911) ♦C3*TT0l-T0?*41011-AI012> 

243  A8*E21*SIN(K*X11-AN6( J) )-E 22* S IN ( K *X l 7-ANG ( J) I 

266  5 R I ( I, J ) • 2.  ♦ A3  ♦ *112(1, J) 

300  IF  (I  .NE.  J ) GO  TO  3 

303  R I ( I , J ) • * I ( I , J ) - TP 

314  3 RJII.J)  • -TP*A4 

322  POT ( I, J ) * TK*A7  ♦ RK56(I.J> 

334  HOW (1,1)  • — TK*P I E*Ao 

342  CO  10  L • 1,3 

344  RK(I,L)  ■ RKd.LI  ♦ POKI.J)  * *N')*N<J,H 
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364 

10 

VllliLI  • RLIIfU  ♦ HOW(I.J)  ♦ RNORN(J.l) 

<•10 

RKII,*)  • RK  ( I.  * 1 - FEINI J)*PQT(I. J)  ♦ 

c I IN  ( J ) *HQif  ( I » J ) 

*36 

RUI#*)  • HlIiO  - FE  INI  J)*H0W(I#  J)  - 

FITNIJ)*POTI I. J) 

♦ 63 

*1 

IF(J-N)2,7.7 

*66 

2 

711*712 

*70 

T21*T22 

*71 

P0l*P02 

*73 

101*702 

*7* 

A9U-A912 

*76 

71011 *11012 

*77 

C11*C12 

501 

Sll*5>12 

507 

E21 *122 

50* 

7 

CONTINUE 

507 

GO  TO  1 

5v7 

9 

IS  • N - I ♦ 1 

512 

HO  12  L • 1#  3 

513 

RMI.l)  * RK(IStl)  ♦ <-1.0>**l 

527 

12 

RLU.Ll  • RL<IS.LI*<-1.0)**L 

5*5 

DO  11  J ■ It N 

5*6 

DO  16  IC  • If? 

5*7 

16 

IF  ( J .EO.  JCIICM  GO  TO  11 

555 

JS  • N - J ♦ 1 

560 

RKIfJ)  • RKISf  JS) 

570 

RJ(lfJ)  • RJIIS.JS) 

601 

RK ( I f * 1 « RK ( I f * ) - F£IN( JIPPOTII*. JS) 

* CTTNI  JMHOWIIS.  JS) 

627 

Rltlf*)  • R L ( I f * ) - FtINU)*HOW(I*f  JS) 

- fttni J)*PCT( IS. JS) 

655 

11 

CONTINUE 

66) 

1 

CONTINUE 

66) 

12  ■ 0 

66* 

32 

CO  22  I ■ 1 f N 

666 

DO  1*  1C  • If  5 

667 

1* 

IF  l L .EO.  JCIICM  GO  TO  22 

67  5 

12  • 12  ♦ 1 

677 

II  • 12  ♦ NC 

701 

00  31  C • If* 

702 

MI2fl>  • RK(I.L) 

712 

31 

B(Il.l)  • Rl(l.L) 

725 

J2  • 0 

726 

DC  22  J ■ 1 f N 

730 

DO  15  IC  • 1.5 

731 

15 

IF  <J  .EO.  JCIICM  GO  TO  22 

737 

J2  * J2  ♦ 1 

7*1 

JN  * J2  ♦ NC 

7*3 

*112.  J?)  • RllltJI 

7 5* 

AII2.JNI  * - RJ(I,J) 

765 

* ( 1 1 . J 2 ) • RJ(t.J) 

776 

AIII.JN)  • RKIfJ) 

100  7 

22 

CONTINUE 

101* 

SC«lE*l. 

1015 

NN  « 2*NC 

1017 

LL*LNEOF(50fNN, *,*,?,  SCALE, ERASE) 

1030 

PRINT  27, SCALE 

1035 

27 

FOR  HAT  I //,5X. *DF TERN  INANT*  ♦. 1PE12.*) 

1035 

12  * 0 

1036 

DO  26  I ■ l.N 
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i<m 

DO  17  IC  • 1.5 

1052 

17 

I F ( I .EU.  JC(  IC)  ) 

1050 

12  • 12  ♦ 1 

lOf  ’ 

II  • 12  ♦ NC 

1055 

DO  35  L ■ 1.5 

1055 

FE(I.l)  « 8(12.1) 

10b5 

15 

Fl<  l.LI  • tllliLI 

1100 

2o 

CONTINUf 

1103 

29 

RETURN 

1105 

END 

Table 
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4 


6 

ft 

ft 

ft 

ft 

ft 

7 

23 

42 

47 

93 

45 

96 

57 

60 

61 

63 

64 
70 

104 

123 


126 

116 

145 


146 

166 

200 

205 

207 

211 

213 

214 

215 

216 
224 
255 
306 
311 
326 

344 

361 


sob*  our  inf  phvscmio 

COMMON  *I12<29,25).R«56(2S.25)»  R3T(*S,»J>,  >Ow ( 2! . 25 t , Ff ( 25, / ) , 
1FI(29,6),  *1  (25. 251.  *1(29.251.  •<(»*,  41.  •1(29.4). 

2R4UI3.3.10).  R6AMI 3.3,lu> , Fid. 10).  0F6e*<3.l0).  MWB (25.9.10). 

3 0614(25.6. 10)  » *06(25.10) 

COM MON /ONE  /*(25)»V(25)»*B(25)»V9<?5)»AN,;(29),DE6<25),VV(i'5) 
1,FE1N(25).  FIIN(25).  RNuRM { 25.3).  JC(5) 

COMMON  /Two/  N.NNW,  NHAVE6,  ISTM,  IS*I»,  NC  . * 1 1 . GAMMA, M. T* . T • 

COMMON/THREE/  WAVE6 ( 10) • WNdOI.  *01(10), T6 

COMMON  /SEVEN/  ARE*.  9,  0.  *0E»  6Ee.  *TfT6E(3),  TlT6E(<t) 

RE*6  K 
DO  3 I • 1 

C«**4*M00E  6 • INCIDENT  ♦ D1FRACTE0  R0TEN7TA6* 

F £ 1 1,6)  • Ff  ( 1.4)  ♦ FE ( X . 5 ) 

3 F I ( I , ft)  • F 1 ( I > 4 1 ♦ Fid, 5) 

FACH  • ( ft4*2 ) /AREA 

FAC6  • FACMASORTIK) 

FACF  • F ACM  ♦ K 

00  1 6 • 1.3 

001  Ml  • 1,6 
RA  • 0.0 

RM  . 0.0 

I F ( Ml  .tO.  4)  GO  TO  1 
1F(M1  .EO.  5)  GO  TO  1 

(NTIGRATt  RRESSURt  COMPONENTS  OVER  ROOT 
00  5 I • l.N 

RM  • RM  ♦ Ff(I.Ml)  * RNORM ( 1,6)  • DF6IT) 

9 «*•**♦  FII(, Ml)  A RNORM (1,6)  • 0/6(1) 

IFIMl  ,GT.  41  SO  TO  * 

C*4***A00E0  MASS  4N0  OAM»(HG  IN  OIAECTION  6 OOE  TO  MOTION  Ml  AT 
(•••••MAVE6ENGTM  H. 

RM0I6.MI,  161  ■ RM*FACM 
R 6 AM ( 6 . Ml , III  • R4AFAC6 
GO  TO  1 

C ♦♦•••NAVE  FORCE  AMPLITUDE  AND  PHASE  IN  niRFCTION  mj  OUE  TO 
C****MNC  IOENT  4 A VE  AT  4AVE6ENGTH  16 
a F5I6,  161  • SORT( RM**2  ♦ RA**2)  • MAC* 

0E6 Fft( 6, 16  > • ATAN2 (-RA.RM) 

1 CONTINUE 

IN  • N ♦ 1 

IMA*  • N ♦ NNM 

DO  30  I ■ IN,  IMA* 

DO  6 6 • 1,4 

C*****COMPUTE  ROTENTI 46  AT  FREE  SURFACE  *OINTS  USING  6REENS  THEORUM 
DO  4 2 • 1.  N 

CO  10  IC  « 1,5 

10  I F ( J .EO.  JCdO)  GO  TJ  4 

FEd.6>  • FE  ( 1 , 6 ) ♦ FE(J,6)*RJ(I,J)  -eT(J.6)*RJ(I»J) 

Fid, 6)  • Fid, 6)  ♦ FE(1.6)*RJ(I.J)  ♦/ T ( J . 6 ) *R  I ( I . J I 

4 CONTINUE 

FE  ( 1, 6 ) • (Fed. 6)  - RK(I»6t)/TP 
6 Fid, 61  ■ (Fid, 6)  - *6  d , 6 ) ) /TP 
C«****MOOE  6 • INC  I Ot  NT  ♦ OIFRACTED  R0TFNTIA6* 

FFlI.b)  ■ Ffc ( I , 4 ) ♦ FE ( I , 5 ) 

Fid, 6)  ■ Fid. 4)  A Fid, 5) 
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375  II  • I - N 

C*****NQn  OIHENTIONALIZE  FREE  SURFACE  POSTTflN  «TTM  WAVELENGTH 
377  XOKII.il>  • X<I>  * 8/  WAVE  l ( II  > 

410  00  2 Ml  - 1,6 

C*«***WAVE  AMPLITUDE  ANO  PHASE  AT  POINT  TI  0U=  TO  NODE  Ml  AT  WAVFl FNGTH 
412  HUB (II.Nl.IL)  • SOR II F£ ( I.Nl ) **2  ♦ FT(T.Nl>**?>  * K 

443  2 DELWCII.Hl.IU  • ATAN2  ( -FI  ( I, HI  > , tf(I,«tlll 

470  30  CONTINUE 

472  7 RETURN 

473  END 
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2 


?. 

2 

2 

2 

2 

2 

2 

2 

2 

10 


10 

12 

32 

77 

180 

221 

223 

224 
241 


324 

326 

327 
334 
337 
353 
355 

436 

440 

441 
460 

534 

542 

544 

561 

632 

633 


SUBROUTINE  POTOUT 

CONNON  RI12(25.25)»RK56<25»25>»  POT  ( > 8,  ■»?  ) , H3w(?5»25)»FE(2r'»6)» 
1FK29.6),  R I ( 25, 25 ) » RJI25.25),  RMPB.4),  *L(25.4), 

2RNU(3. 3,101,  RLAN( 3.3.10) . FB(3.10I.  Oe L *3 ( 3, 10 ) « HWB < 2 5, 5 . 1 0 ) » 

3 DE  LW( 25, 6, 10 ) » X0K2S.10) 

CONNON  /TWO/  N.NNW.  NwAVEl.  ISYN,  I SK  T •«  N'T,  PIE, GANNA, M, TK, T® 

CONHON/THREE/  WAVEL (10),  WN(10),  BPL( IP! , II 

CONHON  /SEVEN/  AREA,  8,  0,  ROE.  GF®.  3TIUF(3),  TITLE(B) 

CONNON  / EIGHT/IBIHU(3»3»3),  LBLAN» 3, 3. 3 ) , IB L FB ( 3 , 3 ) . L »L HW« ( 7, 3 ) , 
ILBL ( 10, 3 ) , DE  G ( 3,10) 

1001  F0RNATI//3X,  3A10,  / <5X,  3A10,  10®n.4)) 

1002  FORNAT  ( / / 3 X.  3All),  / (5X,  3A10,  10*10.4  /5X,  3A10,  1JF10.4/I) 

1003  FORNAT ( 5X,  3A10.  10F10.4  / 5X,  3A1D.  10*10.4  /) 

10U4  FORNAT!  //3X,  3 Aid,  2 X,  1 OF  10. 4 / 3X,  3*10.  ?X,  10*10.4) 

PRINT  2000,  BTITIE 

2000  FOkNAT  (1H1,  POX,  ‘NONOINENSIONAL  POT*NTI*L  COEFFICIENTS*  III  ?*v 
1 * W • SORTIG/B),  W2  • G/B*  / 25X*B  • *,  3410  / 

1 25X*G  • ACCELERATION  OF  GRAVITY*/ 

1 25 X**OE  • NASS  OEnSITY  OF  FLMTP*/ 

22SX  *ETA  • INCIDENT  WAVE  ANPL I TUOE  */PBX,*WAVEL  • INCIDENT  OR  GENE 
3RATE0  WAVE  LENGTH*//) 

00  0 II  ■ 1.10 

9 DEGU.IL)  • 3URT(  (GEE*BOL(Il)  ) /(T®*BI) 

PRINT  1004,(LBL(1,K),K  . 1,3),  ( BOL ( IL)  » IL  • 1 . 1C), 

1 (LBK2.K),  K • 1,3),  (OEG(l.IL),  IL  • 1,10) 

PRINT  1001,  ( LBL  ( 3»K ) » K*  1,3),  ( ( ( LBLNIM  I,  J,K  I , K*  1,3), 

1CRNUII, J . IL  > » IL  • 1.19),  J « 1,3),  I • 1.3) 

PRINT  1C01  , ( LBL ( 4,K ) ,K  • 1,3),  ( I ( LB  LA" ( T , J , K ) , K • 1,3). 
llRLANII, J,IL),  IL*  1.10),  J • 1,3).  T • 1.3) 

DO  1 I • 1,3 
00  1 IL  ■ 1,10 

1 DEGU.IL)  • 57.  298  * OELFB  ( 1 , IL  > 

PRINT  1002,  ( LBL ( 5, K ) , K • 1,3),  ((l*Lc*(I»K)»K*  1,3), 

1 ( FB 1 1, IL  > . IL  « 1.10).  ( LBL ( 6, K ) , K • 1,  3 ) , ( D*G ( I , IL ) » IL • 1 . 1 0 ) 
2,1*  1,3) 

00  2 I ■ l.NNW 
00  8 IL  ■ 1,10 
DEGU.ILI  ■ O.C 

IF  I IL  .GT.  NWAVEL  ) GO  TO  B 
DEGU.IL)  *X0L(t,ILI*8/B0L(IL) 

B CONTINUE 

PRINT  1C02,  (LBLCB.K).  K • 1,3),  (LBL(9.K1.  K * 1.3), 

1 (XOL(I.IL).  IL  • 1,10),  (LBLUC.K),  «*l,3  ),  ( Of  G <1 , IL  ) , H *1 , 10 ) 
00  3 J • 1,3 
DO  3 IL  • 1,10 

3 OEGIJ.IL)  • 57.298  * DELWd.J.IL) 

PRINT  1003.  I (LBLHWBU.K),  K • 1,3),  I HWB  ( I . J . IL  ) . IL  ■ 1.10, 

1 (LBLI7.K),  K • 1,3),  (DEGU.IL).  IL  • 1.10).  J • 1,3) 

IF  (XOUI.l)  ,LT.  0.)  GO  TO  4 
DO  5 IL  ■ 1.19 

5 DEGll.ILI  • 57.298  * 0ELWU.6.IL) 

PRINT  1003,  (LBLHWBU.K),  K*l,3),  ( HWB(  1 , 6,  IL  I , IL  * 1,10), 

1 (LBL(6,K),K*1,3),  (DEGU.IL),  IL  • 1.1N) 

GO  TO  2 

4 DU  7 J ■ 1,3 
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RUNT  VERSION  FtB  74 


17*12  04/23/76 


635 

DO  7 II  . i,lu 

636 

7 

DEC  IJilLI  • >7.298  ♦ D£LV<  l»  J*3»  Ml 

i 

669 

PRINT  lbi.3, 1 llBLHNBI  J,K>,  K • 1,3). 
1 ( LBL 1 6, F ) , K>1,3),  ( 0 1 C ( J-3, IL  > , U 

INrfMT,J,  III,  IL 
■ 1*101,  J*  4,6) 

732 

2 

CONTINUE 

735 

RETURN 

735 

END 

Table  D-2.  Continued 
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'SW#*- 


1 * 1 0 > • 
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12 

12 

1<> 

15 

17 

20 

34 

36 

52 

55 

57 

60 

63 

65 

73 

75 

102 

113 

124 

125 
132 
143 

154 

155 
165 
167 
175 
205 
207 

230 

231 
252 
255 
257 

265 

275 

277 

315 

316 

317 

321 

322 

323 
325 
342 
362 


FUNCTION  INEQF«M,N,NI,A,8»0TRMNT,7» 

C..  SOLVES  SIMULTANEOUS  LINEAR  EQUATIONS  RV  GAUS  S T AN  reduction. 

C..  FORTRAN  TV  EQUIVALENT  OF  LNEOS. 

REAL  A<N,N> .8<N»M»  » Z< M ) . OTR NNT.RN A V . RNC XT , V. 00V 

NN1-N-1 

DO  40  J-1,N«1 

J1*J*1 

C..  FIND  ELEMENT  OF  COL  J.  ROWS  J-N,  WHTCM  MAS  MAX  ABSOLUTE  VALUF. 
LMAX ■ J 

R MA  X • A6  S ( AIJiJI  ) 

00  6 K-Jl.N 
RNEXT«ABS(A(K, J I ) 

IF  UMAX  ,6E.  RNEXT)  GO  TO  6 
RMAX-RNEXT 
LNAX-K 
e CONTINUE 

IF  (LMAX  .N6.  Jt  GO  TO  10 
C..  MAX  ELEMENT  IN  COLUMN  IS  ON  DIAGONAL 
IF  ( A ( J . J ) ) 20>V4.20 

Ct.  MAX  ELtMENT  IS  NOT  JN  DIAGONAL.  EXCMAVGF  ROUS  J ANO  LMAX. 

10  DO  12  L'J.N 
W*A<J,L> 

A( J,L)*A(LMAX,L> 

12  A(LHAX.L>*U 
DO  14  L-l.NI 
VB(J.L) 

B(J.L)*9(LMAX.U  — — L 

14  BtLMAX.D'V 

DTRMNT  • -OTRNNT 

C..  ZERO  COLUMN  J 8EL04  THE  DIAGONAL. 

20  ZUI'l./AIJ.JI 
00  30  K'Jl.N 
IF  ( A < K . J | ) 22.  30.22 
22  W»-Z ( J ) *A < K. J ) 

DO  24  L * J1 » N 

24  A(K.Lt*H*A( J.LI*4(K.L> 

OU  26  L * I . N 1 

26  6(K.U«W*R(J.L)»B(K.L) 

30  CONTINUE 
40  CONTINUE 

IF  (A(N.N) ) 42.74.42 
42  Z(N)«1./A(N.N) 

C..  OBTAIN  SOLUTION  8 Y BACK  SUBSTITUTION. 

00  50  L»1.N1 
50  B<N.L)«Z(N|*S(N,L> 

00  60  K ■ l ,NN I 

J«N-K 

Jl*J»l 

00  58  L • 1 » N1 

w»o. 

DO  56  I * J I » N 
56  6«Al J. I)*8( I.L*»W 
58  81 J.L)  ■<»U,L)-N)  *Z  < J) 

60  CONTINUE 
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WiA 


LNFOF 


RUNT  VERSION  FEB  7*  B 17*12  LR/23/76 


C..  EVALUATE  OETERHINANT. 


36* 

IF  (0TRNNT1  70.7**73 

366 

70 

00  72  J ■ 1 * N 

370 

72 

DTRHNT»DTRNNT*A(J*J) 

377 

7* 

L SE  OF • 1 

♦ Cl 

R f TURN 

C..  SINGULAR  NATRI*.  SET  ERROR 

A 01 

9R 

LNEOF  • 2 

603 

DTRUNT  *0. 

*06 

RETURN 

*03 

END 

Table  D-2.  Continued 
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RUNT  VERSION  F EH  74  B 17il2  04/23/76 


SUBROUTINE  CPV( X.Y.c.Cl.Sl.AD.AlO.Kt 
C ....CAUCHY  PRINCIPAL  VALUE  INTEGRAL. 


13  CONNOR  /TWO / N.NNW,  NVAVEl,  ISYH,  TS<TP.  NC.  P I E » GA NHA, N, TK , TP 

13  COHHON /SIX/XN(5)»CN<5) 

13  REAL  * 

li  IF  (Y  . Gt . 0.01  Y • -1.0E-09 

15  IT-ATAN2JY.X) 

24  TH-PIE/2.*TT 

C FOR  NEGATIVE  X, CORRECTION  TO  RANGE  OF  AT  AN?  . 

27  IF  C X.LT.O. )TH-TH»TP 

32  AA-K*Y 

34  E • E X P ( A A ) 

43  BB»K*X 

45  Cl-COSIBB) 

54  S1*SIN( 8B ) 

63  R-K*S0RTI X**2*Y**2) 

102  SUN1-0. 

103  SUN2-0. 

104  IF( R .GE . 10. I GO  TO  13 

107  SLNll-O. 

110  SUM22-0. 

111  FAC-1.0 

113  SUN  1C- 1 . 

114  SUN2C-1. 

115  SOLTh-0. 

116  CDLTM-O. 

117  ASSIGN  3 TO  LDC 

120  IF(X.E9.0.»ASSIGN  6 TO  LDC 

122  PL-1.0 

124  DO  1 L-..100 

125  OL-L 

126  FAC-FAC*OL 

130  RL-R*RL 

132  OLFAC-FAC*OL 

133  OLTH-Ol*TH 

135  Al-RL/OLFAC 

137  IF(ABSIC0LTH).IE.1.E>07)G0  TO  2 

151  SUN1C • ABS I A1/SUN1) 

157  I F ( SUN1C . LE. 1 . E -O > ) GO  TO  7 

165  2 CDLTH«COS<OLTH> 

171  SUN11-A1*C0LTM 

172  SUN1-SUN1+SUN11 

174  7 GO  TO  L0C.C3,») 

203  B SUN2C-0. 

204  GO  TO  5 

205  3 IF  I A t>  S ( S0LTH),L£.1.F— C7)60  TO  4 

217  SUN2C-ABSIA1/S0N2) 

225  IF  ( SUNJC  . LE  . I «F— vi5 ) GO  TO  5 

233  4 SOLTH- ilN(OLTH) 

237  SUN22- A 1*S0L  TH 

240  SUN2-SUN2-SUN22 

242  5 I F ( SUN  1C . Le . 1 . E-Oi. AND. SU«2C . LE , 1. E“05) GO  TO  fc 

261  1 CONTINUE 

263  6 C«GAMNAPAL0G(RI4SUHl 

DISCONTINUITY  OF  2PIS  IF  X NEGATIVE  TN  E T FUNCTION. 
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It 


RUNT  VERSION  Fee  7*  8 171X2  0*/23/76 


m 

300 
30  2 
30b 
313 

31* 

316 

321 

330 

333 

337 

3*0 

3*3 

3*7 

353 

35* 


1FIX.IT.C.)TH*TN-TP 

S*TH«SUM2 

AS-E*(Cl*OSl*S) 

*lo-t*(-Cl*S*$l*C> 

60  TO  9 

C...  IAGUERRE  OU»DRATORl-FIVE  POINT. 
13  00  1*  1*1.5 
A*XN( I ) *AA 

TERrt*CN(I)/(A*A*BB«88) 
SUNl*TERH*A»iUMl 
1*  SUM2*TERM*SU«12 
F»l. 

IF( X.LT.O. 1F*-1. 
A9«F*PIE*Sl*E-yJNl 
A1C— F*PIE*Cl*E*8B*SU*2 
9 RETURN 
END 
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RUM  VERSION  PER  74  8 1 7 « X » 04/23/76 


2 


2 

2 

2 


2 

2 


2 
2 
2 
2 
2 
3 J 


33 
6 7 


67 

71 

111 


156 

160 

161 

176 


261 

263 

264 
301 

362 

364 

365 
404 

460 

466 

470 

501 

516 

567 


SUBROUTINE  OTNOUT 

COMMON  *112(25, 25). *<56(25, 25),  P0T(25,»5).  HOW ( 2 5, 25 ) , Ff ( 2< , 6 ) . 
1EI<25,6>,  *1(25*25),  RJ(25*25>*  RK(?5,4),  *1(25,4)* 

2RMU (3*3*13) * RL1NI 3* 3*10) • FBI3.10).  Of L FR ( 3, 1C ) * HWR ( 2 5 , 6* 1 0 ) . 

3 OELWI 25*6* 1C)  » X01 (25*10) 

COMMON  /DO/  N.NNV,  SWIVEL.  ISYM*  ISKT*.  NO*  P 1 F , GAMMA , M,  Ttr , T P 
COMMON/  THRE  E / WAVELdO),  wN(10l#  BOldOl.U 

COMMON/FOUR/RAR(3,lJI,OELR(3.lO).  HW*(P<,3,n>.  L EL  W' ( 2 5 , 3 . 1 0 > . 

2 mwT(26.10>,DELWT(25.10)*RKHYD(3*3).  R<Mn*(3.3l.  *<TB(3.3). 

3 XG,  YG»  RMASS.  RINERT,  0AMP(3) 

COMMON  /SEVEN/  AREA.  8,  0,  ROE , GEE,  8TITIE<3),  TITLE  € f > 

COMMON  /NIKE/  lBLRAR  (3*3),  LBLHWR(5.3),  LHR(5*3) 

COMMON  / EIGHT/IBLNUI 3,3,3) . LRLA“(  3. 3. 3 ) , LB  L F8  ( 3 . 3 ) . L 8L  H WR  ( 7,  3 ) . 
1LBL ( 10. 3) • OEGI3.10) 

1001  F0RMATI//3X.  3*10,  / (5X,  3A10*  10F10.4)) 

1002  FORMAT  Unit  3 A10*  / ( 5 X , 3A10.  10*10.4  /5X.  3*10,  10F10.4/)) 

LOO 3 FORMAT!  SX*  3A10,  1JF10.4  / 5X,  3AtO*  13F10.4  /) 

1004  FORMAT!  / / 3 X , 3Alu.  2X,  10F10.4/  3A10,  ?Y,  1JF10.4) 

PRINT  2000 * AREA.  R.  XG,  YG.  RMASS*  R IN = R T. 0* M» (1 > , 0 AMP ( 2 ) . DA MP ( 3 ) 

2000  FORMAT!  IH  l*  20X*DYNAMIC  MOOEL  RESULTS*//*  ARF  A • ♦ F 10 . 3 . 5 X*R  ,*e  1 f> . 3 . 
1 5 X *XG**Fl0.3.5X*Yi«*F10.3,5X*MAS<:.*F1'),3.5Y*INERTIA**F10.3// 

3*  ADDITIONAL  DAMPING  ADDED-  IN  SWAY-*P6.2*  LAM0A11  IN  HF1VF- 

3*F6.2*  LAM0A22  IN  R0LL-*F6.2*  LAM0A33*  //) 

PRINT  2001* ( (RKHYO! I* J) * J»l. 3), 1-1.3) * ( ( RK MO* ( I * J ) * J ■ 1 » 3 ) . I • 1 . 3 ) 

2001  FURMATI*  SPRING  CONSTANTS  KT I <1>  <13  <21 

1 <22  <23  <31  <32  <33*/ 

2*  HYDRO  ST  ATIC*7X»9E10.3  / * MCO* ING*H X. RCi o. 3 / / / ) 

DO  9 1L  • 1*10 

9 CEG(l.IL)  ■ SORT! (GEE*BUL( ID ) /(T**R)) 

PRINT  D04,  IL8LU.M.K  • 1,3),  ( BOL  ( I L ) . Tl  • 1,101, 

1 ( L BL ( 2 *< ) * < • 1,3),  (DEG(l.IL),  TL  • ’.ID) 

DO  1 I • 1,3 
DO  1 IL  • 1.10 

l DEG(I.U)  ■ 57.298  « OeLR(I.IL) 

PRINT  1 002, ( L BLR ( 1,<  > . < • 1.3),((LBLR*»(T.<),  < • 1.3). 

1 IRAKI  I , I L ) , IL  • l.U).  ( L3L  ( 6,  <)  , < • 1.3),  ( DE  G ( I , I L ) . II  d . 1 0 ) , T 

2 • 1.3) 

DO  2 I ■ l.NNW 
DO  8 IL  • 1.10 

8 DEG(l.IL)  *XOL( I, IL )*B/BOL! ILi 

PRINT  lul2,  (LBl(3,<),  K * 1,31,  (LBL(9.<),  < • 1.3',. 

1 (XOL(I.U),  IL  • 1.10,  (LRLIIO.O.  <«1,3  {,  (OFC.  X.  ILI*  IL-1,10) 
CO  3 J ■ 1.3 
00  3 IL  • 1.10 

3 DEGU.IL)  • 57.298  *0E  L <R  ( I . J . II  ) 

PRINT  DO 3,  ( (LdLHWRI J,<).  < • 1,31,  ( HyR ( T , 2 , I L I , IL  • 1.10). 

I ( LBL  (6,0,  < • 1,3),  (OEG(J.IL).  IL  • 1,101.  J • 1,31 
IF  (XOLd.I)  .LT.  J.)  GO  TO  4 
DO  5 IL  • 1.10 

DEG  ( 2,  I L ) • 57.298  • OELWTd.IL) 

5 OtGd.lL)  • 57.298  * DELWII.6.IL) 

PRINT  1003.  (L8LHWB(7»<).  <*1,3I»  ( HWB( T » 4,  IL  I . IL  ■ I.1C). 

1 (LBL  ( 6.K  ) ,K»1,  i ) , (OEGd.IL),  IL  • 1.13) 

PRINT  1003,  (LBLMWR(4»<),<~ " >3).  ( HWT ( I , IL I . I L • 1.12). 

1 (L8K6.D,  < • 1.3),  ( OEG(  2,11),  TL  • 1.10) 


Table  D-2.  Continued 


133 


4m 


dtnout 


RUNT  VERSION  P EB  74  B 17«12  94/23/76 


637 

6 J TO  2 

640 

4 

00  7 IL  • l.n 

642 

DEGtl.lt)  • 57.298  * 0EIW<I,4,II) 

655 

7 

OfcG(Z.lL)  ■ 57.298  * DELWTt I» It ) 

670 

PRINT  10J3.  ( l 8LHWB ( 4 , K ) » K • 1,3).  ( H W« ( 1 . 4 , 1 L ) . I L • 1.10). 
1UBL16.K).  K • 1,3),  (0E6d.IL).  Tl.  • l.  I'))  # ( L BLHWP  ( 5 ,K  ) , *•  ) . 3) , 
2tHWTU.lL).  IL  ■ 1.10).(LBL(6,K).K  • 1,  3 ) . ( 0«=G  ( 2.  IL  ) . 1 1 ■ 1 . 1 C) 

1003 

2 

CONTINUE 

1010 

RETURN 

1010 

END 

Table  U-2.  Continued 
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^CALL  PHYSCL^ 
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POTENTIAL 
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7.  Program  Comments  and  Glossary  of  Terms. 


The  program  listing  contains  many  comments  which  aid  in  following 
the  logic  of  the  program.  Descriptions  of  variables  also  appear  where 
they  are  read  into  the  program . 

8.  Run  Time  and  Memory  Size. 

BRK2D  requires  about  70  seconds  of  central  processor  time  on  the 
CDC  6400  computer  to  compile  and  compute  results  for  10  different  beam 
wavelength  ratios.  A central  memory  of  about  55,000  octal  is  required. 

9.  Run  and  Card  Deck  Setup  Procedures  and  Special  Operation  Instruc- 
tions. 

In  order  to  run  the  FORTRAN  source  program  deck  on  the  University 
of  Washington  CDC  6400,  the  following  deck  is  required: 

BRK2D,CM55000,T100.  Job  card 

ACCOUNT  (Account  No.,  password) 

FORTRAN 

LGO(LC=6000)  LC  = line  count  value 

7/8/9 

FORTRAN  DECK 

7/8/9 

DATA  DECK 
6/7/S/9 

10.  Sample  Output  Data. 

Table  D-3  is  the  output  for  the  Oak  Harbor  breakwater.  The  input 
is  given  in  Table  D-l. 
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Table  D-3.  Example  output  for  program  BRK2D  (Oak  Harbor  breakwater) 


OilCNSI OMAL  JPPStT?  >*SeT?  **fO*niKT$  0>  SCW1S 


l 


#4 


ifi»'lfu'(M'<'0^'f'  X**' 
IM  ^ INi  IM  -f  N « l*M  fS,  l*\  fft. 


- ai  « rw»  o » N ^ <>  ft 

«KMf*ftrjfv.«*ft<\».X  M ft.  >»  ft 


3 _ - ■*'■'-*-*  — 

.■>*..*  r.  } N j ft  N ft. 

3 H Ift  ft  >.3  0 ft  ft  ft 


ft  X C O ft  **  't  Ift  i*l  O * £ ftif.  BTPIftXIft 

< •»!  O **ft  ft-  -•  ft  N C O x N O <*ft  ->  « r <• 

ftx^.  * *'■* 


*-»'**'  X r * **ft  -m  O 

1 1 r r : r r r • • 


■ K ftOI’xOMNCJ’yft  3 *SJ  **•  *"  •rft'  O C-  C O 
30  JOft*<*»  ftCxff^NXX^O  J O J 

*>«.*'*>.*•#***«»*’  <»i>re*i,»'»,**“«r  **•••*''“''* 

I I l I I I I I • I • 


.-.  O ■ J3ftftJift(lJOft«'.iftftwOi('C'lM, 
O *s-  tr  *■  ooft-ft-ftf-r^r.'  >o»',»*i*r*r-ft''^Oft'X'*3ic 
5«Nr'iftift|ft»"VH^ro--N<'  ftftift  «'  ***  *m  *■•  o 

»*  r •*  r »*  •*  i*  f r r r f •*•*  »*»"•*•  • 

c OCC  C c ■ * MM«>«  «i  » N fti  ft(  ® e c-  C O C C 
rrrrr.iftifft**|ift«»i(fftftif''«  c o f.  c. 

• l I l < • i l i l • • 


00>0300U*>33«'P,'C'>3  30  3^0*)  3 
O x -ft  ft  3 O <>  ’ft  O « ft  1 3 ft  »>  3 ift  ft  3 P I*  3 JftO 

0'Vftft»'>0»*l*’ftNN"'0WM>ftft»K0  3 N «ft  M O 

r<  V *1  ftft  Ml  •»>  *ft  »«J  -M  -• 

I I I I I « I I t • « I I I I t I I I I 

o o o a 

o O 300ftXftjftjM)'l)*rft^»iliXOft>OyO 

OCOOO‘ft'>x*'Vfto*s»ftft'iX'>r«wftftgftwftv<ir<x.  0 0 0 3 3 
Xft  Mft  <T>  MS  M>  .f  * **ft  Wt  «*ft  .r  # »ft  ftn  f»ft  •*  ^ ft  ift  1ft  A Ml 

t I I • I I I I I I I I 


N 3 IM  O >ft  N 

r*  5 -«  -n  .r  o 

H O O H IM  n 
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Table  D-3.  Continued 


OT*A*IC  *00*1  8ESUITS 


Si*  ST  AESUl  TAH  MEAVE/ETA  . 1 555  .**1?  .81*5  .1781  .168*  .09)6  .05*2  .023*  .6110  .CJ52 

HASi  «€l  TO  ETA  AT  *0  - 0E6  81.3023  -2*. 2**7  -»*.6?01  -177.808*  -1*2.**C3  -1*3.9711  -139. /1*2  -125.8*29  -107.80*2  -6*. 229* 


0t6  -11*. *9*6  -59. ?9*?  -ir»J. 05?9  -11?. *9*1  -117.11*9  -115. *06*  -111.5926  -U«.6C19  -46.66*9  -69.6965 
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APPENDIX  E 


DERIVATION  OF  PRESSURE  TO  SECOND  ORDER 
FOR  TWO  PROGRESSIVE  WAVES  AT  DIFFERENT  FREQUENCIES 


Consider  the  problem  of  the  nonlinear  interactions  of  waves  at  two 
distinct  frequencies  traveling  in  the  same  direction.  The  complete 
boundary  value  problem  is  well  known. 

The  Laplace  equation, 

= 0,  (E-l) 

applies  throughout  the  fluid  below  the  free  surface. 

The  boundary  condition, 

* 8 !y  * 2V<t>  ‘ v It  + \ V*  ' V (V<t,'V^  = °.  (E  — 2) 

3 1 

must  be  satisfied  on  the  free  surface,  y = n.  The  boundary  condition  on 
the  bottom  is  : 


lim  3$ 
y-*--°"  3y 


(E-3) 


for  an  infinitely  deep  fluid.  In  addition  a radiation  condition  re- 
quiring the  generated  waves  to  travel  away  from  the  body  is  needed  to 
ensure  uniqueness  of  the  solution. 


In  this  formulation  the  x axis  lies  in  the  direction  of  incident 
wave  propagation. 


The  difficulty  in  solving  this  boundary  value  problem  stems  from  the 
nonlinearity  of  the  free-surface  boundary  condition. 


In  order  to  "linearize"  the  free-surface  boundary  condition,  expand 
the  velocity  potential,  $,  in  a Taylor  series  about  the  undisturbed  free 
surface: 


$(x,n,t)  = <Kx,0,t) 


Also  expand  n 
n(x,t)  = 
4>(x,y,t) 


♦ n[Mx>Jr»tI]  + I n2  fWx.y.t), 

3y  2 1 1 3y  J A 

' y=0  1 y=0 

and  <F  in  power  series: 
cn^(x,t)  + e2n^(x,t)  ♦ 0(t3), 

* et(1)(x,y,t)  ♦ eV2)(x.X.t)  ♦ 0(e3). 


♦ 0(n3). 


(H-4) 


(E-5) 
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Substituuig  the  expansion  for  <{>  into  the  free-surface  boundary  con- 
dition : 

324>(1)  (x,y,t)  2 32t  (2)  3<|>(1J  2 3*(2) 

c ■ V7’  - + c -J— ^ ■ + ge  — + ge  -3- — 


nf  ,3<|>(1)  + 34>(1)  -f,  2,  S^2-1  -t  3<t-(2)  fAl 

+ 2[d^d3T  1 + -3T  • } + e { ~^x  1 + ly~  >}]m 


3 -t  3 , , 
[1  H + J 37j 


2 3*<2\ 

-rr  + £ -it-) 


1 r 3<j>  ^ 1 f 3<f>^  t,  2,  34>^2^  -t  34> 

+ 2 le  aT  1 + ^3T  J>  + « < “V  1 + "37“  J}] 


3 ■*  3 nfr  r 3 d>  ^ ^ + 3$  ^ ^ 

37  + J WPlK  1 + -17 


+,  2f34>(2)  f 34>(2)  tn 

-1}  + e {~37‘  1 37”  ^ } 


..(1)  . aA(l)  . ? ,a(2)  . a.(2) 

[ei^-  1 + ^ + E {il33r  1 + T>jJ  - 0(e3)  = 0 (E-6) 

on  y = n . 


Now  use  the  Taylor  expansion  for  <(>(x,n,t)  and  neglect  terms  of  or- 
der C3  in  the  boundary  condition: 

r32$(1)  (x,0,t)  (1),  ^ 33<fr(1\  . 2 324(2) 

e{ — z >—.?—• ’-2-  + cnv  (x,t)  — i— =— } + e — •E— j- 

3t  3y3t  3t 

f34>C1')  (1)  , , 32<t>(1\  2 3(J>(2:) 

♦ + Cx.t)  gy^—  + 8E  ^T 


„ 2r3*('1)  3$^  nr  3, 

+ 2e  i-V  -skr + “V  3t3T“]  + 0Ce  > = °- 


(E-7) 


Grouping  terms  by  order: 


First  Order  e : 

32*<»  a,0 

~7?~  * 8 ^ 


0 on  y = 0. 


(E-8) 


Second  Order  e : 

324> 3<J>(2)  (1)  3 f32<J>(li 

ir'*—  " * i?" 


-}  + 2 


a»»>  »2^‘) 
3x  3x3t 


Using  the  dynamic  boundary  condition  on  the  free  surface,  one  finds: 


n(x,t)  = - - i V4»*V4>}  on  y = n. 

g dt  £ 

Substituting  the  expansions  into  this  equation  yields: 
en{1)(x,t)  + eV2)  + 0(e3)  = - ~ + \ 


y=0 


i 37  {H  * 1 


y=0 


Substituting  for  the  right-hand  side  becomes: 

- 


= - i (eli- 

g 3t 


(2)  2 ..(1)  2 -.(1)  2 
♦ e <■  T ic— Z5J— > * (— rH  ]1 


9y 


2 (1)  .2(1)  . 

e n { > + 0(e3),  on  y = 0. 


g ' 3y9t 
First  Order  e: 
n 


(Dfx  = . 1 

g at 


(E-10) 


(E-ll) 


(E-12) 


Second  Order  e : 

( 2) , n . I ,(1)  9* 

n (x,t)  g 'n  3y3t 


CD  3AC1)  , 

12! } - - — { r 

at  2g  ax 


,3i{i 


CD  2 

■)  ♦ (“ 


a<j> 


CD 

17 


-)  ) 


or 


on  y = 0 

n(2)  f n _ £ ( i. 

n (x,t)  g { g 3t  3y3t 

3*U)  2 

♦ (— 3y~ ) > on  y = 0. 


1 , 1 a<t>(1)  aV1' 


.(2) 


(D 


♦ liLi).  i_ 

at  ' 2g  1 ^ ax  ; 


(E-13) 


Using  the  first-order  relationship  above  in  the  second-order  boun- 
dary condition  on  the  free  surface  (E-9)  , one  finds: 


aV2> 


at 
- 2 


g 


a# 


(2) 


k(D 


ay 


= *f  — 


1 a<t.v*'  a ,a*» 


— {• 


2.(D 


k(D 


g 9t  3y 


a<ji  (1^  aV1!  2 3* 


ax 


3x3t 


t1) 


at' 


9y 


-} 


3y  3yat 


CE— 14) 
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1.  First-Order  Solution  of  the  Boundary  Value  Problem. 


This  solution  results  from  the  superposition  of  the  velocity  poten- 
tials for  individual  waves: 

j(x,y,t)  = ^-=-  e ly  cos  (kjX  - u^t  + 6^ 
gA 

+ ^-=-  eK2y  cos  (k2x  - w2t  + &2)  . (E- 15) 

Check  the  solution: 


vVJ>  = 0. 


Lim  djf 
y + -<  ~ 


(1) 


3y 


32<t>^  34> 

-JL~r  + g — 

3t 


0 because  of  exponential  function. 


(1) 


k y 

gp-  = ~ g“1A1e  V cos  (kxx  - u^t  + 6^ 


k y 

- g“2A2e  2 COS  ^k2X  " W2t  + 62^ 


k.t 


+ g {wjAje  1 cos(kjX  - u>  t + 6j) 


k y 

+ ^>2^2e  27  cosC^2X  " * ^3  ) = 


Therefore,  this  is  a solution. 


Surface  elevation  then  becomes: 

.CD, 


n^(x.t)  = - ^ — t')  = - Aj  sin(kxx  - u^t  + 


A2  sin(k2x  - u2t  ♦ 62) . 


(E-16) 


To  prepare  for  the  second-order  solution,  construct  the  right-hand 
side  of  the  free-surface  boundary  condition  (E- 14) : 


(1) 


rl  3<f>^  3 , 3 2(J) ^ 1 ^ 3d. 

[2  ~3T  37  {-^r  * g % 


-}  - 2 


34C«  3V» 


3x  3x3t 


, tV'K  i 

- 2 -5T  ^y3F“1y=0  = 7 {gAl  sin(klx  " “l*  + 61} 


gA2  sin(k2x  - w 2 1 62)  HO)  - 2{  - u^Aj  sinfkjX  - u^t  + (E-17) 
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3K 


...  v >vr. *-> 


- uj2A7  sin(k2x  - u>2t  + 62)}x{u)1  Aj  cos  (kjX  - 
2 

+ to  2A2  cos(k0x  - u)2t  + 6-,)  ) - 2{u)jAj  cos  (k^x  - Wjt  + 6^) 

2 

+ ^2^2  cos(K->x  " w2t  + 52^X^0i>1  A1  s^n(kix  " “i1  + 

2 

♦ u A2  sin(k2x  = u>2t  + 52)  ) = 0. 

Since  this  condition  is  homogeneous,  the  first-order  potential  is 
the  solution  to  the  second-order  problem. 


2.  Second-Order  Results. 

The  free- surface  elevation  will  be  modified  when  terms  of  second 
order  are  included: 

rn  i .J1)  i ;u(1)  2 ad,(1)  2 

n IX, tj  ^2  t 3t  3y3t  2g  u 3x  ‘ 1 3y  J y=0 


= + (gA^  sin(k^x  - w^t  + 6^)  + gA2  sin(k2x  - w2t  + 50) } * 

g 

2 2 
{A^  sin(kjX  - u^t  + 6^)  + A2  u>2  sin(k2x  - + 62) } 

- ~ {[-  uijAj  sinCk^x  - w^t  ♦ 6^)  - W2A2  sxn^2X  ' <1)2t  + ^2^ 

2 

+ [w^Aj  cos  (k^x  - w^t  + 6^)  ♦ «2A2  cosC^2X  " W2t  + ^2^  ^ 


or 


gn^(x.t)  = u^Aj2  sin2(k1x  - uijt  ♦ 6^ 


♦ uijA^  sxn(k2x  - u)2t  + 62)  sinCk^x  - Ujt  + <5^ 
2 

+ U)2  A^2  sin(kjX  - Wjt  ♦ 6j)  sin(k2x  - u>2t  + 62) 


+ uj72A92  sin2(k2x  - i»2t  + 69) 


2 r\ 2 3 ah  2 2 

j {w^2A^2  sin2(kjX  - uijt  + 6 ^ ) 

♦ 2w1w2A1A2  sinfkjX  - Wjt  + fi^sinCkjX  - u>2t  + 62) 
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2 2 2 

+ u) 2 A2  sin  " u)2t  + ®2^ 

♦ (d^Aj"  cos2(kjXj  - uijt  + 6^) 

+ 2w1u2A1A2  cosCk^x  - uj^t  + 61)cos(k2x  - w2t  + 62) 

2 2 2 

+ (i>2  l<2  cos  C^2X  " w 2 t * 62^’ 

Using  the  trigonometric  relationships: 

2 2 2 2 2 2 2 
gn  (x,t)  = uij  A^  sin  (k^x  - w^t  ♦ 6^)  + w2  A2  sin  (k2x  - u> 2 1 + 62) 

+ ^ wj"aiA2  ^c°s[(k1  - k2)x  - (u^  - w2)t  + - 62] 

- cos[(kj  + k2)x  - (u>j  + w2)t  + 61  + fi^-} 

+ j u22a1A2  ^cosUki  ' k2)x  ’ ^W1  " “2^  + 61  " 62^ 

- cos[(k^  + k2)x  - (w1  + ui2)t  + 6^  + 62J> 

1 2 2 2 2 

- 2 Aj  + u 2 A2  ) - “jw2  A1A2  cost^ki  “ k2^  x 

- (tuj  - w2)t+  6^  - 62J. 

Combining  further: 

gn(2)(x,t)  = - y u12A12  cos[2{k1x  - u^t  + 6^] 

- i .22A22  cos[2{k2x  - w2t  + 62}]  (E-18) 

- (“j2  ♦ “)22^A1A2  C0S  ^kl  + k2^X  " ^“l  + “2^  + 61  + ^2^ 

12  2 

♦ j (“1  - 2wiw2  + °J2  ^A1A2  cost(ki  ‘ k2^x  ‘ (wl  * a)2^t  + 61  ‘ <52^» 

which  is  the  final  form  for  the  second-order  term  for  free-surface  ele- 
vation. 
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Now,  turn  to  the  equation  for  pressure  which  is  necessary  to  com- 
pute the  force  on  the  body. 


Take  the  pressure  to  be  zero  at  the  free  surface.  Then  Bernoulli's 
equation  may  be  written: 


P = - p — - - pV$ • v$  - pgy. 


(E- 19) 


Substituting  the  expansion  for 

P - r *(1)  , ^ , 1 f 2 difi^K2 

p - - p{c  — r ♦ e ♦ 2 ( jx"1 

♦ e2  C^-y-)  ] ♦ gy}  ♦ 0(e3). 


(2) 

Since  <j>  = 0,  we  can  drop  this  term  and  proceed  to  separate  the 

equation  by  order: 

p(1)  = - P-4J-  - pgy  (K - 20) 


at 


and 


(2)  ,.U)  2 zJV  2 


(E-21) 


Substituting  the  velocity  potential  into  the  equation,  one  finds: 
P^1"*  = - Pg(A1ekly  sin(k^x  - Wjt  + 6j) 


k y 

+ A^e  2 7 sin(k7x  - w2t  + 6?)  + y} 


(E-22) 


for  the  first  order,  and 

^ oijA^e^l^  sin(k^x  - m^t  ♦ 6j) 

k y 2 

- “2A2e  2 sin(^2x  " u>2t  * 52^J 

k y 

♦ [w^e  l7  cosCk^x  - ui^t  + 6^) 

+ w2A2ek2y  cos(k2x  - u>2t  + 6 2)  ] 2 } 


for  the  second  order.  Note  that  this  is  identical  to  part  of  the 
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equation  for  surface  elevation.  The  second-order  pressure  may  be  reduced 
to: 

P(2)  = - J {w12A12e2kly  + o)22A22e2k2' 

- 2u)^u)2AjA2e^kl  k2^ycos[(k^  - k2)x  - (u)^  -u^)t  + 6^  - 62J)  (E-23) 

which  indicates  that  the  second-order  pressure  is  composed  of  a component 
independent  of  time  and  at  the  "difference  frequency". 

This  is  surprising  since  the  equation  for  the  free-surface  elevation 
(eq.  18)  includes  terms  at  twice  the  incident  wave  frequencies  and  at  the 
sum  of  these  two  frequencies.  Using  trigonometric  relationships  the 
first  two  terms  in  equation  (E-23)  could  be  expanded  to  yield  terms  at 
twice  the  incident  wave  frequency.  A term  at  the  sum  of  the  two  inci- 
dent wave  frequencies  may  appear  in  the  pressure  computed  using  the 
velocity  potentials  representing  wave  diffraction  or  forced  oscillation. 
It  might  also  appear  if  the  present  analysis  were  carried  to  the  third 
order.  The  derivation  included  here  was  intended  to  reveal  the  presence 
of  a low-frequency  component  in  the  exciting  force  and  has  not  been  used 
to  determine  the  other  velocity  potentials  or  carried  beyond  the  second 
order . 


List  of  Special  Symbols  for  Appendix  E. 


^i  > a2 


= Wave  amplitudes 


kp  ^2 


= Acceleration  of 


gravity 


u>l  W2 

= Wave  numbers,  — , respectively 

= Cartesian  coordinates  (x-directed  parallel  to  the 
direction  of  wave  propagation,  y-directed  vertically 
upward) 


= Wave  phase  angles 


n(x,t) 


= Free-surface  elevation 


<Hx,y,t) 


Velocity  potential 


Vw2 


Wave  circular  frequencies 


APPENDIX  F 


PHYSICAL  PROPERTIES  OF  SEVERAL  FLOATING  BREAKWATERS 

Proposed  Oak  Harbor  Floating  Breakwater  (Davidson,  1971). 

a.  Physical  Properties. 

m = mass  per  unit  length  = 25.1  slug/ft 

2 

I = mass  moment  of  inertia  = 621  slug-ft  /ft 

x = x-coordinate  of  center  of  gravity  = 0.0  ft. 

^ (on  centerline) 

y^  = y-coordinate  of  center  of  gravity  = -2.34  ft  (below  WL) 

KH22  =64.5  lb/ft/ ft 

KH33  = 1,165  ft- lb/ ft 

All  other  KH . . = 0 
il 

b.  Mooring  Line  Tension  Response  (change  per  unit  displacement). 

77  = 1,170  lb/ ft 
AX 

— = 281  lb/ft 

Ay 

|I.  i,710  lb 

c.  Computed  Mooring  Spring  Constants  (depth  = 29.5  feet) 


KMn  = 

119  lb/ft/ft 

™12  - 

-5.24  lb/ft/ft 

m!3  - 

166  lb/ft 

mn  - 

-5.73  lb/ft/ft 

m22  - 

10.2  lb/ ft/ ft 

km23  = 

-3.37  lb/ft 

km3i  = 

160  lb/ ft 

**32  - 

2.06  lb/ft 
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KM,,  = 282.  ft-lb/ft 
3a 

Rectangular  Breakwater  Tested  by  Nece  and  Richey  (1972). 

Physical  Properties  (at  prototype  scale) . The  cross  section  is  a 
"rectanglcTof  beam  10  feet  and  draft  5 feet . 

m =.100  slugs/ft 

I = 2,740  slug- ft2/ ft 

x = 0.0  ft  (on  centerline) 

g 

y = -1 .0  ft  (below  WL) 

g 

KH22  = 640  lb/ft/ft 
KH33  = 5,340  ft-lb/ft 
All  other  KH^  = 0 


All  KM^ 


0. 


Rectangular  Breakwater  Tested  by  Sutko  and  Haden  (1974) . 

Physical  Properties  of  Model.  The  cross  section  is  a rectangle 
of  beam  0.333  feet  and  draft  0.222  feet. 

m = 0.143  slug/ft 

I = 0.023  s lug- ft2/ ft 

x = 0.0  ft  (on  centerline) 
g 

yg  = -0.123  ft  (below  WL) 

KH22  = 20.7  lb/ft/ft 
KH33  = 0.244  ft-lb/ft 
All  other  K1L  ^ = 0 
All  KM^  = 0 

Alaska- Type  Breakwater, 
a.  PhysicaL  Properties, 
m =62.3  slug/ft 
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I = 4,234 slug-ft/ft 

x = 0.0  ft 

& 

y = -1.3  ft  (below  WL) 

S 

KH2,  = 528  lb/ft/ft 
KH„  = 32,885 ft-lb/ft 
All  other  KH. . = 0 

b.  Mooring  Line  Tension  Response  (change  per  unit  displacement). 

= 97.0  Ib/ft 

Ax 

— = 90.5  lb/ ft 

Ay 

fr  * -572  lb 

c.  Computed  Mooring  Spring  Constants  (tide  = +7.0  feet). 

KMU  = 3.0  lb/ft/ft 

KM12  = 0.245  lb/ft/ft 
KM13=-9.23  lb/ft 
KM21  = 0.302  lb/ft/ft 
KM22  = 1.91  lb/ft/ft 
KM23  =-2.68  lb/ft 
KM31  =-9.52  lb/ft 
KMj2  =-2.82  lb/ft 
KM33  = 88.9  ft-lb/ft 

5.  Friday  Harbor  Breakwater, 
a.  Physical  Properties, 
m = 61.02  slugs/ft 
I 3 4,160  slugs-ft3/ft 
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x = 0.0  ft  (on  centerline) 

g 

y =*'0.49  ft  (below  WL) 

KH21  = 884  lb/ft/ft 

KH53  = 55,610  ft- lb/ft 

All  other  KH. . = 0 
ij 

b.  Mooring  Line  Tension  Response. 

= 222  lb/ ft 
Ax 

~ = 25.0  lb/ft 
Ay 

H ■ 657  ib 

c.  Computed  Mooring  Spring  Constants  (tide  = +5.33  feet). 


KMn 

= 

6.46 

lb/ft/ft 

mi2 

= 

0.510  lb/ft/ ft 

KMl3 

= 

18.5 

ft- lb/ft/ft 

km2i 

= 

0.510  lb/ft/ft 

KM22 

* 

0.390  lb/ft/ ft 

km23 

X 

1.71 

ft- lb/ft/ ft 

km3i 

■ 

18.6 

lb/ft 

^32 

X 

1.71 

lb/ft 

^33 

X 

64.6 

ft- lb/ft 

159 


APPENDIX  G 


DATA  SUMMARY  SHEETS  FOR  FRIDAY  HARBOR 
FLOATING  BREAKWATER  (WINTER  1975) 


Appendix  G contains  a summary  of  all  the  data  recorded  at  the  Fri- 
day Harbor  breakwater  during  the  winter  season  of  1975.  Seven  tapes  were 
recorded  during  this  period,  with  a total  of  95  records.  The  tapes  are 
numbered  in  sequence  from  FH7-1  through  FH13-8.  The  date  of  each  tape 
is  given  along  with  the  pertinent  statistical  data  for  each  record  in 
the  tapes.  The  number  of  days  and  hours  given  for  each  record  begins 
with  the  day  and  hour  given  for  that  particular  tape. 

All  minimum  and  maximum  values  are  measured  from  zero  mean.  The 
transmitted  wave  data  were  digitally  high-pass  filtered  (cutoff  fre- 
quency was  0.05  hertz)  before  these  calculations  to  remove  tidal  draft. 


SUMMARY  OF  STATISTICAL  DATA  FOR  FRIDAY  HARROR  FLOATING  BREAKWATER  (FH7  - 1330  - IZ/30/7AI 

(MAI.  AMO  MIN.  VAIJES  MEASURED  FROM  ZERO  MEANI 
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APPENDIX  H 

INCIDENT  AND  TRANSMITTED  WAVE  SPECTRAL  PLOTS 


Appendix  H contains  the  incident  and  transmitted  wave  spectral 
plots  along  with  the  corresponding  transmission  response  curve  for  11 
representative  records.  The  data  for  the  first  nine  were  recorded  at 
Friday  Harbor,  Washington,  during  the  winter  of  1975.  Figures  H-ll  and 
il-12  were  computed  from  similar  data  collected  in  Alaska  during  the 
winters  of  1974  and  1975. 

The  original  time  series  were  high-pass  filtered  at  a cutoff  fre- 
quency of  0.05  hertz  to  remove  tidal  drift.  Each  series  consisted  of 
2,048  data  samples  and  were  sampled  at  a period  of  0.5  second  for  the 
Friday  Harbor  data  and  0.44  second  for  the  Alaska  data. 

The  standard  deviations  and  corresponding  overall  transmission 
coefficients  for  each  of  the  Friday  Harbor  plots  are  given  in  Appendix 
G. 


FREQUENCY  (Hz) 

Incident  and  transmitted  wave  spectral 


FREQUENCY  (Hz) 

Incident  and  transmitted  wave  s 
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Figure  H-4.  Incident  and  transmitted  wave  spectral  data  (FH7-9) . 
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Figure  H-8.  Incident  and  transmitted  wave  spectral  data  (FH11-11). 


FREQUENCY  (Hz) 

Figure  H-9.  incident  and  transmitted  wave  spectral 


INCIDENT  WAVE  SPECTRA 


Figure  H-10.  Incident  and  transmitted  wave  spectral  data  (TK7-1),  Tenakee  Springs,  Alaska. 
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Figure  H-U.  Incident  and  transmitted  wave  spectral  data  (SK4-10),  Sitka,  Alaska. 


APPENDIX  I 


LOW-FREQUENCY  SPECTRAL  ANALYSIS  OF  FORCE  DATA 


Appendix  I contains  the  low-frequency  autospectral  and  cross-spec 
tral  plots  for  record  FH7-8.  The  data  were  recorded  at  Friday  Harbor, 
Washington,  on  6 January  1975  at  0030  hours. 

The  original  time  series  were  low-pass  filtered  at  a cutoff  fre- 
quency of  0.2  hertz  and  every  eighth  data  point  used  to  generate  a new 
time  series.  This  gives  256  points  with  a sampling  period  of  4 seconds 


FREQUENCY  (Hz) 

. Coherency  between  northeast  and  southeast  forces 


FREQUENCY  (Hz) 

Figure  1-6.  Phase  between  the  southeast  and  southwest  forces. 
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Figure  1-7.  Coherency  between  the  southeast  and  southwest  forces. 


APPENDIX  J 


HIGH-FREQUENCY  SPECTRAL  ANALYSIS  OF  FORCE  AND  MOTION  DATA 


Appendix  J contains  the  incident  wave  spectral  plot  along  with  the 
autospectral  and  cross-spectral  plots  for  the  force  and  motion  data  for 
record  FH7-8.  The  data  was  recorded  at  Friday  Harbor,  Washington,  on  6 
January  1975  at  0030  hours. 

The  incident  wave  spectra  was  unfiltered.  All  the  force  and  motion 
spectral  data  were  digitally  high-pass  filtered  at  a cutoff  frequency  of 
0.1  hertz.  The  autospectral  data  is  plotted  as  a percent  of  the  vari- 
ance, i.e.,the  total  area  under  the  spectra.  Wave  heights,  forces,  and 
motions  were  measured  in  feet,  pounds,  and  feet  per  second  square,  re- 
spectively. 

All  spectra  were  computed  from  2,048  data  points  sampled  at  0.5- 
second  intervals. 


Figure  J-2.  boutheast  force  spectra  (FH7-8) . 


Figure  J-3.  Northeast  force  spectra  (FH7-8). 


J-4.  Southeast  force  spectra 
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Figure  J-7.  Phase  between  the  southeast  and  southwest  forces. 


Figure  J-9.  Heave  motion  spectra 


Figure  J-10.  Sway  motion  spectra  (FH7-8) . 


Figure  J - 1 1 . Roll  motion  spectra 


Figure  J-I2.  Phase  between  heave  and  roll 
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Figure  J-15.  Coherency  between  sway  and  roll. 


Figure  J-17.  Coherency  between  sway  and  heave. 


APPENDIX  K 
WAVE  MEASUREMENT 

1 . Wave  Staff  Design. 

A block  diagram  of  the  wave  staff  and  associated  electronic  cir- 
cuits is  shown  below: 


The  wave  staff  itself  consists  of  a length  of  PVC  tubing  which  is 
spirally  wound  with  a resistance  wire,  such  that  when  it  is  immersed 
in  seawater,  the  electrical  resistance  varies  in  direct  proportion  to 
the  length  of  the  exposed  staff. 

The  electronic  circuits  driving  the  wave  staff  consist  of  a fixed 
frequency  square  wave  oscillator  (having  a precisely  controlled  output 
amplitude)  driving  a precision  bilateral  current  source  with  an  output 
current  directly  proportional  to  the  input  voltage.  Thus,  the  wave 
staff  is  driven  by  a current  source  of  constant  magnitude,  but  one  which 
changes  direction  with  each  one-half  cycle  of  the  square  wave  oscilla- 
tor. The  output  of  the  wave  staff  then  is  a square  wave  voltage  with  a 
magnitude  (peak  to  peak)  that  is  directly  proportional  to  the  length  of 
the  exposed  wave  staff.  This  output  is  fed  to  a high  input  impedance 
voltage  follower  circuit  which  serves  as  a buffer  between  the  wave  staff 
and  the  ac  detector  circuit.  The  precision  ac  detector  circuit  uses 
two  operational  amplifiers  in  conjunction  with  two  diodes  to  form  a 
precision  full-wave  rectifier  circuit  that  is  capable  of  operating  at 
very  low  input  voltages.  Ordinary  diode  detector  circuits  cannot 
operate  on  ac  signals  of  peak  magnitude  less  than  the  forward  voltage 
drop  of  the  diodes  and  produce  large  conversion  errors  unless  the  sig- 
nal magnitude  is  large  with  respect  to  the  diode  voltage  drop.  A gain 
control  has  been  incorporated  in  the  detector  circuit  so  that  full-scale 
output  can  be  set  at  any  positive  value  up  to  +10  volts  with  a wave 
staff  resistance  of  300  ohms  up  to  3,000  ohms. 

Alternating  current  is  used  to  drive  the  wave  staff  to  avoid  both 
the  corrosion  effects  that  would  occur  if  direct  current  were  used  and 
the  dc  offset  which  occurs  as  a result  of  the  use  of  dissimilar  metals 
in  a conducting  solution.  The  latter  is  eliminated  by  use  of  ac  coup- 
ling in  the  output  from  the  wave  staff. 

Bench  tests  of  the  wave  staff  electronic  circuits  were  made  using 
a 1,000-ohm  variable  precision  resistor  in  place  of  the  wave  staff.  The 
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circuit  was  adjusted  to  produce  an  output  range  of  0 to  10  volts  with 
the  resistor  varied  from  0 to  1,000  ohms.  Linearity  was  determined  to 
be  0.1  percent  of  full  scale  over  this  range. 

Tests  were  also  made  to  determine  the  effect  of  temperature  on 
sensitivity  and  zero  drift.  A decrease  in  sensitivity  was  noted  with 
decreasing  temperature  of  about  0.03  percent  of  reading  per  °Celsius  over 
the  temperature  range  of  0 to  24°Celsius.  A zero  drift  of  2 millivolts 
was  also  noted  over  the  same  temperature  range.  A +10  percent  change  in 
supply  voltage  from  the  nominal  +^15  volts  produced  no  observable  change 
in  output.  If  we  assume  an  operating  temperature  range  of  +5°Celsius, 
the  maximum  error  in  the  wave  staff  electronics  due  to  the  combined  ef- 
fects of  nonlinearity  and  sensitivity  variations  with  temperature  is 
+0.2  percent  of  reading.  Since  the  primary  interest  is  in  a dynamic  mea- 
surement of  waves,  the  zero  drift  noted  will  have  negligible  effect  on 
the  experiment  since  temperature  variations  of  any  appreciable  magnitude 
will  only  occur  over  long  periods  of  time  compared  to  the  wave  periods. 

Further  calibration  tests  were  conducted  using  actual  wave  staffs 
of  1-inch  diameter  and  20-foot  lengths,  and  3.5-inch  diameter  and  8-foot 
lengths  at  various  depths  of  immersion  in  saltwater.  These  tests  were 
conducted  from  a dock  at  Shilshole  Bay  on  Puget  Sound.  Because  of 
ripples  and  waves  on  the  water  of  the  order  of  1 inch  (peak-to-val ley) 
it  was  difficult  to  obtain  a highly  precise  measurement.  The  output 
was  recorded  on  a strip  chart  recorder  and  it  was  therefore  possible  to 
average  these  variations  to  some  degree.  The  readout  resolution  of  the 
strip  chart  (and  accuracy)  is  about  +^/4  of  a minor  division.  Full 
scale  across  the  chart  is  50  minor  divisions  and,  thus,  the  resolution 
is  about  0.5  percent  of  full  scale.  Some  nonlinearity  is  noted  near  full 
immersion  (see  calibration  curve).  Some  offset  was  expected  because  of 
the  finite  resistance  of  the  saltwater  path  in  the  ground  return  which 
is  not  taken  into  account  during  initial  calibration  of  the  wave  staff 
unit.  The  initial  calibration  is  made  with  the  wave  staff  on  the  dock 
where  full  scale  and  zero  are  set  by  making  actual  contact  between  the 
ground  wire  and  the  wave  staff  resistance  element  at  the  corresponding 
ends.  However,  measurements  were  made  of  the  resistance  of  the  salt- 
water path  to  ground  in  the  same  location  where  the  wave  staffs  were 
immersed  and  the  value  of  resistance  measured  (on  the  order  of  10  ohms) 
does  not  account  for  the  offset  observed  at  full  immersion.  In  addition, 
the  offset  should  occur  at  all  readings  and  it  does  not.  Therefore,  it 
is  believed  that  the  nonlinearity  observed  is  a result  of  some  other 
phenomenon  as  yet  undetermined.  Both  units  produced  highest  accuracy 
near  center  scale  with  decreasing  accuracy  toward  either  end.  Overall 
accuracy  including  end  points  is  about  ^3  percent.  If  the  rq^ge  of  oper- 
ation is  reduced  so  as  not  to  use  the  last  1 foot  on  each  end  of  the 
wave  staff,  the  accuracy  is  improved  to  about  +^1  percent. 

The  output  from  the  wave  staff  electronic  circuit  is  fed  directly 
into  a voltage  to  frequency  converter;  the  frequency  output  is  then 
counted  and  stored  on  separate  storage  registers,  once  every  50  milli- 
seconds. If  an  8-bit  register  is  used  for  the  wave  staff  measurement. 
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the  maximum  count  that  can  be  stored  is  255;  therefore,  the  sample  time 
must  be  on  the  order  of  25.5  milliseconds  (maximum  count  divided  by  maxi- 
mum frequency  output  from  voltage  to  frequency  converter) . The  wave 
buoys  use  an  8-bit  register  with  a 32.5-millisecond  sample  time 
while  the  wave  staffs  use  a 16-bit  register  with  a 250-millisecond  sam- 
ple time. 

The  error  due  to  gain  instability  and  nonlinearity  of  the  voltage 
to  frequency  converter  is  of  such  low  magnitude  that  it  can  be  neglected 
and  the  overall  accuracy  of  the  recording  is  essentially  the  same  as 
given  for  the  wave  staff  unit  by  itself  (i.e.,  between  +1  and  <5  percent 
depending  on  the  range  of  operation  on  wave  staff). 


2.  Spar  Buoy  Design. 

Spar  buoys  were  used  at  two  of  the  sites  because  of  their  advantage 
in  handling  and  transport  and  because  they  minimized  the  placement  diffi- 
culties due  to  navigational  hazards,  water  depth,  and  tidal  conditions. 

The  spar  buoys  were  made  of  two  PVC  pipes  coupled  together  near  the  cen- 
ter of  the  buoy.  The  lower  section  is  a 15  foot  by  6 inch  pipe  filled 
with  styrofoam.  The  top  section  is  12  feet  by  3 inches  wherein  the  upper 
8 feet  is  wound  with  a resistance  wire  which  measures  wave  elevation. 

The  wave  staff  electronics  are  mounted  inside  the  top  section,  above  the 
waterline,  with  the  remainder  being  filled  with  a wood  core  to  add  stiff- 
ness. The  buoys  also  have  a 2.5-foot  diameter  damping  plate  mounted 
on  the  bottom  and  are  anchored  using  a dual  point  mooring  system  with 
the  anchor  lines  attached  at  the  center  of  drag  on  the  buoy  to  prevent 
it  from  being  pulled  underwater  in  strong  currents.  One  of  these  buoys 
was  tested  in  the  Puget  Sound  just  north  of  Seattle.  Its  performance 
exceeded  expectations  both  in  terms  of  minimized  response  to  the  waves 
and  accuracy  of  wave  height  measurement.  Figure  K-l  gives  a sample  of 
the  output  from  the  buoy's  wave  staff  in  saltwater  for  a plus  and  minus 
1 foot  excitation  of  the  buoy  in  heave.  This  was  accomplished  by  push- 
ing the  buoy  up  and  down  by  hand.  Some  distortion  results  from  this 
approach  which  shows  up  in  the  output  of  the  accelerometer  mounted  at 
the  center  of  the  response  of  the  buoy  in  heave  and  roll  in  calm  water. 
The  natural  periods  for  heave  and  roll  taken  from  these  plots  are  approxi- 
mately 18  and  14  seconds,  respectively.  These  are  well  out  of  the  range 
of  the  3-to-4-second  wave  periods  expected  at  the  site.  Visual  obser- 
vations of  the  buoy  in  waves  in  excess  of  1.5  feet  indicated  no  ob- 
servable heave  or  roll  motion,  but  some  yaw  about  the  anchor  line  caused 
by  the  current  and  wind.  This  motion  resulted  in  less  than  a 1 foot 
variation  from  the  buoy's  horizontal  position  in  calm  water  and  appeared 
to  have  periods  in  excess  of  30  to  60  seconds.  For  comparative  mea- 
surements, the  buoy  was  located  about  30  feet  from  an  existing  four-gage 
array  of  1-inch  diameter  Oceanographic  Services,  Inc.  resistance  wire 
wave  staffs.  A comparison  of  simultaneous  output  from  the  two  wave 
staffs  (buoy  mounted  and  stationary)  is  shown  in  Figure  K-4.  The  auto- 
spectras  computed  from  data  obtained  from  one  of  the  stationary  wave 
staffs  and  from  the  spar  buoy,  in  a 25-miles  per  hour  storm  with 
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*0.4  cm/acala  div. 
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Figure  K-l.  Wave  and  acceleration  data  for  par  uoy. 


maximum  wave  heights  in  excess  of  i.5  feet  are  shown  in  Figure  K-5. 
These  spectra  were  computed  from  simultaneous  records  of  20  minutes  in 
length. 


! 


I 


I 


NORMALIZED  SMOOTHED  SPECTRAL  ESTIMATE  (Hz) 

0.76  1.53  2.29  3.06  3.82  4. 
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Figure  K-5.  Wave  spectra  from  spar  buoy  and  stationary  staff. 
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